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Abstract 



This thesis describes the development of some basic mathematical tools of 
wide relevance to mathematical physics. Transmission and reflection coeffi- 
cients are associated with quantum tunneling phenomena, while Bogoliubov 
coefficients are associated with the mathematically related problem of exci- 
tations of a parametric oscillator. While many approximation techniques for 
these quantities are known, very little is known about rigorous upper and 
lower bounds. 

In this thesis four separate problems relating to rigorous bounds on trans- 
mission, reflection and Bogoliubov coefficients are considered, divided into 
four separate themes: 

• Bounding the Bogoliubov coefficients; 

• Bounding the greybody factors for Schwarzschild black holes; 

• Transformation probabilities and the Miller-Good transformation; 

• Analytic bounds on transmission probabilities. 
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Preface 



This thesis looks at a number of problems related to the derivation of rig- 
orous bounds on transmission, reflection, and Bogoliubov coefficients: To 
set the stage, we shall first briefly describe the general ideas underlying the 
Schrodinger equation, and the concept of the WKB approximation for barrier 
penetration probability. In addition, we shall present a discussion of some 
general features of scattering theory in one space dimension. By considering 
one-dimensional problems involving an incident beam of particles, we shall 
derive an important connection between reflection and transmission ampli- 
tudes. Furthermore, we shall collect several known analytic results, and show 
how they relate to the general results presented in this thesis. We shall also 
review and concisely describe the concept of quasinormal modes, and see how 
most of the concepts introduced here are important tools for comparing the 
bounds derived in the body of the thesis with known analytic results. 

The technical heart of the thesis is this: We shall rewrite the second-order 
Schrodinger equation as a set of two coupled first-order linear differential 
equations (for which bounds can relatively easily be established). Systems of 
differential equations of this type are often referred to as Shabat-Zakharov 
systems or Zhakarov-Shabat systems. After this initial investigation, we 
shall use this system of ODEs to derive our first bound, and then continue by 
finding several slightly different ways of recasting the Schrodinger equation 
as a Ist-order Shabat-Zakharov system, in this way deriving a number of 
slightly different bounds. 

Regarding the chapter "Bounding the Bogoliubov coefficients", we have 
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developed a distinct method for deriving general bounds on the Bogoliubov 
coefficients, providing a largely independent derivation of the key results; a 
seperate derivation that short-circuits much of the technical discussion. 

Proceeding further along this branch of our investigation, we shall con- 
sider the Regge-Wheeler equation for excitations of a scalar field defined on 
a Schwarzschild spacetime, and adapt the general analysis of the previous 
chapters to this specific case. We shall demonstrate that rigorous and ex- 
plicit analytic bounds are indeed achievable. While these bounds may not 
answer all the physical questions one might legitimately wish to ask, they 
are definitely a solid step in the right direction. 

We shall then use the Miller-Good transformation (which maps an initial 
Schrodinger equation to a final Schrodinger equation for a different potential) 
to significantly generalize the previous bound. Moreover, we shall then use 
the Miller-Good transformation to generalize the bound to make it more 
efficient. 

Finally we shall consider analytic bounds on the transmission probabilities 
obtained by comparing a simple "known" potential with a more complicated 
"unknown" one. In this case we shall obtain yet another (distinct) Shabat- 
Zakharov system and use it to (partially and formally) "solve" the scattering 
problem. In this case we can derive both upper and lower bounds on the 
transmission coefficients and related Bogoliubov coefficients. 

Chapter by chapter outhne 

This thesis is divided into twelve main chapters. The first chapter is devoted 
to describing the general ideas of the Schrddinger equation and the concepts 
of WKB approximation for barrier penetration probability. Furthermore, we 
introduce the concept of the classical turning point, which is one of the key 
ideas in developing the WKB estimate. Moreover, these general concepts are 
important for understanding the bounds we will derive on transmission and 
refiection in Bogoliubov coefficients. 



iv 



In chapter [2| we shall introduce scattering theory in one space dimen- 
sion. This is an elegant topic that is mathematically simple and physically 
transparent. We shall apply the Schrodinger equation to a generic system 
to identify the potential-energy function. Furthermore, we shall derive a 
significant relationship between reflection and transmission amplitudes by 
considering one-dimensional problems with an incident beam of particles. 

In chapter |3} we shall concentrate our attention on collecting several 
known analytic results, and show how they relate to the general results 
presented in this thesis. We shall review and briefly describe the concept 
of quasinormal modes, and see how most of the concepts introduced here 
are important for comparing the bounds we shall derive with known ana- 
lytic results. By taking specific cases of these bounds and related results it 
is possible to reproduce many analytically known results, such as those for 
the delta-function potential, double-delta-function potential, square poten- 
tial barrier, tanh potential, sech^ potential, asymmetric square- well potential, 
the Poeschl-Teller potential and its variants, and finally the general Eckart- 
Rosen-Morse-Poeschl-Teller potential. 

The next two chapters, chapter |4] and chapter [5} can be seen as two 
deeply interconnected chapters. The key idea in chapter |4] is to recast the 
Schrodinger equation as a Ist-order Shabat-Zakharov system. In chapter [5} 
we shall use the Shabat-Zakharov system of ODEs to derive our first bound 
on the transmission, reflection, and Bogoliubov coefficients. 

In chapter [6} we shall deal with some specific cases of these bounds and 
develop a number of interesting specializations. We shall collect together 
a large number of results that otherwise appear quite unrelated, including 
reflection above and below the barrier. In addition, we have divided the 
special case bounds we consider into five special cases: special cases 1 — 4, and 
"future directions" . At the end of this chapter, we take further specific cases 
of these bounds and related result to reproduce many analytically known 
results. 

In chapter [7| we shall re-cast and represent these bounds in terms of 
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the mathematical structure of parametric oscillations. This time-dependent 
problem is closely related to the spatial properties of the time-independent 
Schrodinger equation. 

In chapterjs} we shall re-assess the general bounds on the Bogoliubov coef- 
ficients developed in [SH], providing a new and largely independent derivation 
of the key results, one that short-circuits much of the technical discussion 
in [HH]. 

In chapter |9| we shall develop a complementary set of results — we shall 
derive several rigorous analytic bounds that can be placed on the greybody 
factors. Furthermore, we shall consider the greybody factors in black hole 
physics, which modify the naive Planckian spectrum that is predicted for 
Hawking radiation when working in the limit of geometrical optics. 



In chapter 10 , we shall use the Miller-Good transformation (which maps 
an initial Schrodinger equation to a final Schrodinger equation for a different 
potential) to significantly generalize the previous bound. At the end of this 
chapter, we shall discuss the possibility of using the Miller-Good transforma- 
tion to derive generalized special-case bounds to make them more efficient. 



In chapter 11, we shall develop a new set of techniques that are more 
amenable to the development of both upper and lower bounds. Moreover, we 
shall derive significantly different results (a number of rigorous bounds on 
transmission probabilities for one dimensional scattering problems), of both 
theoretical and practical interest. 



In chapter 12 , we finally conclude with a brief discussion of lessons learned 
from these rigorous bounds on transmission, refiection, and Bogoliubov co- 
efficients. 



Structure of the thesis 

This thesis has been written with the goal of being accessible to people with a 
basic background in non-relativistic quantum physics, especially in transmis- 
sion, refiection, and Bogoliubov coefficients. Mathematically, the key feature 
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is an analytic study of the properties of second-order linear differential equa- 
tions, and the derivation of analytic bounds on the growth of solutions of 
these equations. 

This thesis is made up of twelve chapters and five appendices. Four of the 
appendices are papers published or submitted on work relating to this thesis. 
All of them were produced in collaboration with my supervisor, Professor 
Matt Visser. At the time of writing three papers have been published [891 
[501 m\ , and the latest has been submitted for refereeing [22] • 

Use of references 

Regarding referencing — For completely non controversial background infor- 
mation (and only for completely noncontroversial items) we will often just 
reference Wikipedia or similar reasonably definitive web resources. For more 
technical information, especially recent research, we will always directly cite 
the appropriate scientific literature. 
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Chapter 1 

General introduction 



1.1 Introduction 

This chapter is an introduction to the topic of developing rigorous bounds on 
transmission, reflection, and Bogoliubov coefficients. We shall introduce the 
basic ideas underlying the Schrodinger equation, and its application to the 
wave-function that describes the wavelike properties of a subatomic system. 

We shall also review the concept of the WKB approximation, which is 
an important and significant method to derive approximate solutions for the 
wave function. For instance, as we shall show, the WKB approach can be 
used as a "basis" for formally writing down the exact solutions. Most physi- 
cists, and many mathematicians, have seen how important the WKB approx- 
imation is for estimating barrier penetration probability. Unfortunately, the 
WKB approximation is an example of an uncontrolled approximation, and 
we do not know if the resulting estimate is high or low. As part of the main 
work reported in this thesis, we modify, improve, and extend the approach 
originally developed by Visser [88J . 

We shall derive a number of rigourous bounds on transmission proba- 
bilities (and reflection probabilities, and Bogoliubov coefficients) for one- 
dimensional scattering problems. The derivation of these bounds generally 
proceeds by rewriting the Schrodinger equation in terms of some equivalent 
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system of first-order equations, and then analytically bounding the growth 
of certain quantities related to the net flux of particles as one sweeps across 
the potential. 

While over the last century or more considerable effort has been put into 
the problem of flnding approximate solutions for wave equations in general, 
and quantum mechanical problems in particular, it appears that as yet rela- 
tively little work seems to have been put into the complementary problem of 
establishing rigorous bounds on the exact solutions. We have in mind either 
bounds on parametric amplification and the related quantum phenomenon 
of particle production (as encoded in the Bogoliubov coefficients), or bounds 
on transmission and refiection coefficients. 

In this thesis, we introduce and prove several rigorous bounds on the 
Bogoliubov coefficients associated with a time-dependent potential, and also 
derive several rigorous analytic bounds that can be placed on barrier trans- 
mission probabilities. As a specific application, we shall then explore grey- 
body factors in black hole physics, which modify the naive Planckian spec- 
trum that is predicted for Hawking radiation when working in the limit of 
geometrical optics. 

Additionally, we will extend these ideas to address topics of considerable 
general interest in quantum physics, such as transmission through a poten- 
tial barrier, and the related issue of particle production from a parametric 
resonance. This is an example of finding new physics (and new mathematics) 
in an old and apparently well-understood area. 

To begin with, we need to briefiy describe the concept of the Schrddinger 
equation, and the rationale behind the WKB estimate for barrier penetration 
probability, otherwise the rigorous bounds on transmission, refiection, and 
Bogoliubov coefficients will be difficult to understand. 
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1 . 2. THE SCHRODINGER EQ UATION 



1.2 The Schrodinger Equation 

The Schrodinger equation was discovered by the Austrian physicist Erwin 
Schrodinger in 1925, it describes the space - and time - dependence of the 
quantum amphtude that characterizes quantum mechanical systems [T] . 

Both Erwin Schrodinger and Werner Heisenberg independently developed 
different versions of the "modern" quantum theory. Schrodinger's method 
relates to partial differential equations, whereas Heisenberg's method uses 
infinite-dimensional matrices. 

However, both methods were soon shown to be mathematically equiv- 
alent. Furthermore, from the modern viewpoint it seems very clear that 
Schrodinger's equation has a clearer physical interpretation via the classical 
wave equation. Indeed, the Schrodinger equation can be shown to be a form 
of the wave equation applied to matter waves |2]. 

It is apparent that this equation defines the behaviour of the wave function 
that describes the wavelike properties of a subatomic system. Furthermore, it 
deals with the kinetic energy and potential energy, both of which contribute 
to the total energy. It is solved to derive the different energy levels of the 
system. More generally, Schrodinger applied the equation to the hydrogen 
atom, and its properties can be predicted with remarkable precision. It 
should be remarked that the equation is applied widely in atomic, nuclear, 
and solid-state physics [5]. 

Actually there are two slightly different equations which go by Schro- 
dinger's name as follows: 

1.2.1 The time- independent Schrodinger equation 

We start with the one-dimensional classical wave equation |2|, 



d'^u 1 d'^u 




Let us consider the separation of variables 
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which then leads to 

dx^ dt 



m'^ = -,^p(-y-^- (1.2.3) 



When we introduce one of the standard wave equation solutions f{t) such as 
exp{iut), we easily obtain 

^-^^•<^)' (^-) 

It is now easy to "derive" (in the sense of a physicist's plausibility argument) 
an ordinary differential equation describing the spatial amplitude of the mat- 
ter wave as a function of position. We note that the energy of a particle is 
the sum of kinetic and potential parts 

2 

E = ^ + V{x), (1.2.5) 
m 

which can be solved for the momentum, p, to obtain 

p = {2m[E -V{x)]y/^. (1.2.6) 

We now see that it is convenient to use the de Broglie formula to get an 
expression for the (position dependent) wavelength 

, h h 



If we recall uj = ^tiu and v\ = v, then the term uj'^/v'^ in equation (1.2.4) 
can be rewritten in terms of A: 



2 



47r2z/2 47r2 2m[E - V{x)] 



y2 y2 ^2 



1.2.8) 



Additionally, when this result is substituted into equation (1.2.4), we also 
"derive" the well-known time-independent Schrodinger equation, 

^ + |^[E-F(x)Mx) = 0. (1.2.9) 
Let us now rewrite the above equation in a more standardized form, we get 



2m dx 



2 



+ V{x)^Ij{x) = Ei/j{x). (1.2.10) 
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We now have all the important information about our system. Moreover, 
this single-particle one-dimensional equation can clearly be extended to the 
case of three dimensions, where it becomes 

V'^'^ir) + V(r)ilj(r) = Eiph). (1.2.11) 

2m 

A two-body problem can also be treated by this equation if the mass m is 
replaced by the reduced mass /x: 

I ^ I ■ (1-2-12) 

mi ~'~ 1712 

Nevertheless, it is important to point out that this analogy with the 
classical wave equation only goes so far. We cannot, for example, "derive" 
the time-dependent Schrodinger equation in an analogous fashion, at least 
not without several additional hypotheses. (For instance, the time-dependent 
Schrodinger equation involves the partial first derivative with respect to time 
instead of the partial second derivative.) Finally, we would like to comment 
that historically, Schrodinger presented his time- independent equation first, 
and then went back and postulated the more general time-dependent equa- 
tion [2]. 



1.2.2 The time- dependent Schrodinger Equation 

In this section we now present the time-dependent version of the Schrodinger 
equation. Although we were able to "derive" the single-particle time-indepen- 
dent Schrodinger equation starting from the classical wave equation and the 
de Broglie relation, the time-dependent Schrodinger equation cannot be "de- 
rived" using elementary methods, and is generally given as a postulate of 
quantum mechanics [2]. 

In other words, we shall postulate the single-particle three-dimensional 
time-dependent Schrodinger equation as 

if&^ = -i^vV(r,t) + V{r)^{r,t). (1.2.13) 
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We now focus on the case where V is assumed to be a real function, which 
represents the potential energy of the system. It is very easy to see that 
the time- dependent equation can be used to derive the time-independent 
equation. If we write the wavefunction as a product of spatial and temporal 
terms, ^jir^t) = ip{r) f(t), then equation (1.2.13) becomes 



il){r)ih 



d/(t) 
dt 



fit) 



or 



1 



ih df{t) _ 



2m 



+ V{r) 



■ip{r). 



[1.2.U) 



;i.2.15) 



It is easy to see that the left and right hand sides must each equal a constant 
E when the left-hand side is a function of t only and the right hand side is a 
function of r only. (This is just the usual separation of variables technique.) 

Alternatively, if we appropriately denote this separation constant by E, 
(since the right-hand side clearly must have the dimensions of energy), then 
we extract two ordinary differential equations, specifically 



1 d/(t) 
/(t) dt 



iE 



and 



2m 



V^V^(r) + V{r)ip{r) = Eip^r). 



[1.2.16) 



;i.2.17) 



The equation (1.2.17) is once again the time- independent Schrodinger equa- 



tion. Furthermore equation (1.2.16) is easily solved to yield 



f{t) = exp{-iEt/h). 



1.2.18) 



Most generally, we can show that the Hamiltonian in equation (1.2.17) is a 



Hermitian operator, and that the eigenvalues of a Hermitian operator must be 
real, so E is real. This implies that the solutions /(t) are purely oscillatory, 
since f{t) never changes in magnitude [recall Euler's formula exp{±i6) = 
cos 6 ±i sin 6] . In the following if we set 



i'irjt) = ■ip{r) exp{—iEt/h), 



;i.2.19) 
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then the total wave function ip{r,t) differs from iplr) only by a phase factor 
of constant magnitude. 

We can easily show that the quantity |'?/'(r, t)p is time independent as 
follows: 

|^(r,t)|2 =ilj*{r,t)ilj{r,t) = exp{iEt/h)'^*{r)exp{-iEt/h)-^{r) = V^*(r)V^(r). 

(1.2.20) 

Furthermore, if ^^{r^t) satisfies (1.2.19), then the expectation value for any 
time-independent operator is also time-independent. Thus it is easy to see 
that 



{A)= / ^*(r,t)i^(r,t) 



il)*{r)A%l){r) 



^1.2.21) 



Wave functions of the form (1.2.19) are called stationary states. The state 
'?/'(r, t) is "stationary" , but the particle it describes is not. It is now easy to see 
that equation (1.2.19) represents a particular solution to equation (1.2.13). 
The general solution to equation (1.2.13) will be a linear combination of these 
particular solutions [2] 



ip{r,t) = Cj exp{-iEit/h)4)i 



;i.2.22) 



In the next section, we shall introduce an important technique, the WKB 
approximation, which will be used several times in the body of this thesis. 



1.3 WKB approximation 

The WKB (Wentzel-Kramers-Brillouin) approximation is also known as the 
WKB J (Wentzel-Kramers-Brillouin- Jeffreys) approximation, or sometimes 
the JWKB approximation. The basic idea is it estimates a real Schrodinger 
wave function by a sinusoidal vibration whose phase is presented by the 
space integral of the classical momentum, the phase integral, and whose 
amplitude varies inversely as the fourth root of the classical momentum. In 
fact, in its original 1800's incarnation as the Jeffreys approximation, the WKB 
approximation was already a meaningful expression for the physical waves of 
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optics, acoustics, and hydrodynamics. After 1925, this approximation was 
rapidly apphed to the new Schrodinger probability waves [9j. 

The WKB approximation is an important method to derive approximate 
solutions and estimates for many physical problems. For instance, it is mainly 
applicable to problems of wave propagation in which the frequency of the 
wave is very high, or equivalently, the wavelength of the wave is very short 
(compared to the typical distance over which the potential varies). Despite 
the fact that the WKB solutions are approximate solutions, sometimes they 
are amazingly accurate [TO] . 

Let us begin with the one-dimensional time- independent Schrodinger equa- 
tion [g 

- + V{x)^{x) = E^{x), (1.3.1) 

which can be rewritten as 

d^ _ 2m 



where d'^ip{x)/dx'^ = second derivative with respect to x, iplx) = Schrodinger 
wave function, E = energy and V = potential energy. 

We can now write the wavefunction in terms of the exponential function 
by putting it in the form [U] 

^{x) = A{x) exp{iS{x)/h). (1.3.3) 



Substituting ^(x) into equation (1.3.2), we derive 

A{x) S'{xf - ihA{x) S"{x) - 2ihA'{x) S'{x) - A\x) = 2m{E- V) A{x). 

(1.3.4) 

By comparing the first two terms, we expect that the quasi-classical region 
is given by 

S'ixf ^hS"{x). (1.3.5) 



We take the real and imaginary parts of equation (1.3.4): 

S'{xf = 2m{E-V) + h^A"{x)/A{x), (1.3.6) 
-S"{x) = 2S'{x){d\n{A'{x))/dx). (1.3.7) 
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Now we are considering only one- dimensional problems, which of course also 
include radial motion in central potentials. One can then express equation 



(1.3.7) in the form 



d /I , d^(x) , , , , 



and one finds 



A(x) 



C 



:i.3.8) 



;i.3.9) 



In equation (1.3.6), let us neglect the term h? A" {x) / A{x) compared to S"(x)^. 
(This is where the approximation is made.) The resulting equation 



S\xY = 2m{E -V{x)), 



can then easily be integrated: 



S{x) = ± dx' ^/2m{E -V{x')). 



;i.3.io) 



:i.3.11) 



Substituting (1.3.9) and (1.3.11) into (1.3.3), one finds 



^(x) = V exp j ±i [ dxp{x)/h\ (1.3.12) 

± Vp(x) I J } 



with momentum 



p{x) = ^j2m{E - V{x)) 
Now we introduce the notation 



k{x^ 



2 2m\E - V{x)\ 



;i.3.i3) 



;i.3.i4) 



By the JWKB approximation, we derive 



^ _ ^ exp[z /A;(x)] ^ ^ exp[-i / k{x)] 



;i.3.15) 



This shows that the JWKB approximation is a fruitful method of calculation, 
that can be used to develop a perturbation theory. 
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1.4 Classical turning points 

We can start by considering one of the most interesting aspects regarding 
the WKB approximation] that being what happens at the classical turning 
points where V{x) = E. In fact it is easy to realize that as long as we 
keep away from these points, the approximation works very well indeed. To 
get near or pass through a turning point one has to go beyond the WKB 
approximation. The most straightforward way to do so is by using a linear 
approximation to the Taylor series expansion of the potential in the vicinity 
of the classical turning point. The exact solution to this approximate problem 
is given in terms of an Airy function. {Bessel function of order |.) Using this, 
the standard approach is now to derive a specific way of patching the wave 
functions on either side of the turning point — this leads to the so-called 
"connection conditions". Finally, it is interesting to note that historically 
the WKB approach to barrier penetration application very quickly yielded 
significant achievements in terms of understanding alpha decay lifetimes [3]. 

1.5 Discussion 

In this chapter, we introduced the Schrodinger equation, which is a specific 
partial differential equation used in the development of the "new" (1925) 
quantum theory. The Schrodinger equation was discovered by the Austrian 
physicist Erwin Schrodinger in 1925, and describes the space -and time- de- 
pendence of quantum mechanical systems pL]. In addition, physicists quickly 
applied the WKB approximation to the new Schrodinger probability waves. 
The WKB approximation is generally applicable to problems of wave prop- 
agation in which the frequency of the wave is very high, or equivalently, the 
wavelength of the wave is very short. 

The problem of finding approximate solutions for wave equations in gen- 
eral, and quantum mechanical problems in particular, has been extensively 
considered over the last century or two. However, it appears that as yet rel- 
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atively little work seems to have been put into the complementary problem 
of establishing rigourous bounds on the exact solutions. 

As the theory of the WKB approximation, and the concept of the time- 
independent Schrodinger equation, both underlie all our subsequent analyses, 
we have presented a very general introduction to these concepts first — so 
that the bounds we will soon derive on transmission, reflection, and Bogoli- 
ubov coefficients will be easier to understand. 

Finally we believe that this introduction has provided sufficient context 
for the reader to appreciate the role played by the various topics to be dis- 
cussed in this thesis, and to place them into a wider perspective. In brief, 
quantum mechanics is a generic tool for addressing empirical reality, and in 
this thesis we are probing the complementary problem of establishing rigor- 
ous bounds on the exact solutions. 
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Chapter 2 

Scattering problems 



2.1 Introduction 

In this chapter we shall present quantum scattering theory in one space 
dimension. It is a beautiful subject that is mathematically simple and phys- 
ically transparent. Moreover, it still contains various important results [17]. 

One-dimensional scattering problems appear in a vast variety of physical 
contexts. For instance, in acoustics one might be interested in the propaga- 
tion of sounds waves down a long pipe, while in electromagnetism one might 
be interested in the physics of wave-guides. Another important context which 
we want to stress in this chapter is that in quantum physics the canonical 
examples related to one-dimensional scattering theory are barrier penetra- 
tion and reflection. In contrast, in classical physics an equivalent problem is 
the analysis of parametric resonances . 

Furthermore, when considering the basic ideas of "reflection and trans- 
mission probabilities" , we shall introduce a useful technique to derive a con- 
nection between reflection and transmission coefficients, showing that they 
are related via a conceptually simple formalism. This technique will be used 
several times in the main part of this thesis. 

In particular, at the end of this chapter we shall (purely as an example) 
illustrate how to derive either transmitted or reflected probability waves as 
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a result of scattering of an object in the delta-potential well. More generally, 
we are specifically interested in the Schrodinger equation as shown below 



in equation (2.2.1) in conditions where the potential V{x) is zero outside 
of a finite interval — mathematically we are most interested in considering 
potentials of compact support. (Though much of what we will have to say 
will also apply to potentials with suitably rapid falloff properties as one moves 
to spatial infinity.) 

2.2 Reflection and Transmission Probabilities 

Let us consider the one-dimensional time-independent Schrodinger equa- 
tion ISZ]-I5I] 

h'^ d^ 

ipix) + V{x) ipix) = E ipix). (2.2.1) 



2m dx^ 

If the potential asymptotes to a constant, 

V{x ±oo) ^ V±oo, (2.2.2) 

then in each of the two asymptotic regions there are two independent solu- 
tions to the Schrodinger equation 

I r ^ \ exp{±ik±^x) 

i/j±{x ±oo) ^ j== . (2.2.3) 

Here the ± distinguishes right-moving modes e"'"*'^^' from left-moving modes 
^-ikx^ while the ±oo specifies which of the asymptotic regions we are in. 
Furthermore 

To even begin to set up a scattering problem the minimum requirements are 
that potential asymptote to some constant, and this assumption will be made 
henceforth. The so-called Jost solutions are exact solutions J^±{x) of the 
Schrodinger equation that satisfy 
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^+(^^+00)^"+ 7^ +/?+ 7^ , (2.2.6) 

V "^+00 V "^+00 

and 

^ +00) ^ "^P(-!^°°^) ^ (2.2.7) 

J_(a; ^ -00) — > q;_ h P- , . (2.2.8) 

Y "^-00 Y "^-00 

Identifying the reflection and transmission coefficients. 
There are unfortunately at least four distinct sets of conventions in common 
use, depending on whether or not one absorbs factors of a/a5±oo into r and 
t respectively, and on whether one chooses to focus on left-moving or right- 
moving waves as being primary. Let us, for the current section, adopt the 
convention of not absorbing the factors of A/A^ioo into r and t. (We shall 
discuss the other convention a little later in this chapter). We start by 
introducing a minor variant of Messiah's notation [51] 



—00) — > t+ exp{+ik^oox) , (2.2.9) 

J+{x +00) ^ exp{+ik+^x) + ^^ ^^P(-^°o x)^ (2.2.10) 

Y ^+00 

By comparing these two different forms for the asymptotic form of the Jost 
function we see that in this situation the ratios of the amplitudes are given 
by 

+ - ^+ =t+:l:r+. (2.2.11) 



Thus we obtain 



We also derive (in this set of conventions) 



r+ = = — • (2-2.12) 



00 «+ «H 



1 



''+00 



(2.2.13) 

□ 
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Thus we have demonstrated that a+ and the (right-moving) Bogoh- 
ubov coefficients, are related to the (left-moving) reflection and transmission 
amplitudes by 

= ^; U = -. (2.2.14) 

Without further calculation we can also deduce 

The explicit occurrence of /c+oo and /c_oo is an annoyance, which is why many 
authors adopt the alternative normalization we shall discuss later on in this 
chapter. 

In Bogoliubov language these conventions correspond to an incoming flux 
of right-moving particles (incident from the left) being amplified to amplitude 

at a cost of a backfiow of amplitude In scattering language one should 
consider the complex conjugate J'^ — this is equivalent to an incoming flux of 
left-moving particles (incident from the right) of amplitude being partially 
transmitted (amplitude unity) and partially scattered (amplitude P*^). If the 
potential has even parity, then the left-moving Bogoliubov coefficients are 
just the complex conjugates of the right-moving coefficients, however if the 
potential is asymmetric a more subtle analysis is called for. 

The second interesting issue is that we can deal exclusively with a+ and 
dropping the suffix for brevity — if information about a_ and is 
desired simply work with the reflected potential V{—x). It should also be 
borne in mind that the phases of (3 and (3* are physically meaningless in that 
they can be arbitrarily changed simply by moving the origin of coordinates (or 
equivalently, physically moving the location of the potential). The phases of a 
and a* on the other hand do contain real and signiflcant physical information. 

For completely arbitrary potentials, with no parity restriction (so the 
potential is neither even nor odd), a Wronskian analysis yields (see for ex- 
ample [ni], noting that an overall minus sign between Messiah and the con- 




(2.2.15) 
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ventions above neatly cancels): 

A;_oo[l-|r+|'] = A;+oo|t+r; (2.2.16) 

A;_oo|^-r = A;+oo[l-|r-|']; (2.2.17) 

k.oot- =k+^t+; (2.2.18) 

A;_oor+i; = -A;+oor_f ; (2.2.19) 
with equivalent relations for a and /3. Then 

= = = r_ (2.2.20) 

and barrier transmission is independent of direction. We also have 

phase {t+) = phase (i_) (2.2.21) 

and 

phase (r+/t+) = tt - phase (r_/t_) (2.2.22) 

with equivalent relations for a and 

If we now adopt the (to our minds) more useful convention, by absorbing 

factors of fc+oo and k^^o into the definitions of r and t then things simplify 
considerably: We restart the calculation by now defining 

^ ( ^ r^^\ + exp(+iA;_ooa;) , . 

J+{x^-oo)^t+ -== , (2.2.23) 

Mx ^ +OC) ^ 'M+^k^^^) + ^^eM-^k^^ ^^.2.24) 

Y k+oo Y ^"+00 

By comparing these two different forms for the asymptotic form of the Jost 
function we see that in this situation the ratios of the amplitudes are given 
by the much simpler formulae 

l:a+:p+^t+:l:r+. (2.2.25) 

We now have 

r+ = ^, (2.2.26) 
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and 

t+ = —. (2.2.27) 

We see that by putting the factors of ^yk±ao into the asymptotic form of the 
Jost functions, where they really belong, the formulae for r and t are suitably 
simplified. 

For completely arbitrary potentials, with no parity restriction (so the 
potential is neither even nor odd), a modified Wronskian analysis now yields 
(in analogy with that reported by Messiah 



= 1 - |r+|2; (2.2.28) 

|t_|2 = 1 _ |r^|2. (2.2.29) 

t_ = t+; (2.2.30) 

r+tl = -r^t*_; (2.2.31) 
with equivalent relations for a and (3. Then 

T+ = = = T_ (2.2.32) 

and barrier transmission is independent of direction. Because they are in- 
dependent of any overall scaling by a real number, also retain the previous 
results 

phase (t_|_) = phase (t_) (2.2.33) 

and 

phase (r+/t+) = n - phase (r_/t_) (2.2.34) 

with equivalent relations for a and (3. It is this modified set of conventions, 
because they have much nicer normalization properties, that we shall prefer 
for the bulk of the thesis. 

We shall now derive some very general bounds on |a| and which also 
lead to general bounds on the reflection and transmission probabilities 

i?=|r|2; r=|tp. (2.2.35) 
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2.3 Probability currents 

The expressions for reflection and transmission coefficients were based on the 
assumption that the intensity of a beam is the product of the speed of its 
particles and their linear number density. In classical physics, the assumption 
seems very natural, however, we should always be careful about carrying over 
classical concepts into quantum physics 



Definition (Unbound state): Provided V±oo = 0, the Schrodinger 



equation (1.3.1) can be solved for any positive value of energy, when 
E > 0. In addition, the positive energies can be shown to define a 
continuous spectrum. Nevertheless, the corresponding eigenfunctions 
do not vanish at infinity; their asymptotic behavior is analogous to 
that of the plane wave exp{ikx). More accurately, the absolute value of 
wave functions approaches a non-zero constant when x — > oo. 

Otherwise, the absolute value oscillates indefinitely between limits, 
one of which at least is not zero. It is clear that the particle does 
not remain localized in any finite region. This type of wave function is 
commonly applied to collision problems; the usual language is that one 
is dealing with an unbound state, or stationary state of collision 



2.4 Reflection and Transmission of Waves in 
unbound states 

The Schrodinger equation also can be analyzed in terms of the functions 
u and V, as defined by Messiah [51], and their complex conjugates u* and 
V*. Moreover, the Wronskian of any two such solutions is independent of x; 
especially, it takes on the same value in the two asymptotic regions. Actually 
our approach can be seen as equating these two values; we now derive a 
relation between the coefficients r+,t+,r_,t_, or their complex conjugates. 
Six such relations can be formed with the four functions u, v, u* and v* . From 
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what we have seen earher it is clear that they are very basic relations which 
must be maintained whatever the form of the potential function V{x) ^T\ . 
Specifically, we derive (in Messiah-like conventions) 

^) = A;_oo|t+r; (2.4.1) 

') = k+^\t.\^; (2.4.2) 

■U; (2.4.3) 

k-oorlU. (2.4.4) 



'-W{u,u*) 


— ^+oo(l ~ 


'-Wiv,v*) 


= fc_oo(l-r 




~ k-\-oot^ = k 







The equations ( 2.4.1[ ) and (2.4.2) are called the relations of conservation 



of flux. They should always be true, and this should be verified in spe- 
cial cases. This name comes from the following statements regarding the 
wave function ip of an unbound state in the asymptotic region. We let 
A exp{ikx) + B exp{—ikx) be the expression of the wave function ip in one 
of the asymptotic regions, for — oo case. 

The total flux of particles when passing a given point is the difference 
between the flux {hk/m)\A\'^ of particles traveling in the positive sense, and 
the flux {hk/m)\B\'^ of particles traveling in the negative sense. This flux is 
equal, to within a constant, to the Wronskian W{ilj,ip*) (5lj : 

-[\A\'-m = l-W{^,r) (2.4.5) 
m 2 m 

The equality of the Wronskian Wlipjtp*) at both ends of the interval 
(— oo,+oo), denotes that the number of particles entering the interaction 
region per unit time is equal to the number which leave it. In accordance 



with this interpretation, one or the other of equation (2.4.1) and (2.4.2) can 
be written as: 

incident flux — reflected flux = transmitted flux. (2.4.6) 

Considering the same interpretation, we now can define the transmission 
coefficient (transmission probability) T as follows: 

^ transmitted flux 

T = n ■ 2.4.7 

mcident nux 
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STATES 

We have in particular 

T+ = ^|t+ooP, T_ = ^|t_p. (2.4.8) 



This resuh shows that the absolute values of the two sides of equation (2.4.3) 
are equal, and one obtains the equality 

r_ = T+. (2.4.9) 



Thus the transmission coefficient of a wave at a given energy is independent 
of the direction of travel. This is the reciprocity property of the transmission 
coefficient. It is just as hard to traverse a potential barrier in one direction 
as in the other. 



The equality of the absolute values of the two sides of equation (2.4.4) 



coupled with the conservation relations (2.4.1) and (2.4.2), again yields the 



reciprocity relation (2.4.7) we also obtain relations between the phases of the 



reflection and transmission amplitudes: 



phase(t+) = phase(t_); 
phase = TT — phase ^—^ 



The most interesting point for these relations is the fact (not further inves- 
tigated in this thesis) that the phases are related to "retardation" effects in 
the propagation of the wave packets, with equivalent relations for a and f3. 

As previously, we can re-scale r and t by absorbing appropriate factors 
of \/k±oo, and so simplify the discussion as in the previous section. (We 
will not repeat the details of the analysis, as it is straightforward.) We shall 
now generalize some very general bounds on |a| and which also lead to 
general bounds on the reflection and transmission probabilities 

i?=|r|2; T=\t\^. (2.4.10) 
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Definition (Bound states): Provided V±oo = 0, when E < 0, the 



Schrodinger equation (1.3.1) has solutions only for certain particular 
values of energy forming a discrete spectrum. The eigenfunction iIj{x) 
corresponding to it — or each of the eigenfunctions when several ex- 
ist — vanishes at infinity. More accurately, the integral J \ tp{r)\'^ dr 
extended over the whole configuration space is convergent. There is a 
vanishing probability of finding the particle at infinity and the particle 
remains practically localized in a finite region. The particle can now 
be defined to be in a bound state fSTl. 



2.5 Bogoliubov transformation 

Definition (Bogoliubov transformation): This is a unitary trans- 
formation from a unitary representation of some canonical commuta- 
tion relation algebra or canonical anticommutation relation algebra 
into another unitary representation |T]. 



To see the import of this definition, let us consider the canonical commutation 
relation for bosonic creation and annihilation operators in the harmonic basis 



a, a' 



1. 



(2.5.1) 



Using this method we can derive a new pair of operators. 



b 



ua + va^; 
u a' + V t 



(2.5.2) 
(2.5.3) 



where the equation (2.5.2) is the hermitian conjugate of the equation (2.5.3). 



This transformation is a canonical transformation of these operators. It 
is easy to find the conditions on the constants u and v. For instance, the 
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transformation remains canonical by extending the commutator. 



[b, = [ua + va\ u*a^ + v*a] = 




(2.5.4) 



It can be seen that 

\u\^-\v\'^ = l (2.5.5) 

is the condition for which the transformation is canonical (This normahza- 
tion condition will occur and re-occur many times in the calculations which 
follow.) Finally we note that since the form of this condition is reminiscent 
of the hyperbolic identity 

cosh^ r — sinh^ r = 1 (2.5.6) 

between cosh and sinh, the constants u and v are usually parameterized as 

u = exp(i^) coshr; (2.5.7) 
V = exp(i^) sinhr. (2.5.8) 



2.6 Transfer matrix representation 

We can also investigate quantum mechanical tunneling by the so-called "trans- 
fer matrix method" or "transfer matrix representation" . Ultimately, of course, 
this is still equivalent to extracting the transmission coefficient from the so- 
lution to the one-dimensional, time-independent Schrodinger equation. As 
before, the transmission coefficient is the ratio of the flux of particles that 
penetrate a potential barrier to the flux of particles incident on the barrier. It 
is related to the probability that tunneling will occur [15j. We again consider 
a one-dimensional problem which is characterized by an incident beam of 
particles that is either transmitted or reflected as a result of scattering from 
an object |2T]. For current purposes it is easiest to work with potentials of 
compact support, where V{x) = except in some finite region [a,b]. 

As long as the potential V{x) is of compact support, it splits the space in 
three parts {x < a,x E [a,b],x > b). In both (— oo, a] and [b, oo) the potential 
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energy is zero. Moreover, in each of these two regions the solution of the 
Schrodinger equation can be presented as a superposition of exponentials by 



'ipiix) = Arexp{ikx) + Aiexp{—ikx) 
iPr{x) = Brexp{ikx) + Biexp{—ikx) 



X < a, and (2.6.1) 
x>b, (2.6.2) 



where Ai/r and Bi/r are at this stage unspecified, and k = \/2mE/h. But 
because ipL and ipR are solutions to the Schrodinger equation that can be 
extended to the entire real line, and because the Schrodinger equation is a 
second-order differential equation so that its solution space is two-dimensional, 
there must be some linear relation between the coefficients appearing in ipL 
and ipR — specifically, there must be a 2 x 2 matrix M such that 



(2.6.3) 



The 2x2 matrix M depends, in a complicated way, on the potential V{x) 
in the region [a, h]. In the transfer matrix approach we shall seek to extract 
as much information as possible without explicitly calculating M. 

To now derive amplitudes for reflection and transmission for incidence 
from the left, we put A^. = 1 (incoming particles), Ai = r (reflection), Bi = 
(no incoming particle from the right) and Br = t (transmission) in equations 









= M 




Br 







(2.6.1) and (2.6.2). We now derive 

iPl{,x) = exp{ikx) + rLexp{—ikx) 



(2.6.4) 



where is the left-moving reflection amplitude and on the right of the 
potential 

4'r{x) = t^exp^ikx). (2.6.5) 



where ti is the left-moving transmission amplitude. This tells us that 



tl 


= M 


1 












(2.6.6) 
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A 








r 


/ \ K 










1 = \tf + |7-|^ 



Figure 2.1: This shows an incoming flux of particles from the left, being 
partially transmitted to the right (with amplitude t), and partially reflected 
back to the left (with amplitude r) 



But since the Schrodinger equation (1.2.19) is real, the complex conjugate of 



any solution is also a solution. Therefore the solution which on the left has 
the form 

V'l = exp(— z/cx) + r\ exp(+zA;x) , (2.6.7) 
must on the right have the form 



^i?(a;) = tiexp(-ifca;) 



(2.6.8) 



and so we also have 



M 



(2.6.9) 
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These two matrix equations now imply 
M ' 



1 - rlvL 

But by conservation of flux we must have 

\tL\' + ItlI' = I. 



tl -tLTl 
tlTL tl 



(2.6.10) 



(2.6.11) 



We just have seen an important connection between reflection and trans- 
mission amplitudes. In addition, it is interesting to show how to derive the 
above equation by following. 

From the equation (2.6.4), we can see that this corresponds to a flux in 
the positive x direction. For x < a this is of magnitude 

h 



2mi 

h 

2mi 



d^/j dip* 

^ 7^ 

ox ox 

(exp(— z/cx) + r2 exp{+ikx)) (ik exp{ikx) 



-rM-^k.)) -complex conjugate 

h ( 

2ik - 2ik\rL\ 



2mi 

hk 

m 



1 - kLl 



(2.6.12) 



and for x > 6, we similarly derive from equation (2.6.5) the fact that we can 
write the flux corresponding to this equation is 

h ( 

(tlexp{—ikx) X ik{tLexp{ikx))) 



2mi 

— (tLexp{ikx) X —ik(tlexp{—ikx))) 

h 



2mi 

hk 

m 



ikltil + iklt^l 



\tL\'' 



(2.6.13) 



26 



2.6. TRANSFER MATRIX REPRESENTATION 



Definition: The probability current ^ of the wave function is 
defined as 

, d'\\) dip* 



h 



r 



(2.6.14) 



2mi \^ dx dx 
in the position basis and satisfies the quantum mechanical continuity 
equation 

^P(^,^) + ^^(^,i) = 0, (2.6.15) 
where p{x,t) is probabihty density ^n\. 



Since there is no time dependence in the problem, the conservation law 



in equation (2.6.14) implies that ^{x) is independent of x. Hence the flux 
on the left must be equal to the flux on the right, that is, we expect that 



hk 



m 



1 - {tlI 
1 - kil 



hk 



m 



\tL\'- 



\tL\'. 



therefore, 



so 



whence 



\tL\' + \rL\' 



1 — r\rL 1 — Ir^l 



M 



1 



tL -tirl 
-tlrL tl 



Utl -rl/tl 
-ri/tL l/tL 



(2.6.16) 
(2.6.17) 

(2.6.18) 



Similary, consider a wave moving in from the right 

exp(— z/cx) 



(2.6.19) 



which then hits the potential, is partially reflected and partially transmitted. 
In this case, on the right of the potential we have 



i^nix) = exp{—ikx) + rRexp{+ikx) , 



(2.6.20) 
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where is the right-moving reflection amphtude and on the left of the 
potential 

Mx) ^ tReM-^kx) , (2.6.21) 
where tR is the left-moving transmission amplitude. This tells us that 



(2.6.22) 



Again, since the Schrodinger equation is real, the complex conjugate of any 
solution is also a solution. Therefore a related interesting solution which on 
the left can be cast in the form 





= M 







1 







iPl{x) = t*Rexp{+ikx) , 
must on the right have the form 

i^nix) — ex.p{+ikx) + r*Rexp{—ikx) , 

whence 



1 


= M 












(2.6.23) 



(2.6.24) 



(2.6.25) 



But now these two matrix equations imply 



r-R/tR 
ryt*R 1/tR 



(2.6.26) 



Combining the information from left moving and right moving cases we have 
first that 

tL = tR. (2.6.27) 

So we again derive the equality of the transmission amplitudes. 
Similarly we see that 

(2.6.28) 



tR t*r 



implying 



_ _ * ^_L_ 



(2.6.29) 
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Note that we cannot in general deduce = tr. Indeed, in general this is 
false. 

So for any potential we have 

T ^ \tL\' ^ \tR\'; 



(2.6.30) 



implying (in the same manner as the previous argument) that the transmis- 
sion and reflection coefficients are independent on whether or not the particle 
is incident from the left or the right — and we have not made any assumption 
here about any symmetry for the potential V{x) itself. We conclude 



1/t* -rl/t* 
-TL/t l/t 



l/t* rn/t 
ryt* l/t 



(2.6.31) 



Note the key step in this general derivation: In any region where the 
potential is zero we simply need to solve 

d2 



il){x) — Eip{x), 



2m dx'^ 

for which the two independent solutions are 



exp{±ikx); 



V2mE 



h 



or more explicitly 



^/2mE 
exp ± I — - — X 
n 



To the left of the potential we have 

= a exp{ikx) + b exp{—ikx) , 
while to the right of the potential we have 

iPr{x) — c ex.p{ikx) + d exp{—ikx). 



(2.6.32) 

(2.6.33) 
(2.6.34) 

(2.6.35) 

(2.6.36) 



Even without knowing anything more about the potential V{x), the linearity 
of the Schrodinger ODE guarantees that there will be some 2x2 transfer 
matrix M such that 

(2.6.37) 



c 




a 




= M 




d 




b 
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This transfer matrix relates the situation to the left of the potential with the 
wave-function to the right of the potential. For this reason we shall now use 
this formalism, for instance, to think about the propagation of electrons down 
a wire (approximately one-dimensional) with V{x) used to describe various 
barriers placed in the path of the electron. Moreover, similar matrices also 
occur in optics, where they are referred to as "Jones matrices" . 

The components of the transfer matrix M will be some horrible nonlinear 
function of the potential V{x), but by linearity of the Schrodinger ODE these 
matrix components must be independent of the parameters a, b, c, and d. In 
some particularly simple situations we may be able to calculate the matrix 
M explicitly (see in particular the next chapter), but in general it will be a 
complicated mess. 

From the above discussion we now understand, from at least two differ- 
ent points of view, the basic concepts of transmission and reflection. The 
probability that a given incident particle is reflected is called the "reflection 
coefficient", R. While the probability that it is transmitted is called the 
"transmission coefficient", T [2T] . 

2.7 Discussion 

In this chapter, we have presented basic aspects of scattering theory in one 
dimension. For a one- dimensional model, only one of the three coordinates of 
3-dimensional physical space is explicitly involved. Specifically, we considered 
potentials of compact support, when the potential V{x) is mathematically 
zero outside of a finite interval. The situation where the potential is zero 
is referred to as "the free particle". These one-dimensional models provide 
solid examples exhibiting all the basic features and ideas needed to derive 
the properties of quantum states of definite energy E. 

A further step in our investigation is that we have just seen an important 
connection between refiection and transmission amplitudes. In particular, it 
is interesting to show how to derive them directly by using scattering theory. 
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Furthermore, we have now introduced the concept of transmission and re- 
flection. Wc called the probability that a given incident particle is reflected 
as the "reflection coefiicient" . While the probability that it is transmitted is 
called the "transmission coeflicient" . 

More importantly, we introduced the probability current to express the 
reflection and transmission coeflicients. The probability current is based on 
the assumption that the intensity of a beam is the product of the speed of its 
particles and their linear number density. It is then a mathematical theorem 
that this probability current is conserved. We then introduced important 
ideas of reflection and transmission of waves in both unbound and bound 
states. By considering reflection and transmission of waves in unbound states, 
we have seen that in principle they are completely specified by the potential 
function V{x). 

For instance, the linearity of the Schrodinger ODE guarantees that there 
will be some 2x2 transfer matrix. Moreover, this transfer matrix can be 
represented by investigating quantum mechanical tunneling by extracting 
the transmission coefficient from the solution to the one-dimensional, time- 
independent Schrodinger equation. 
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Chapter 3 

Known analytic results 



3.1 Introduction 

In this chapter we shall collect a number of known analytic results in a form 
amenable to comparison with the general results presented in subsequent 
chapters. We shall review and briefly describe the concept of quasinormal 
modes, and see how most of the concepts introduced here are important tools 
for comparing the bounds with known analytic results. Furthermore, we shall 
reproduce many analytically known results, such as the tunnelling probabil- 
ities and quasinormal modes [QNM] of the delta-function potential, double- 
delta-function potential, square potential barrier, tanh potential, sech^ po- 
tential, asymmetric square-well potential, the Poeschl-Teller potential and 
its variants, and finally the Eckart-Rosen-Morse-Poeschl-Teller potential. 

In the following, we shall first introduce the quasinormal modes, which 
are the modes of energy dissipation of a perturbed object or field. In partic- 
ular, the most outstanding and well-known example is the perturbation of a 
wine glass with a knife: the glass begins to ring, it rings with a set, or super- 
position, of its natural frequencies - its modes of sonic energy dissipation. In 
the absence of any damping, when the glass goes on ringing forever, we can 
call these modes normal. In the presence of damping, when the amplitude 
of oscillation decays in time, we call the modes quasi-normal [H]. 
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To a very high degree of accuracy, quasinormal ringing can be approxi- 
mated by 

^p{t) ^ exp{-iu"t) cos(cu't), (3.1.1) 

where ip{t) is the amphtude of oscillation, u' is the frequency and u" is the 
decay rate. We can express the quasinormal frequency in two numbers, 

uj = {iu',uj"), (3.1.2) 

or more compactly 

^(t) ^ exp(icut), (3.1.3) 

iu = uj' + ico", (3.1.4) 

where for ip{t) we are to understand that we are only interested in the real 
part. In our explanation here, u is generally referred to as the quasinormal 
mode frequency. The most interesting point is that it is a complex number 
with two pieces. One of them is a real part which describes the tempo- 
ral oscillation, and the other part is an imaginary part which describes the 
temporal exponential decay. Formally, quasinormal modes are most easily 
found by looking for complex frequencies where the transmission amplitude 
becomes infinite. 

In theoretical physics, a quasinormal mode is a formal solution of some 
linear differential equations with a complex eigenvalue. In black hole physics 
these linear differential equations typically come from linearizing the full Ein- 
stein equations. It is important to note that black holes have many quasinor- 
mal modes that express the exponential decrease of asymmetry of the black 
hole in time as it evolves towards the perfect spherical shape [14J . Experience 
obtained from black hole physics and related fields has shown that it is quite 
common for QNM to be approximately of the form 

LJn = a + inb + 0{l/n), (3.1.5) 

where a is a complex number called the "offset" and 6 is a real number known 
as the "gap" [98] • 
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3.2 Delta— function potential 

As a first approach to the comparison of the bounds with known analytic 
results we can start by studying in detail the concept of the delta-function 
potential. It is important to understand that the delta-function potential 
is one limiting case of a square well. It is a very narrow deep well, which 
can adequately be approximated by a mathematical delta function when the 
range of variation of the wave function is much greater than the range of the 
potential [2i] . 

The time-independent Schrodinger equation for the wave function iIj{x) 

is 



HiP{x) 



d2 



+ Vix) 



ip{x) = Eip{x), (3.2.1) 



2m dx^ 

where H is the Hamiltonian, h is the (reduced) Planck constant, m is the 
mass, E is the energy of the particle and the potential V{x) is the delta func- 
tion well with strength a < concentrated at the origin. Without changing 
the physical results, any other shifted position is also possible [25], i.e. 
For a delta function potential take 

V{x) = aS{x). (3.2.2) 

In this case the transmission coefficient is well known to be (see, for in- 
stance, [381 139] ) 

T= — ^. (3.2.3) 

ma 

1 + 



2Eh^ 



Quasinormal modes: T = oo when 

ma^ = -2Eh^, (3.2.4) 



that is 



ma^ , ima 



Note that there is only one pair of complex conjugate QNM. Because the 
width of the delta function is zero, the "gap" is infinite, and the other QNM 
are driven off to imaginary infinity. 
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Deriving the amplitudes: 

For completeness we will explicitly provide the calculations required to deal 
with the delta potential barrier — this is a textbook problem of quantum 
mechanics. Generally, the problem consists of solving the time-independent 
Schrodinger equation for a particle in a delta function potential in one di- 
mension |25j . 

We now consider particles entering from the left traveling to the right 
with E > encountering a potential of the form [18] 



V{x) = a5{x). (3.2.6) 
We look for solutions of the time-independent Schrodinger equation 

- ^^J^ + a6{xMx) = E^lj{x), (3.2.7) 
im dx^ 

with \l'(x,t) = il){x) exp{—iEt/h). For E > in the region x 7^ we have 
d^ipix) 2mE 



dx^ h? 



iIj{x) = -k^ij{x), (3.2.8) 



2mE „ „ 

with 1^ = \ ^ or E = (where n = —k ). 

The most general solution is 

i^Lix) = Aiexp{+ikx) + BLexp{—ikx), x<0, and (3.2.9) 
i^R{.x) = Ajiexp{+ikx) + Bjiexp{—ikx), x > 0. (3.2.10) 



Boundary Conditions: 

We must require ipi^) to be continuous everywhere, and d-?/'(x)/dx to be 
continuous (except possibly where V is infinite). Thus, at x = 0, we have 
'0-l(O) = '?/'_r(0) which implies that Al + Bl = Ar + Br. 
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Figure 3.1: This diagram describes the scattering process at a delta-function 
potential of strength a. The amplitudes and direction of left and right moving 
waves are indicated. We indicate, in red, the specific waves used for the 
derivation of the reflection and transmission amplitudes [19] . 



37 



CHAPTER 3. KNOWN ANALYTIC RESULTS 



Also notice that for e 



we have 



da; 



2m 
2m 



+e 



{V{x) -E)ij{x)dx, 
{a6{x) — E)iIj{x) dx 



2ma 



(3.2.11) 



Thus, 



dx 



dx 



2ma 



2ma 



V^l(O), 



(3.2.12) 



which imphes that [ikAn — ikBn] — [ikAi — ikBi] 



2ma 



(Al + Bl). 



Let us consider the second of these equations, which follows from inte- 
grating the Schrodinger equation with respect to x. The boundary conditions 
thus give the following restrictions on the coefficients |T9j 



AL + BL = An + Br- 
2ma 



ikAji — ikBpi — ikAi + ikBi 



{Al + Bl). 



(3.2.13) 
(3.2.14) 



Reflection and Transmission: 

We shall see how to find the amplitudes for reflection and transmission for 



incidence from the left, by putting in the equations (3.2.13) and (3.2.14); 
Al = 1 (incoming particle), Bl = r (reflection), Br = (no incoming 
particle from the right) and Ar = t (transmission). We obtain 

l + r = t, (3.2.15) 

and 



ik{l - t 



2ma 
1^ 



t. 



(3.2.16) 



To find the amplitudes for reflection and transmission for incidence from the 
left, we now solve for r and t: 

1 



t 



ima 



(3.2.17) 
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and 



-1 



1 + 



ikfi' 



2 • 



ma 



Now we can derive the probability for transmission and reflection 
the transmission coefficient and the reflection coefficient) 



T = \t\ 



1 



and 



R 



1 1 1 
ma 

1 



(3.2.18) 
(given by 

(3.2.19) 



1 + 



m^a^ 



This confirms the result we previously quoted, and by looking for 
the transmission amplitudes, confirms the locations of the QNM. 



(3.2.20) 

poles of 



3.3 Double-delta-function potential 

For the double delta function 

V{x) = a{6{x - L/2) + 6{x + L/2)}, 
the transmission coefficient is known to be |17] 



1 + 



2ma 



cos (kL) 



1 / 2ma 
2\¥k 



sin {kL) 



2 • 



(3.3.1) 



(3.3.2) 



It is an easy exercise to check that this satisfies the bounds (6.1.6) and (6.1.8) 
that we shall derive later on in this thesis. 



Quasinormal modes: T = oo when 



2ma If 2ma^ ^ 



n?k ' ' 2 \ t9k 



cos {kL) + - -— sin {kL) = ±i, (3.3.3) 



39 



CHAPTER 3. KNOWN ANALYTIC RESULTS 



which leads to quite horrible algebra, so that there is no explicit formula for 
the QNM. Implicitly, working with the transmission amplitude: 

t = oo {k - ikof + kl exp{2ikL) = 0, (3.3.4) 

so that 

ik^ 



ex.p{2ikL) = + — j . (3.3.5) 
Deriving the amplitudes 

To derive the transmission and reflection amplitudes for the double-delta- 
function potential, we now start by considering the potential [E] 

V{x) = a{6{x - L/2) + 6{x + L/2)}. (3.3.6) 

For a particle incident from the left we now have 

{exp{+ikx) + r exp{—ikx) (at x < —L/2); 
Aexp{+ikx) + B exp{-ikx) (at - L/2 < x < L/2); 
t exp{+ikx) (at x > L/2). 

(3.3.7) 

Note that now we also have to explicitly consider the region between the two 
delta-functions. Applying the same sort of boundary conditions we now have 
four equations. From continuity at x = —L/2 we have 

exp{-ikL/2)+rexp{+ikL/2) = Aexp{-ikL/2) + B exp{+ikL/2), (3.3.8) 

while continuity at x = +L/2 implies 

Aexp{+ikL/2) + B exp{-ikL/2) = texp{ikL/2). (3.3.9) 

Integrating across the delta functions leads to 

— {ik[exp{-ikL/2) - r exp{+ikL/2)] - ik[Aexp{-ikL/2) - B exp{+ikL/2)]} 
+a [A exp{-ikL/2) + B exp{+ikL/2)] = 0, (3.3. 10) 
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and 

—{ik[exp{-ikL/2) - r exp{+ikL/2)] - ik[Aexp{-ikL/2) - B exp{+ikL/2)]} 
= +a[Aexp{-ikL/2) + B ex.-p{+ikL/2)]. (3.3.11) 

We rearrange this to obtain 

{ik[exp{-ikL/2) - r exp{+ikL/2)] - ik[Aexp{-ikL/2) - B exp{+ikL/2)]} 
= -— [Aex^i-ikL/2) + 5exp(+iA;L/2)], (3.3.12) 

and 

{ik[ex^{-ikL/2) - r exp(+iA;L/2)] - ik[Aex^{-ikL/2) - B ex^{+ikL / 2)]} 
= 2kQ [Aex^{-ikL/2) + 5 exp(+iA;L/2)]. (3.3.13) 

Some further rearrangements lead to 

-—{ik[AexT^i+ikL/2) - B exp(-ikL/2)] - ik[texp(+ikL/2)]} 
2m 

+« [t exp(+iA;L/2)] = 0, (3.3.14) 

and 

— {ik[AexY>{+ikL/2) - B exp{-ikL/2)] - ik[texp{+ikL/2)]} 
2m 

^a[t exp(+iA;L/2)] . (3.3.15) 

Finally we have 

{iA;[Aexp(+iA;L/2) - Sexp(-iA;L/2)] - ik[tex^{+ikL/2)\} 

= 4^[tcxp(+^fcL/2)], (3.3.16) 

It is very useful to reduce clutter by defining 

A:o = (3.3.17) 



then 



{ik[Aex^{+ikL/2) - B exp{-ikL/2)] - ik[t exp{+ikL /2)]} 

^2ko[texp{+ikL/2)]. (3.3.18) 

41 



CHAPTER 3. KNOWN ANALYTIC RESULTS 



in which case our four boundary conditions become: 
ex-p{-ikL/2) + r exp{+ikL/2) = Aexp{-ikL/2) + B exp{ikL/2), (3.3.19) 
Aex.p{+ikL/2) + B ex.p{-ikL/2) = tex.p{+ikL/2), (3.3.20) 

{ik[ex.p{-ikL/2) - rex.p{+ikL/2)] 

-ik[Aex.p{-ikL/2) - B ex.p{+ikL/2)]} 
= 2ko[Aex-p{-ikL/2) + B exp{+ikL/2)], (3.3.21) 

{ik[Aexp{+ikL/2) - B exp{-ikL/2) - ik[texp{+ikL/2)]]} 

^2ko[texp{+ikL/2)]. (3.3.22) 

These are four simultaneous linear equations for four unknowns: r, A, B, 
and t (in terms of the known quantities k, ko, and a). These can be solved, 
either by direct calculation or by Maple or something similar. A httle work 
then leads to 

k^ 

^ " {k - ikof + kl exp{2tkL) ' (3.3.23) 

and 

^ 2iko[kcos{kL) - kosm{kL)] 

(k-iko^ + kieM'^ikL) ■ ^ ^ 

Note that r/t is pure imaginary, in agreement with our general argument 
regarding definite parity potentials. Furthermore note that 

R=\r\^ oi[k cos{2ka) - ko sin(2A;a)]^ (3.3.25) 

and so R — whenever 

k 

tan(2A;a) = — . (3.3.26) 

That is, the system exhibits "transmission resonances" where T ^ 1 and 
R 0. If we work at fixed energy then these resonances occur at equally 
spaced spatial separation for the two delta functions, namely: 

Oresonance = ^{^^^"^ / kg) + mr}; U & Z. (3.3.27) 
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If we hold a fixed and vary k then the location of the resonances is determined 
by the transcendental equation 



tan(2/ca) 



k 

ko' 



(3.3.28) 



The existence of "transmission resonances" in one-dimensional scattering is 
in fact widespread, it is not specific to this particular example. A brief 
computation leads to the explicit transmission coefficient 

k^ 



T 



k^ + 4k'^[kcos{2ka) - kosm{2ka)Y 



or 



T 



1 + 



2ma 
~Wk 



cos (kL) 



2ma 
~h?k 



sin (kL) 



2 ' 



(3.3.29) 



(3.3.30) 



which agrees with the preceding argument on the location of the transmission 
resonances. Finally note T{k — > 0) — > and T{k oo) 1. After we look 
at one more illustrative example, we will turn to the issue of obtaining some 
general theorems governing one-dimensional scattering. 



3.4 Square barrier 

Let us now start by introducing another useful potential barrier — the square 
barrier |llj : 

VW=(^° 'f" "^^^^'^ (3.4.1) 
I (otherwise). 

In our particular case, we set Vq > 0. Tunneling over a square barrier is an 
elementary problem which however is not always discussed in the textbooks. 
In contrast, tunneling under a square barrier is much more popular. The 
exact transmission coefficient is known to be 

T ^ (3.4.2) 

E{E - Ve) + sin2(^2m(E - Ve)L/h) 
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For more details see (for example) Landau and Lifshitz [37], or Schiff [15] . 
We can re-write this as 

^ (3.4.3) 



mV^L^ sm\^2m{E -Ve)L/h) ' 
^ ^ 2En? 2m{E -Ve)L'^/n? 

Quasinormal modes: T = oo when the numerator is nonzero and 

sin(V2m(E-K)V^) = ±2^E{E - K)/K , (3.4.4) 

which leads to hopeless algebra. Although E = Ve is one solution, it corre- 
sponds to the numerator vanishing and is not a transmission pole. There is 
no simple explicit formula for the QNM. 

Deriving the amplitudes: 

We shall now consider a particle of mass m and energy E > Q interacting 
with the simple square potential barrier. Let us consider ip{x) in the regions 
to the left and to the right of the barrier, we have 

- (3.4.5) 

where k"^ = 2mE/h^ . We choose the following solution of the above equation 
to the left of the barrier (i.e., a; < 0) 

iIj{x) = exp{ikx) +rexp{—ikx). (3.4.6) 

The composition of this solution is a plane-wave of unit amplitude traveling 
to the right [since the time dependent wave function is multiplied by a factor 
exp{—iEt/h)], and a plane wave of complex amplitude r traveling to the 
left. Moreover it should be stressed that the first plane wave is incoming 
particle. Indeed, the second plane wave is a particle reflected by the potential 
barrier. Therefore |rp is the probability of reflection. This can also be seen 
by calculating the probability current in the region (x < 0), which takes the 
form 

^i = ^i^-\rn (3.4.7) 

fit 
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We choose the following solution to equation (3.4.5 ) to the right of the barrier, 
that is, for x > L: 

ip^x) = texp{ikx). (3.4.8) 

We have seen that this solution consists of a plane wave of complex amplitude 
t traveling to the right. This implies that this solution can be interpreted as a 
particle transmitted through the barrier. Consequently, |tp is the probability 
of transmission. In fact we can write the probability current in the region 

X > L as 

A = -\t\"- (3.4.9) 
m 

If we set = then we derive 

|r|2 + |t|2 = i. (3.4.10) 

At this point is easy to see that inside the barrier (i.e.,0 < a; < L), the 
wavefunction '^{x) satisfies 



d^^ 
dx^ 

where 



(3.4.11) 



, 2m(E-1/o) 
= \, (3.4.12) 

We consider the case where E > Vq. In addition, the general solution to 



equation (3.4.11) inside the barrier takes the form 

ip^x) = Aexp{+iqx) + B exp{—iqx), (3.4.13) 

where q = ^j2m{E — Vq)/^?. From the continuity of ip and dip/dx at the 
left edge of the barrier {i.e.,x = 0) we derive 

1 + r = A + B, (3.4.14) 
fc(l-r) = q{A~B). (3.4.15) 

Moreover, continuity of ip and dip/dx at the right edge of the barrier, for 
{x = L) gives 

Aexp(+igL) + Sexp(-igL) = texp(+iA;L), (3.4.16) 
g(v4exp(+igL) - 5exp(-zgL)) = ktexp{+ikL). (3.4.17) 
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It is now relatively easy to see that, (after considerable algebra), the above 
four equations yield 

, ,2 {e -qysm''{qL) 



R 
T 



4k^q^ + {k^ - q^ysin\qL)' 



(3.4.18) 
(3.4.19) 



3.5 Tanh potential 

For a smoothed step function of the form 



V{x) 



+ 



X 



tanhf ^ 



(3.5.1) 



2 2 

the reflection coefficient is known analytically to be (see, for instance, |37j): 



( sinh[7r(A;_oo - k+^)L/2] \ ^ 
Vsinh[7r(A;_oo + Koo)L/2] J 



(3.5.2) 



This certainly satisfies the general bounds (6.5.2)-(6.5.3) that we shall derive 



later on in this thesis, see chapter |6| and as L — approaches and saturates 
the bound. 



Quasinormal modes: T = oo when 

sinh[7r(A;_oo + k+^)L/2] = . 



that is 



This leads to 



that is 



so that 



sin[i7r(A;_oo + k+oo)L/2] = . 
i7r(A;„oo + k+oo)L/2 = mr , 



k+oo + \l kl^ + AmV+oo/h'^ = -i2n/L , 
+mV+ooL , n 



h? n ~^ L 



(3.5.3) 

(3.5.4) 
(3.5.5) 
(3.5.6) 
(3.5.7) 
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Equivalent ly 



h?n L 



Note the asympototic spacing as n — > oo: 



n 



(3.5.8) 



(3.5.9) 



Note that as L — > all the QNM are driven to imaginary infinity — this is 
compatible with the behaviour of the step potential for which there are no 
QNM. 



3.6 Sech^ potential 

For a sech^ potential of the form 

V{x) = VeSech^{x/L) , (3.6.1) 
the transmission coefficient is known analytically to be (see for example j37j): 

siYih^[K^/2mEL/h] + cos2[|7r^l - ^mVeL'^/h] ' 

provided SmVeL"^ < h^. This satisfies the general bounds derived later in 

this thesis, both the bound T > sech^ ^ dx^ , and the separate bound 

T > (4/i;+oo^-oo)/(^+oo + k^ooY- (Though proving this is somewhat tedious.) 
Start by noting that for this sech potential 



T > ianh^ [K\/2mEL/h] 
and use the inequality (x > 0) 



tanh X > 



X 



1 + x' 



> sech^ (1/^) • 



Then 



T > sech^[h/{7iV2mEL)] , 
sech 




vr V h Sm\Ve\L'^ 
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Provided that the extremum is a peak, l^eak > we can use the bound 
SmVpeakL'^ < to deduce 



T > sech 



m 2L|Vpeak| 

2E h 



(3.6.7) 



This is the particularization of the bound T > tanh^^ {} dx^ to the 
present case. If T4 < we need a different analysis. 



Quasinormal modes: T — oo when 



sinh^lW 2mEL/h] + cos 

which leads to 

— sin^ [inV 2mEL / }t\ + cos^ 



\^:^Jl-%mVeL'^/h^ 



= 0, (3.6.8) 



But then 



sm[i^:^/2mEL/h] = ±cos[^7rVl - %mVeL'^/h\ 



0. (3.6.9) 



(3.6.10) 



and so 



cos[m\/2mEL/h — 7r/2] = ±cos 



Therefore 



i7rV2mEL/h - 7r/2 = ±^7r y/ 1 - SmVeL'^ /h? + ttu, 



leading to 



Finally 



(3.6.11) 



(3.6.12) 



(3.6.13) 



(3.6.14) 
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Again note the asymptotic spacing as n — >■ cx) 



K^i-. (3.6.15) 



Note that if Vg is big, the offset term becomes real 

k = ±^^%mV,Lyh? - 1 + . (3.6.16) 

Finally note what happens as L becomes small, 

kn = T 2imV,L/h^ + + 0{L^) . (3.6.17) 

In the limit only one pair of QNM survive 

k = 2im\im[VeL]/h^ , (3.6.18) 

the others being driven off to infinity. This agrees with the result for the 
single-delta-function potential. 

Derivation of the amplitudes (sketch): 

Let us sketch how to determine the transmission coefficient for a potential 
barrier defined by the formula 



V{x) = Ve sech2(x/L). (3.6.19) 

Comparing this with the analogous bound state computation, it is necessary 
merely to alter the sign of Ve and to regard the energy E now as positive. A 
calculation similar to that used for deriving the bound states when E < 
now gives the solution 



-ik/L - s, -ik/L + s + 1, -ik/L + 1, ^(1 - 0] 



(3.6.20) 
where 

i = tanh(a;/L), (3.6.21) 
49 



CHAPTER 3. KNOWN ANALYTIC RESULTS 



k = ^/{2mE)lh , (3.6.22) 

(3.6.23) 



- 1 + 



8mV, 

1 - 



and F[__, __, __, __] is a hypergeometric function. This solution satisfies the 
condition that, as a; ^ oo (i.e. as ^ ^ 1, (1 — .^) ~ 2exp(— x/L)), the wave 
function should include only the transmitted wave (~ exp{ikx/L)). The 
asymptotic form of the wave function as x ^ — oo, (^ —1) is found by 
transforming the hypergeometric function with the aid of formula 

, ^T(ikL)T(l-ikL) , T(-ikL)T(l-ikL) 

ip ~ exp(— iKXj— -— — — r h exp(«A;x): 



r(-s)r(l + s) 'T{-ikL- s)T{-ikL + s + l) 

(3.6.24) 

Taking the squared modulus of the ratio of coefficients in this function, we 
obtain the following expression for the transmission coefficient T = 1 — R: 

y ^ smh\7r.yj2^L/h) ^ ^5) 

smh\7i^{2mE)L/h) + cos2[i7rv/(l - SmVeLyh^)] ' 

when SmVeL'^ / < 1. 



3.7 Asymmetric Square- well potential 

For the asymmetric square well 

{Vi, X < a; 
^2, a<x<h] (3.7.1) 
V"3, h<x. 

We now define ki = ^J^miE — V.;)/h. The transmission coefficient is (see for 
example [51]): 

Akik^k^ 

" {ki + A;3)2A;i + [/tf/tf + kl{kl - kf - k^)] sin\k2L) ' 



Then 



T > ,,2 ' \ ,2 • 3.7.3 

[kl + kik^Y 
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Similarly to the case for the symmetric square well, the transmission prob- 



ability for the asymmetric square well oscillates between the bound (6.6.10) 
that we shall subsequently derive, and the unitarity limit T = 1. For certain 
values of the width of the well [k2L = (2n + l)vr/2] the transmission coeffi- 
cient saturates the bound thus showing that this bound cannot be improved 
unless additional hypotheses are made. Because V^^c 7^ V+oo the bound 



T > sech^ ^ ?9 dx ) , is not applicable, at least not without modification 



oo 



from its original form. 

Derivation of the amplitudes (sketch): 

As in the problem of the potential step, we build the eigenfunction of the 
form |5T] 

{S exp{—ikix) X < a; 
P exp{—ik2x) + Qexp{ik2x) a < x < b; (3.7.4) 
exp(— ifcax) + Rexp{ik3x) b < x. 

The continuity conditions at points a and b give the values of R, Q, P, and 
S. Without entering into the specific details of the calculation, we simply list 
the results concerning the quantities R and S. We use the following notation 
and conventions: 

a = -L, 6 = 0, K= 

h h ^ ki 



We now derive 



^ _ k2{k'i - ki)cos{k2L) + i{kl - k^ki)sm.{k2L) ^ 
^2(^3 + ki)cos{k2L) — i(A;| + A;3A;i)sin(A;2i^) ' 

S = exp{-ik,L) ^^^^^ ^ fc^)cos(A;2L) - t{kl + k2k,)sm{k2L) ' ^^'^-^^ 
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The wave is in general only partially transmitted, and we can define a trans- 
mission coefficient (note that we are now using Messiah-like conventions) 

^=7:^1^1' = .... , r.M. . ..^,.■.2.. (3-7.7) 



A;3 ' ' kiiks + fci)2cos2(fc2L) + {ki + kskiysm\k2L) 



Akiklk-i 



{ki + k^fkl + [klkl + kl{kl -kl- kl)] sm\k2L) ■ 

3.8 Poeschl— Teller potential 

The Poeschl-Teller potential is most commonly written (see, e.g., [59] ) 

V{x) = Vq cosh^/i{tanh[(x — fiL)/L] + tanh/i}^ , (3.8.1) 

we have 

V.oo = Voe-^>'; Kxtremum = 0; F+oo = K)e+'^ (3.8.2) 

The transmission coefficient is [59] 

2sinh[(7rA;_oo-^) (^rk+ooL)] 



T 



cosh[7r(A;_oo + k+oo)] + cos 



7rJl + ^^^^cosh> 



(3.8.3) 



It is now a straightforward if tedious exercise to check this analytic result 
against all the bounds derived in this thesis. 

This Morse-Feshbach presentation of this potential [59], as given above, is 
rather difficult to interpret — it is much easier to ffist translate the potential 
(x — i> X + L) to obtain 

V{x) = Vo cosh^ /i{tanh [x/L] + tanh fi}^ , (3.8.4) 

expand 

V{x) = Vo cosh^ /i {tanh^[x/L] + 2 tanh tanh [x/L] +tanh^^} , (3.8.5) 

and then re-group terms as 

V{x) = Vocosh^/i |— sech^[x/L] + 2 tanh tanh [x/L] + 2 — sech^/i} , 

(3.8.6) 



52 



3.8. POESCHL-TELLER POTENTIAL 



to see that this is simply a Hnear combination of the "sech^", "tanh", and 
"constant" potentials. So without loss of generality we can re-write the 
Poeschl-Teller potential (with new definition for Vq = — cosh^ fj, [Vojoid) as- 

V{x) = Vo sech^ [x/L] + tanh [x/L] , (3.8.7) 

with 

V{-oo) = -Ko; ^(0) = Vo; V{+oo) = +V^;, (3.8.8) 

in terms of which the analytically known transmission probability is 

2 sinh (irk-ooL) sinh (tt/c+oo-^) 



cosh [7r(/c_oo + k^oo)L] + cosh 



TT./l-^ 



(3.8.9) 



Of course we have already seen that this has at least two much simpler limits: 
the sech^ and tanh potentials. In particular if we let Voo — ^ and play with a 
few hyperbolic identities we recover the results for the sech^ potential, while 
if we let Vo — and play with a few hyperbolic identities we recover the 
results of the tanh potential. 



Quasinormal modes: T = oo when 



cos[i7r(/c_oo + k+oo)L] + cos 



SmVoL^ 



0. 



That is 



whence 



i7r(A;_oo + k+oo)L = ±7rW 1 - 



+ {2n + l)7r, 



, , , .1 A SmVoL^ .2n + l 
(A;_oo + /c+oo) = ^^Xy^ p — X — ■ 



We now rearrange this as 



(3.8.10) 



(3.8.11) 



(3.8.12) 



SmVoL'^ .2n + l 
2n + 1 ± y/l-8mVoLyh^^ 



(3.8.13) 
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and note that it has appropriate hmits for the tanh and sech^ potentials. 
Finally, because this is a simple quadratic equation, we solve for k{n) to 
obtain 



k{n) 



2mV^L/h^ 



(2n + 1) ± v^l - SmVoLyn^ 



+ 



2n + l±Jl- 8mVoL'^/h^ 



2L 



By shifting n by ±1 as appropriate we can re-write this as 



k{n) 



2mVooL/h^ 



2n± 



2n± 



^ + 



l-8mVoL^/h'^ -1 



(3.8.14) 
\ 



2L 



J 



^l-8mVoL^/h^ - 1 

(3.8.15) 

This now has the appropriate limits to reproduce both tanh and sech^ quasi- 
normal modes. Note that asymptotically 



k{n) ^ -^-^ + 0(l/n) 



(3.8.16) 



in accordance with the general suspicions based on black hole QNMs [98]. 



3 . 9 Eckart— Rosen— Morse— Poeschl— Teller 
potential 

Many of the potentials commonly encountered in the literature are actually 
the same quantity in disguise. To start with, consider the following three 
potentials: 
Eckart (1930): 

V(x) = -^-^^^; e = -exp(2x/a). (3.9.1) 

Rosen-Morse (1932): 

V{x) = Btanh{x/d) - Csech^{x/d). (3.9.2) 
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Poeschl-Teller (1933): 

V{x) = VocoshV{tanh([x -/xL]/L) +tanh/i}l (3.9.3) 

To see that all three of these potentials are actually the same, note that: 

4^ _ 4 _ 4 

(1-0^ ~ (^-1/2 + e+i/2)2 - [exp(-x/a) + exp(a;/a)]2 ' 

= -—^— = sech\x/a). (3.9.4) 
cosh [x/a) 

Furthermore 

2^ 1 + ^ l-exp(2x/a) 



^^1-^ 1-^ 1 + exp(2x/a) 

exp(— x/a) — exp(a;/a) , , , . 

' ' i^v / ; -tanh(x/a). (3.9.5) 



exp(— x/a) + exp(a;/a) 
This is enough to show 

(Eckart) <^==^ (Rosen-Morse) 

In fact in the article of Rosen and Morse ^7], they cite Eckart |69], and de- 
scribe Eckart 's potential as begining "somewhat like" their own, but without 
noticing that the two potentials are in fact identical up to trivial redefinitions 
of the parameters. 

Now, for the Poeschl-Teller potential, note that by a trivial shift x — > 
X + /iL we have 

V{x) Vq coshV{tanh(a;/L) + tanh /i}^ (3.9.6) 

which we can without loss of generality relabel as 

V{x) {tanh(x/L) + 

= Vi{taiih^{x/L) + 2Dtanh(a;/L) + D^}, 

= 1/i{-sech2(x/L) + 2L)tanh(a;/L) + + 1}, 

= V2 sech2(x/L) + 1/3 tanh(x/L) + V^. (3.9.7) 
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This is enough to show 

(Poeschl-Teller) <^==^ (Rosen-Morse) 

and so all three potentials are completely identical up to trivial relabeling of 
the parameters and a shift in the zero of energy. In fact, including the offset, 
all three of these can be written in any one of the four general forms below: 



V{x) = A + Btajih{x/a + 9) + Ctanh.^{x/a + 

= Ao + [Bo + Co tanh(x/a + 6*)]^ 
Bi + Ci tanh(x/a) 



A 



Ao + 



i?2 + C2 tanh(x / a 
El + Fi exp(-2x/a) 
E2 + F2 exp(-2x/a) 



(3.9.8) 



Note that there is some redundancy here, but it is a useful redundancy. 
It makes it clear that the Eckart-Rosen-Morse-Poeschl-Teller potential is 
generally a Mobius function of the variable exp(— 2a;/a). Thus implies that 
without loss of generality we can set E1F2 — E2F1 as some convenient con- 
stant. 

The "best" of these equivalent forms is arguably the Mobius form: 



V{x) = Vo + Vi 



A + B exp(-2a;/a) 
C + D exp(-2x/a) 



(3.9.9) 



Comment: Many authors seem to use the phrase "Poeschl-Teller potential" 
only to refer to the special case 

V{x) = \/o/cosh2(a;/a) = Vo sech^ {x / a) . (3.9.10) 

This is historically inaccurate [;65j, and we will have more to say on this later. 
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3.10 Mobius potential 



Overall, the "best" general version of the potentials considered above is prob- 
ably the Mobius form 



_C + D exp(-2x/a) 
If you want to solve the Schrodinger equation 



(3.10.1) 



2m 



+ Vip = Eip, 



or 



(3.10.2) 



(3.10.3) 



then, absorbing the fi^ /2m into a redefinition of E, Vq and Vi, what you need 
to do is to solve 

2- 



i^" + {E-Vq-Vi 



A + B ex-p{-2x/a) 



V' = 0. 



(3.10.4) 



_C + D exp(-2x/a) 

This can either be solved "by hand" , or with the aid of symbolic manipu- 
lation packages. For instance, Maple still needs a little help. Without loss 
of generality, rescale x — > xa/2, and rescale Vq — > 0,Vi — > Vq, and demand 
AD — BC — 1. This makes the new E and Vq dimensionless and we have: 



'iI^"+{E-Vq 



A^B exp(-a;) 



= 0. 



(3.10.5) 



C + D exp(-x) 

This has an explicit solution in terms of hypergeonietric functions. Maple 
(after a little bit of convincing) gives 

^l){x) = (Cexp(a;) + D) 



X < Ao exp ( \/Vo^ - Ex] 2^1(^1,^2; As; exp(a;) 



'^2 



D 



+ S.exp|-,/Vog-^.),F,(i3„i3.;i33;-^^ 



(3.10.6) 
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where the suppressed arguments on the hypergeometric functions are: 



A, 



and 



(3.10.7) 



2 + V^°^-^-2V'+C^ + V^°S^-^' 



, 52 



(3.10.8) 



(3.10.9) 



1/52 1 4Vo A^ 



(3.10.10) 



^ 2 V 2V C2D2 V C2 
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(3.10.11) 



(3.10.12) 



By looking up tabulated properties of the hypergeometric function one can 
now determine the bound state eigenvalues, reflection and transmission co- 
efficients, etc. 



In general, the Mobius potential exhibits both bound and free states. (As it 
must, since after all the way we have derived it is by showing that it is equiv- 
alent to any of the Manning-Rosen, Poeschl-Teller, or Eckart potentials.) 
When the potential energy has a minimum but goes asymptotically to some 
higher finite value at x — +00 and —00, then some of the allowed energies 
will be discrete values, corresponding to states for which the particle is bound 
in the potential valley. For other ranges of energy, higher than the minimum 
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of the two asymptotic values, all energies will be allowed, the particle being 
free to travel to infinity. 

For instance, consider the wave functions and allowed energies for a parti- 
cle of mass m in a Mobius potential that is written in Poeschl-Teller form [59J 

V{x) = VoCOshV{tanh[(x-/iL)/L] +tanh(/i)}l (3.10.13) 

For Vo positive, this potential field has its minimum value = 0) at x = 0. 
As X is increased positive, the potential increases to an asymptotic value 
Voexp(2/i) for x +oo; as x is made negative, V also rises to an asymptotic 
value Voexp(— 2/i), for x — oo. Classically, since in this form the potential 
is constrained to be positive semidefinite, the particle could not have a nega- 
tive energy; for energies between zero and Voexp(— 2/i), (yU > 0), the particle 
would oscillate back and forth in the potential valley; for energies between 
Vbexp(— 2/i) and Voexp(2/t), the particle could come from — oo, be reflected 
by the potential rise to the right of the minimum, and go back to — oo and 
for energies greater than Voexp(2/i), the particle could move from — oo to 
-|-oo or from +oo to — oo [59] . 

3.11 Other potentials: 

Furthermore, we would like to at least mention some other potentials: 
Morse (1929) 

V{x) = Vo{l- exp(-[x - xo]/a)f. (3.11.1) 

The Morse potential is actually a somewhat odd limit of the Mobius potential 
as various parameters go to unity or zero. In terms of the Mobius potential 
above we need Vq ^ 0,Vi ^ Vq, A = B 1,C 1, D 0,x x — Xq. 

Manning Rosen (1933) 

V{x) = Bcoth{x/d) -Ccosech^ix/d). (3.11.2) 
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Warning: The relevant citation [26] is only an abstract in a report of a 
conference. To find it with online tools such as PRO LA look up Phys. Rev. 44 
(1933) 951, and then manually scan for abstract ^ 10. The form actually 
given in the abstract is 

Y(^) _ A exp{-2x/b) exp{-x/b) 

exp(-x/6)]2 + ^ 1 - exp(-x/6) ' ^^-^^'^^ 

which you can manipulate into the form above by noting 

exp(— a;/6) l+exp(— 



1 + 2 



1 — exp(— 1 — exp(— 

exp(a;/26) + exp(— x/26) 



coth(a;/26), 



exp(a;/26) — exp(— a;/26) 

(3.11.4) 



and 



coth^z = 1 + cosech^^. (3.11.5) 

Note that the Manning-Rosen potential can be obtained from the Eckart 
potential by the substitution 

^ = — exp(2x/a) +exp(— 2x/a). (3.11.6) 

In particular, Manning-Rosen can be written in the form 

V{x) = A + Bcoth{x/a) + Ccoth^{x/a) , 

= Ao + [Bo + Co coth{x/a)]^. (3.11.7) 

We can get this from the general Mobius form of the Eckart potential by 
appropriately choosing the parameters. 

Hulthen (1942) 

N exp(— x/a) 

V{x) = Vo- , y 3.11.8 

1 — exp(— x/a) 

The Hulthen potential [TO] is actually a special case of the Manning-Rosen 
potential {A = 0). We can also get this from the general Mobius form of the 
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Eckart potential by appropriately choosing the parameters. 
Tietz (1963) One version of the Tietz potential [71] is: 

^ ^ ° (^{sinh, cosh, exp}(x/a) ) ^ ^ 

We can get this from the general Mobius form of the Eckart potential by 
appropriately choosing the parameters. 

Hua (1990) 

V{x) = vJ ^--M-2x/a) y^ (3_,,_,o^ 
^ ^ ^\l-q exp{-2x/a) J ^ ^ 

We can get this |117j from the general Mobius form of the Eckart potential 

by appropriately choosing the parameters. We note 

exp(x/a) — exp(— x/a) 



Vix) = V, 



exp(x/a) — q exp(— x/a) 



yl smh{x/a) (3.11.11) 

(1 + g)sinh(x/a) + (1 — g)cosh(x/a) / 



If g > define 1 + q = cosh^^ and 1 — q = sinh6'. 
If g < define 1 — q = cosh6' and 1 + q = sinh6'. 
Then we see 

\|smh, cosh|(a;/a + 0) J 
sinh(x/a — 9) 



Vi 

{sinh, cosh, exp}(x/a) 
(Tietz potential), 

A sinh(x /a) + B cosh(x / a) 



Vi , 

{sinh, cosh}(a;/a) 
= (Eckart potential or Manning-Rosen potential as appropriate). 

(3.11.12) 

So all of these potentials are either identical to the Mobius potential, 
or special cases of the Mobius potential. To be historically accurate we 



61 



CHAPTER 3. KNOWN ANALYTIC RESULTS 



should really just call this whole collection of potentials the Eckart potential, 
as Eckhart seems to have been the first author to have given the general 
form. Unfortunately other names are now in such common use that historical 
accuracy is difficult (if not impossible) to recover. 

3.12 Discussion 

Let us summarize the results that we have obtained from this chapter. 

In this chapter, we collected many known analytic results in a form 
amenable to comparison with the general results we shall soon derive. In 
addition, we introduced the concept of quasinormal modes. We shall use 
these tools for comparing the bounds with known analytic results. More- 
over, we reproduced some of the analytically known results, and showed (or 
at least sketched) how to derive their scattering amplitudes, and so calcu- 
late quantities such as the tunnelling probabilities and quasinormal modes 
[QNM] . We did this explicitly for the delta-function potential, double-delta- 
function potential, square potential barrier, tanh potential, sech^ potential, 
asymmetric square-well potential, the Poeschl-Teller potential and its vari- 
ants, and finally the Eckart-Rosen-Morse-Poeschl-Teller potential. 

In addition, we are able to gain some deeper understanding by realizing 
that the Eckart-Rosen-Morse-Poeschl-Teller potential is generally a Mobius 
function of the variable exp(— 2x/a). Furthermore, the Morse potential is 
actually a specific limit of the Mobius potential as various parameters go to 
unity or zero. We also demonstrate that the Hulthen potential is actually a 
special case of the Manning-Rosen potential {A = 0). 

As previously discussed, we have seen that a many of the "exactly solv- 
able" potentials commonly encountered in the literature are actually the same 
quantity in disguise. For instance, we devoted that all of these potentials are 
either identical to the Mobius potential, or special cases of the Mobius po- 
tential. Moreover, we should really just call this the Eckart potential, as he 
seems to have been the first author to have given the general form. Unfor- 



62 



3.12. DISCUSSION 



tunately other names are now in such common use that historical accuracy 
is difficult (if not impossible) to recover. 

Some potentials for which the transmission probability is explicitly known 



Name 


Potential V{x) 


Transmission coefficient T = |ip 


Delta-function potential 


aS{x) 


1 

1 + {ma^/2Eh^) 


Double-delt a-function 


a{6{x - L/2) + 6{x + L/2)} 


1 


1 + 


[^cos(feL) + i(^)2sin(fcL)]' 


Square potential barrier 


\ 


Vo (for 0<x<L) 
(otherwise) 


4A;V + (fc2-gTsin^(9^) 


Tanh potential 


y+oo tanh(x/L) 


1- 


( sinh[7r(fc_oo - k^oo)L/2\ \ ^ 
i^sinh[7r(A;_oo + fc+oo)V2]y 


Sech^ potential 


Vesedi^ix/L) 


>sech2 


Asymmetric Square- well 


< 


Vi, X < a; 
V2, a < X < b; 
. V3, b<x. 


4A;ifc|A:3 
- {kl + kik^Y 



Table 3.1: This table shows some of the potentials for which exact ana- 
lytic results are known, and summarizes key properties of the transmission 
probabilities. 
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Inter-relationships between various "exactly solvable" potentials 



Name 


Potential V{x) 


Properties 


Eckart (1930) 


where ^ = — exp(2x/a) 


Eckart -i^ Rosen-Morse 


T?r\G<=»n l\/rr\T*G<=i (^Q'^0\ 
riObcIl IVlOrbe yliyO^J 


Stanh(x/d) - Csech^(a;/d) 


l-'i^/iic'/^ri 1 1 1 1 i!»T* ^ — ^ l-< i^o/nvi \ /I /~\"r»c' 

Jr^OcbCili ±ciicr v=p- rLObcIl IVlOrbc 


Poeschl-Teller (1933) 


Vb cosh^/x 

x{tanh([x — nL\/L) +tanh;u}^ 


Rosen-Morse <^ Poeschl-Teller 


Eckart-Rosen-Morse- 
-Poescm-ieller 


A + BtsM\i{x/a + e) 
+Ctanh^(x/a + (9) 


This is generally a Mobius function 
oi the variable exp(— 2a;/a) 


Mobius 




A+B exp(-2z/a) 
C+D cxp(-2a;/a) 


2 


The "best" of these equivalent forms 


Morse (1929) 




^1 — cxp(— [a; — 


lol/o))^ 


Specific limit of the Mobius potential 


Manning-Rosen (1933) 


Bcoih{x/d) — Ccosech^(a;/d) 


Obtained from Eckart by specific 
substitution 


Hulthen (1942) 


Vb 


cxp(— a;/a) 




A special case of Manning-Rosen 
(^ = 0) 


1— exp(— x/a) 




Tietz (1963) 


Vo 


sinh([a;— xo]/a) 




Obtained from the general Mobius form 


{sinh,cosh,exp}(a;/a) j 










of Eckart by appropriately choosing 
the parameters 


Hua (1990) 




( l-exp(-2a::/a) \ 
1—q exp(— 2a:/a) i 


2 


Eckart or Manning-Rosen 



Table 3.2: This table shows the inter-connections between many "exactly 
solvable" potentials. Many of these potentials are identical to each other, 
though this may not always be obvious at first glance. 
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Chapter 4 

Shabat— Zakharov systems 



4.1 Introduction 

In this chapter we shall present the general concept of the so-called "Shabat- 
Zakharov systems", (sometimes called "Zakharov-Shabat" systems). We 
shall re- write the second-order Schrodinger equation as a particular set of two 
coupled first order differential equations for which bounds can be (relatively) 
easily established. In addition, we shall introduce the idea of the probability 
current and demonstrate how to obtain the probability current density. We 
shall then use the probability current (or probability flux) to describe the 
flow of probability density. 

Moreover, we shall present an "auxiliary condition" or "gauge condition" 
that is used to relate two complex amplitudes a{x) and b{x) that we shall soon 
introduce, and to eliminate da/dx in favour of db/dx. This allows us to write 
(Pip/dx'^ in either of two equivalent forms, which is key to developing a 2 x 2 
matrix formalism. We shall represent the wave function in an inner product 
form which is the explicit general (but formal) solution to the Schrodinger 
equation. The general solution depends on three arbitrarily chosen functions 
(fi{x), A(x), and x(a;), and a path-ordered exponential matrix. 

We shall consider path ordering as an "elementary" process to derive 
the holy grail of ODE theory (complete quadrature, albeit formal, of the 
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second-order linear ODE). We shall then use the freedom to independently 
choose ip{x), A{x), and x{^) to simplify the Bogoliubov coefficients (both 
the relevant ODEs and the general bounds) as much as possible. 



4.2 Ansatz 

Consider the one-dimensional time-independent Schrodinger equation [37]- 

m 

d^ 

ij{x) + V{x) ij{x) = E ij{x). (4.2.1) 



2m dx^ 
Introduce the notation 



^2 

So we are really just trying to solve 

d2 



kixr ^ ?!Ii^IlZM. (4.2.2) 



dx"^ 

or equivalently in the time domain 

d^ 



ip{x) + k{xy ^{x) =0, (4.2.3) 



^^^m+^itym = o. (4.2.4) 

Motivated by the JWKB approximation, 

^^^eM^mx)]^^eM-^Ikix)]^ ^42.5) 
Vk{x) Vfc(x) 

the key idea is to re-write the second-order Schrodinger equation as a set 
of two coupled first-order linear differential equations (for which bounds can 
relatively easily be established). Systems of differential equations of this type 
are often referred to as Shabat-Zakharov [53] systems. A similar represen- 
tation of the Schrodinger equation is briefly discussed by Peierls [|54j and 
related representations are well-known, often being used without giving an 
explicit reference (see e.g. [HS])- However an exhaustive search has not un- 
covered prior use of the particular representation presented here, (apart, of 
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WKB (or JWKB) approximation: In general, WKB (or JWKB) 
theory is a method for estimating the solution of a differential equation 
whose highest derivative is multiplied by a small parameter e. The 
method of approximation is as follows [H]: 
For a differential equation 

e + a{x) + ... + k{x) -f- + m{x)y = 0, (4.2.6) 

assume a solution of the form of an asymptotic series expansion as 



y\x) ~ exp 



^ oo 

-5^e"S„(x) 

n=0 



(4.2.7) 



In the limit e — 0, substitution of the above ansatz into the differential 
equation and canceling out the exponential terms allows one to solve 
for an arbitrary number of terms Sn[x) in the expansion. 
The most common application of this general formalism is to the 
second-order differential equation presented in the standard form: 



dx 



2+m(x)?/ = 0. (4.2.8) 



Keeping only the first two terms in the WKB approximation yields 

exp s ±^ / a/ m[x) dx \ 
y{x) ^ ^ ^. (4.2.9) 

It is this "standard" form of the WKB approximation that we will be 
using extensively in this thesis. 



course, from [88]), nor the idea of using the representation to place bounds 
on one-dimensional scattering. 

We will start by introducing two arbitrary auxiliary functions '-pix) and 
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A{x) which may be either real or complex, though we do demand that ip'{x) ^ 
0, and then defining 

ij(x) = a(x) I— + 6 a; — — ^ -. 4.2.10 

This representation effectively seeks to use quantities resembling the "phase 
integral" wavefunctions as a basis for the true wavefunction [21]. We will 
ultimately want to interpret a{x) and h{x) as "position-dependent WKB-like 
coefficients"; in a scattering problem they can be thought of as "position- 
dependent Bogoliubov coefficients". This representation is of course ex- 
tremely highly redundant, since one complex number -^(x) has been traded 
for two complex numbers a{x) and 6(x), plus two essentially arbitrary auxil- 
iary functions ^{x) and A (a;). 

To reduce this freedom, we introduce an "auxiliary condition" (or "aux- 
iliary constraint", or "gauge condition"): 

d ^exp(.A)^ ^ d AexpHA) \ ^ ^(^^ ^^^y 2.II) 



dx \ / dx \ 

Here x{^) is yet a third arbitrary function of position. It is allowed to be 
complex, and may be zero. The original analysis, published in [88] corre- 
sponds to the special case A(x) = and = 0. Subject to this gauge 
condition, 

dV' 



iA/V?' {'^(^) 6xp(+iv5 + iA) — 6(a;) exp(— — zA)} + X (4.2.12) 
ax 

Repeated differentiation of this equation will soon lead to our desired result. 



4.3 Probability current 

We use the probability current (or probability flux) to describe the flow 
of probability density. The following equation is use to describe the flow 
of a fluid, a process which occurs all the way through physics and apphed 
mathematics [13] 

liPM = (4.3.1) 
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Definition (Zakharov— Shabat): We shall now turn to the study 
of the non self adjoint Zakharov-Shabat problem. In a rather general 
context, this may be viewed as the first order system of differential 
equations in 71 defined by [7]: 



h 




u[{x, z) 

U2{X,Z) 

—iz A{x) exp{iS{x) / h) 
exp(— ■jS'(x)//i) iz 



where z is a complex eigenvalue parameter, and the prime denotes 
differentiation with respect to x. This system appears (for instance) 
when one wants to solve the non linear initial value problem given by 

ihdtil) + —dl^lj + \ij\^ij = Q, (4.2.14) 

V'|t=o = A{x) ex^{iS{x)/h) , (4.2.15) 
with the inverse scattering method 

The same sort of system of ODEs will also arise in our particular 
problem: We are of course interested in the somewhat simpler problem 
of solving the linear Schrodinger equation. 



The above equation tells us that the rate of change in density is equal to the 
negative of the difference between the amount of "stuff" (be it water, air, or 
"probability") flowing into the point and the amount flowing out. Now in 
the context of the Schrodinger equation, we can write the probability density 
in the following form 

p(x,t) = |7/>(x,t)|2 = ^>. (4.3.2) 
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In quantum mechanics we can certainly describe the movement of particles, 
however we also have an additional difficulty because our particles are not 
classical. Therefore we can only talk about the probability of the particle 
being at a certain place in time, now we can talk about a probability current 
or flux. Consider the time-dependent Schrodinger Equation and its complex 
conjugate: 



Now let us differentiate the probability density with respect to time 

h fd'^ij* , d'^i)\ vij*i/j vij*ip 

- TT^ ^ + 



2m i \ dx"^ dx^ ) ih ih ' 

' (4.3.5) 



Apply 



dx \ 2m i\ dx dx 



|p(^>^) + j^^(^,^) = 0. (4.3.6) 



So we obtain the probability current density 



We can re-write this as 

^(-■'> = ^-(^*^)- 

Here (h/m) is just a normalization (that is often set 1 for convenience). 
There is nothing really important in this normalization (unless we want to 
calculate experimental numbers), so we might as well set 

Jix,t)^lm(^r^y (4.3.9) 
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Now at this stage (/^(.t) and x{^) completely arbitrary possibly complex 
functions subject only to the constraint 7^ 0. For future use, compute the 
probability current using our WKB-based ansatz in terms of a{x) and b{x): 



= Im < ip 



i^/(p'{a{x) exp{+i(p) — b{x) exp{—i(p)} + xi^ 



(4.3.10) 




Re < A —[a{x) exp{+iip) — b{x) exp{—i(p)] 



[a{xy exp{-i(p*) + b{x)* exp{+iip*)] } + Im{x} V'*^, (4.3.11) 



Re 




[|a|2Re{e+'(*'-^*)} - |6|'Re{e- 



}] 




+lm\ Wm{a6*e^('^+'^*)} +Im{x}^*V', (4.3.12) 



= Re 




a|2Re{e+2^'"(^)} - \b\^Re{e-^'^^^^] 



'{v^} 



+Im<^ Um{a6*e2^^*^('^)} + Im{x}V'*'0■ 



(4.3.13) 



4.4 SDE as a first order system 

We now re-write the Schrodinger equation in terms of two coupled first- 
order differential equations for these position-dependent WKB/Bogoliubov 
coefficients a{x) and b{x). To do this we evaluate d'^ip/dx'^ in two different 
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ways, making repeated use of the gauge condition. First 

d'^i/j d 



dx 



= {ae+''^+'^ - be-''^-'^} + x V , (4-4.1) 



'\2 



d f ae 



dx 



lA 



d /be 



e ' ^ - 



-iA 



dx \ 



-i(p 



a e 



iA 



dx \ 



But then 



dx'2 



a/v^' dx 



+X 



ivV{«(a;)e+^'^+^^-6(a;)e"*'^"'^} • (4-4.4) 



So finally 



d^^ 
dx"^ 



[x^ + x'-(^')^]^ + ^^e--^^- 



(4.4.2) 



(4.4.3) 



(4.4.5) 

Now use the gauge condition to eliminate da/dx in favour of d6/da;. This 
permits us to write d^il^/dx"^ in either of the two equivalent forms 



dx^ 



-2v/^A'e-^'^-^^6 + i [if" + 2x¥^'] a- 
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and 



dx ^/^p' 



dx'^ 



-iip—iA 



(4.4.7) 



Now insert these formulae into the Schrodinger equation written in the form 

d^^ 



dx^ 

to deduce the first-order system: 



+ k{xY = 0, 



(4.4.8) 



da 
dx 



+ {y + 2x^'] + I [k\x) + + X' - iv'?]) e-'^^-'^^ h 1,(4.4.9) 



£ = +^ ^ (b" + 2x^1 - ^ [k\x) + + X' - (^')'] ) 



-z[A;2(x) + x' + x'-(^')'-2^'A'] h) . 



(4.4.10) 



It is easy to verify that this first-order system is compatible with the "gauge 



condition" (4.2.11), and that by iterating the system twice (subject to this 
gauge condition) one recovers exactly the original Schrodinger equation. 
These equations hold for arbitrary ^{x)., A(x), and x(x), real or complex. 

When written in 2 x 2 matrix form, these equations exhibit a deep connec- 
tion with the transfer matrix formalism [56]. Let us define quantities pi{x) 
and P2(a^), not necessarily real, as 



pi = + 2x<^'; P2 = k\x) + x' + X- {vY- 



(4.4.11) 



We can then re-write the Shabat-Zakharov system in matrix form as 



d 

dx 



1 



i[p2 - 20' A'] {pi + ip2} exp(-2i(p - 2iA) 

{pi - ip2] exp(+2i(^ + 2zA) -i[p2 - 20' A] 

(4.4.12) 
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This has the formal solution 

a{xf) 
b{xf) 



a{xi 
b{xi 



(4.4.13) 



in terms of a generalized position-dependent "transfer matrix'' 
E{xf,Xi) = 



V exp 



1 



{Pl-^P2} e+2*^+2''^ -z[p2-2^'A'] 



dx 
(4.4.14) 



where the symbol V denotes "path ordering" 



Path ordering: In theoretical physics, path-ordering is the proce- 
dure (or meta-operator V) of ordering a product of many operators 
according to the value of one chosen parameter [29] : 

P[0i(cri)02(0-2)...0„((T„)] = 0p(i)(0-p(i))0p(2)((Tp(2))...0p( 

(4.4.15) 

Here p : {1,2, ...,n} {1,2, ...,n} is a permutation that orders the 
parameters: crp(i) < (7p(2) < ••• < 0"p(„). For instance 

P[Oi(4)02(2)03(3)04(l)] := 04(l)03(3)02(2)Oi(4). (4.4.16) 



Equivalently if we were to be working in the time domain we would have 
Eitj,U) = 



T exp 



{pi - tp2} e+2«^+2^A _ A'] 



dt 
(4.4.17) 



where T would now be the well-known "time ordering" operator (more usu- 
ally encountered in quantum field theory) and we would now define 

pi = ip" + 2x^'; P2 = uj\t) + - X' - {^')\ (4.4.18) 



74 



4.5. BOUNDING THE COEFFICIENTS A{X) AND B{X) 



with If, A, and x iiow being arbitrary functions of t rather than x. 
We can now write the wave function in inner product form 

1 



ijj{x) 



exp(+i(p + iA); exp(—i(p — iA) 



a[Xi) 
b{xi) 



(4.4.19) 



to yield a formal but completely general solution for the Schrodinger equation 



jp{x) 



1 



exp{+i(f + iA); exp(— i</7 — iA) 



E{x,Xo) 



a{xo) 
b{xo) 



(4.4.20) 



Explicitly 






ip(x) = 


1 


e> 








T'exp 




1 





exp(+i</7 + iA); exp(— i</7 + iA) 



(4.4.21) 





a{xQ) 




b{xo)_ 



This is the explicit general solution to the Schrodinger equation. It depends 
on the three arbitrarily chosen functions <f{x), A{x), and x(^)) and a path- 
ordered exponential matrix. If you consider path ordering to be an "elemen- 
tary" process, then this is indeed the holy grail of ODE theory (complete 
quadrature, albeit formal, of the second-order linear ODE). 



4.5 Bounding the coefficients a{x) and b{x) 



Prom 

da 
dx 



db 

dx 



+ ^,^t[e{x)+x' + x'-{vr-2^'^'] a 

+ (b" + 2x^'] + I [k\x) +x' + X- (^T]) e-^^--^^^ ftj, 

(4.5.1) 



(4.5.2) 
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we see 

+ {[ip" + 2x^'] + I [k\x) + + x' - e"''""''^ 

(4.5.3) 

Therefore (now assuming for the rest of this section that A, x are all real 
and (f' > 0) 

, da da* _ Re{{[(p" + 2x^)'\ + % \k^{x) + x^ + x' - {v'f\) e-^'^-^'^ a* b} 
dx dx (fi' 

(4.5.4) 

This implies 

d|ap _ Re{{[ip" + 2x(p'] + i [k\x) + + x' - i^T]) e"^^^"^^^ a* b} 
dx ip' 

(4.5.5) 

But Bje{A) < \A\, so we have 

d|aP W + 2x^'] + ^ [k\x) + x' + x' - {^'f]) e"''^"'^^ h\ 
ax </?' 

implying 

^, ,d|a| ^ |b-- + 2xy^1 + z[fe^(x) + x^ + X^-((^0^]| |a| |6| 

2|a|— < . (4.5.7) 



Thus 



d|a| ^ I y' + 2x^']-ri[k\x) + x^-rx!-{^'?] I 1^1 (4 5 



dx - 2lp' 

and so 



Now if (as per our current assumption) 93, A, x are all real, then it is easy 
to check that 

J = Im(^V) = \a? - l&r, (4.5.10) 
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so current conservation implies 

|a|2-|6|2 = l. (4.5.11) 

Ultimately, it is this equation that allows us to interpret a{x) and b{x) as 
"position dependent Bogohubov coefficients" . 

In view of the relation between a{x) and b{x) we have \b\ — ^/\a^~^^, so 
that we can deduce 

d|a| ^ [ip" + 2xipf + mx) +x'+x'- i^'rf ^- , 

< ^ — Vl«l - 1- (4.5.12) 



dx - 2^' 
But this inequality can now be integrated. For convenience let us define 

Jl^- + ixf'Y + + x= + X' - (v'ff 

" = V ■ 



Then 



But now 



so that 



d|a| 
dx 

1 d|a 



V|a|2-1 dx 



< -d - 1- (4.5.14) 

dx< j 'd dx, (4.5.15) 



{cosh-^|a|}^^ < /" ^-(^dx. (4.5.16) 

Now apply the boundary conditions: as — > — oo we have chosen to set 
things up so that we have a pure transmitted wave, so |&(— oo)| = and 
|a(— oo)| = 1. On the other hand as — > +oo we have chosen to set things 
up so that a{x) and b{x) tend to a and (3, the Bogohubov coefficients we 
are interested in calculating. Thus taking the double limit Xi — > — oo and 
X/ — > +00 we see: 

r+oo 

cosh"^ \a\< 1^ dx. (4.5.17) 



— oo 



That is 

"+00 



a\ < cosh < i3 dx}. (4.5.18) 
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This is the central result of this thesis — it can be modified and rearranged 
in a number of ways, and related inequalities can be derived under slightly 
different hypotheses, but all the applications we are interested in will reduce 
in one way or another to an application of this inequality or one of its close 
variants. 

For notational convenience, we often find it is useful to adopt the short- 
hand 

/r+oo 
= J , (4.5.19) 

since then 

|a| < cosh \ (b ^dx\ . (4.5.20) 



From the normalization condition (4.5.11) we immediately deduce 



|/3| < sinh |y ?9 dxj . (4.5.21) 

When translated into equivalent statements about transmission an refiection 
probabilities, we find 

T > sech^ \ (p^dxl , (4.5.22) 



and 

< tanh^ |y ^9 dxj . (4.5.23) 

Now one of the the points of the exercise (and of this thesis) is to use 
the freedom to independently choose (p{x), A(x), and x{^) to simplify life 
as much as possible. Here are a few special cases, chosen for their simplicity 
and the lessons they teach us. 

4.5.1 Case: A' = p2/{2(p') 

No one can prevent us from choosing 

A' = 0. (4,5,24) 
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that is 

A'=*!M±X!±X:^_ (4.5.25) 

which imphes 

^^jm±^^±^^^, (4.5.26) 
Doing that greatly simphfies hfe since now the system of ODEs becomes 

da 1 

dx 2ip' 



(4.5.27) 

A 

(4.5.28) 



That is 



d 

dx 



1 



{pi + ip2} exp(— 2i(/7 — 2iA) 

{pi — ip2} exp(+2i(/7 + 2iA) 

(4.5.29) 

If you now let (p and x be real, then the matrix above is Hermitian; unfor- 
tunately this does not tell us all that much about the evolution operator. 
(If we had i times a Hermitian operator appearing above, then the evolution 
operator would have been unitary.) Using a specialization of the previous 
argument, we can easily deduce 



d|a| _ \/[^" + 2x^? + mx) + x' + X'-{^'rf ^^^g^^ 



dt 2\ip'\ 
We also have the same constraint 

|a|2-|6|2 = l, (4.5.31) 

and so deduce 

d|a| ^ ^J[^" + 2x^f + [k^jx) + x' + X'- {^Tf r-r^ 

dt 2|(^'| ■ ^ ■ ■ ■> 
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Defining 
we have 

d|a| 



1? V|a|2-1. (4.5.34) 

This inequahty can now be integrated in the manner discussed above, though 
doing so gives us no additional information. 



4.5.2 Case: A = —(p 

No one can prevent us from choosing 

A{x) = -(p{x), (4.5.35) 

in which case 

+ {[<f" + 2x<f']+i[e{x)+x' + x'-{<fr]) (4.5.36) 

£ = +^^{[^" + '^Xf']-imx)+x' + x'-{^r]) a 

-^[k^^x)+x' + X+{^'?] &|. (4.5.37) 
The comphcated phase structure has gone away, and we have 

a*^ = +^,l^^[k^x)+x' + x'+{^r] a* a 

+ {[<p" + 2x<p'] + I [k'ix) + x' + x'- i<p'f]) &|, (4.5.38) 
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whence 



,da da- _ Re{{y" + 2x^'] + i [k\x) + + ^ - {^'f]) 

(X '\ Qj — 

dx dx 2ip' 

(4.5.39) 

Perhaps more to the point, we can derive ODEs for the sums and differences: 



(4.5.40) 



d(a + h) 


1 


da; 


V 


d(a - h) 


1 


dx 





that is 

d[i(a + h)] 
dx 

d[i(a - h)] 
dx 



+ -\[^" + 2xv']{a + h) + 2{ip'ma-h)] 



(4.5.41) 



(4.5.42) 



+ -| [k\x) + x' + x!] (a + 6) - [if" + 2xv'] [i{a - 6)] j. 

(4.5.43) 



While this system of ODEs is somewhat simpler than those derived above, 
we have not been able to extract any significant improvement on or previous 



results in equations (4.5.13) and (4.5.18). 



4.5.3 Case: A = 

We include this case for historical reasons, as it was the first generalization we 
obtained of the original result published in [88]. The derivation is somewhat 
simpler than the full A 7^ discussion presented above, and the ultimate 
bound we extract is no weaker. 

We introduce an arbitrary auxiliary function ^{x) which may be either 
real or complex, however we demand that ^'{x) 7^ 0, and then define 



^(x) = a(x)^^P^ + Ka:)^^^^, 



(4.5.44) 
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Again, to trim down the number of degrees of freedom it is useful to impose 
what can be thought of as a "gauge condition" ( "auxihary condition" ) 



d 

da; 



dx 



- e-'^ = x(x)^(x). 



(4.5.45) 



Here x is some arbitrary function of position. The original analysis in 
corresponds to the special case x{^) = 0- Subject to this gauge condition, 



dV' 
da; 



i^/(f'{a{x) exp{+i(f) — b{x) exp(— + x4'- 



(4.5.46) 



Now the Schrodinger equation can be rewritten in terms of two coupled 
first-order differential equations for these position-dependent Bogoliubov co- 
efficients. We have to calculate d'^ip/dx'^ making repeated use of the gauge 
condition: 



dV 
da;^ 



± (^^^{ae+'^ - be-'^} + X^P 



(4.5.47) 



/\2 



, r d / a 



dx \^/^ 



Htfi 



vV 

/2 C A 

={ae^'^ + be-'^} + iip' {2 



dx VVV 



-tip 



dx 



-lip 



+X 



2i(p' da 



Hp 



i^be-'^ - i^'x'ip + XV 



ip' dx ^ip' 



\r + x - + ^ 



Lp' dx 



(4.5.48) 



(4.5.49) 

(4.5.50) 
(4.5.51) 



Now use the gauge condition to eliminate dajdx in favour of d6/da;. Finally, 
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this permits us to write d^ip/dx"^ in either of the two equivalent forms 



dx 



dx 



(4.5.52) 



-tip 



dx V yj' 



(4.5.53) 



Now insert these formulae into the Schrodinger equation written in the form 



to deduce 



dx"^ 



-k{xYip 



2m[E - V{x)] 



(4.5.54) 



da 
dx 



db 
dx 



+t [k'ix) + x' + X - i^T] (« + bexp{-2tif)) )■ , 



+ ^{ [^" + 2xip']aexp{+2iip) 



(4.5.55) 



-t [e{x) + x +X'- iv'n {b + aexp(+2zv?)) 



(4.5.56) 



It is (again) easy to verify that this first-order system is compatible with the 



"gauge condition" (4.5.45), and that by iterating the system twice (subject to 



this gauge condition) one recovers exactly the orginal Schrodinger equation. 
These equations hold for arbitrary ip and x^ ^^^al or complex, and when 
written in matrix form, exhibit a deep connection with the transfer matrix 
formalism |56] . Let us define quantities pi (x) and p2 (x) , not necessarily real. 



as 



pi = ip" + 2xf'; p2 = k {x) + x +X - if') 



i\2 



(4.5.57) 
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We can then re-write the Shabat-Zakharov system in matrix form as 



d 


a 


1 


ip2 


{pi + ip2}exp(-2zv9) 




a 


dx 


b 


~ V 


{Pi - ip2} exp{+2i(f) 


-ip2 




b 



This has the formal solution 

a{xf) 
b{xf) 



E{xf, 



Xi 



a[Xi 
b{xi 



(4.5.58) 



(4.5.59) 



in terms of a generalized position-dependent "transfer matrix" [56 
E{xf,Xi) = 



V exp 



ip2 {pi + ip2} exp(-2iv5) 

{Pi - ip2] exp(+2i(/)) -ip2 



da; 
(4.5.60) 



where the symbol V denotes "path ordering". (See the boxed text earlier 
in this chapter for details.) Now we can write the wave function in inner 
product form 



ip^x) 



1^' 



exp(+iy9); exp(— Z(y9) 



a\Xi 
b{xi 



and where 

Therefore 

'?/'(x) = 



a{x) 
b{x) 



E{x, xq) 



a{xo) 
b{xo) 



(4.5.61; 



(4.5.62) 



exp{+iip)] exp{—iLp) E{x,xo 



a{xo) 
b{xo) 



(4.5.63) 



to yield a formal but completely general solution for the Schrodinger equation 



(p'ip{x) = exp{+i(f); exp{—iif) E{x,Xo) 



a{xo) 
b{xo) 



(4.5.64) 
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exp{+i(f); exp{—iip) 



ip2 

{pi - ip2}exp(+2iv9) 



{pi + iip2}exp(-2zv2) 



(4.5.65) 



This is the exphcit general solution to the Schrodinger equation. It depends 
on the two arbitrarily chosen functions '^[x) and ^ind a path-ordered 

exponential matrix. If you consider path ordering to be an "elementary" pro- 
cess, then this is (again) the holy grail of ODE theory (complete quadrature, 
albeit formal, of the second-order linear ODE). 

The development of bounds automatically follows as in the previous dis- 



cussion, and we can without further calculation assert that equations (4.5.13) 
and (4.5.18) hold in this situation as well. 



4.6 Discussion 



There are several ways to derive a number of rigourous bounds on transmis- 
sion probabilities (and reflection probabilities and Bogoliubov coefficients) 
for one-dimensional scattering problems. The derivation of these bounds 
generally proceeds by rewriting the Schrodinger equation in terms of some 
equivalent system of first-order equations, and then analytically bounding 
the growth of certain quantities related to the net fiux of particles as one 
sweeps across the potential. In this chapter we obtained a number of sig- 
nificant bounds, considerably stronger than those in [88], of both theoretical 
and practical interest. 

Even though the calculations we have presented are sometimes somewhat 
tedious, we feel however, they are more than worth the effort — since there 
is a fundamental lesson to be learnt from them. Technically, we demon- 
strated that the Schrodinger equation can be written as a Shabat-Zakharov 
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system, which can then be re-written in 2 x 2 matrix form. We rearranged 
this formation in terms of a generahzed position-dependent "transfer ma- 
trix" involving the symbol V which denotes "path ordering". Therefore the 
wavefunction t/;(x) can be written in inner product form. This is the explicit 
general solution to the Schrodinger equation. It depends on the three arbi- 
trarily chosen functions (p{x), A{x), and x{x) and a path-ordered exponential 
matrix. If one considers path ordering to be an "elementary" process, then 
this is the holy grail of ODE theory (complete quadrature, albeit formal, of 
the second-order linear ODE). We have seen that it is often convenient to use 
the freedom to independently choose (p{x), A{x), and xi^) to simplify life 
as much as possible. Furthermore, we have considered a few special cases. 
For instance, case A' = P2/ i^ip'), case A = ip, and case A = 0. The bounds 
that we have derived on the Bogoliubov coefficients a and (3, and on the 
transmission and reflection probabilities T and R, are the key results of this 
thesis — the next few chapters will be devoted to developing several variants 
of these bounds, developing independent proofs that might ultimately lead 
to new bounds, and developing various applications of these bounds. 
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Chapter 5 

First derivation of the bounds 



5.1 Introduction 

In this chapter we shall review the analysis of [HE], developing various tech- 
niques for estimating the scattering properties. We shall review and briefly 
describe some very general bounds for reflection and transmission coefficients 
for one-dimensional potential scattering, and then indicate how the results 
of this thesis extend and expand on the earlier results. Equivalently, these 
results may be phrased as general bounds on the Bogoliubov coefficients, or 
statements about the transfer matrix [88] . 

Finally, we shall re-demonstrate the use of Shabat-Zakharov system of 
ODEs (now in a greatly simplified context) to derive a first elementary bound 
on the transmission, reflection, and Bogoliubov coefficients. 



5.2 Shabat— Zakharov systems 

Consider the one-dimensional time-independent Schrodinger equation [37]- 



2m dx' 



;iIj{x) + V{x) ipix) = E ipix). (5.2.1) 
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If the potential asymptotes to a constant, 

V{x ±00) ^ V±oo, (5.2.2) 

then in each of the two asymptotic regions there are two independent solu- 
tions to the Schrodinger equation 

ij{±i; ±00; x) ^ -. (5.2.3) 

Here the ±z distinguishes right-moving modes e^*'^^' from left-moving modes 
g-«fcx^ while the ±00 specifies which of the asymptotic regions we are in. 
Furthermore 

^ ^?^(^E^. (5.2.4) 

To even begin to set up a scattering problem the minimum requirements are 
that the potential asymptote to some constant, and this assumption will be 
made henceforth. 

The so-called Jost solutions [52j are exact solutions J7± (x) of the Schrodinger 
equation that satisfy 

Mx - +00) - -M+^h^, (5.2.11) 

Y rC-00 

and 

J+{x -00) a h p , , (5.2.13) 

V "^-00 V "^-00 

:r_(^ - +00) - «*!^P(Z£±^ + ^.exp(+£wO_ ^^_2.14) 

Here a and /3 are the (right-moving) Bogoliubov coefficients, which are re- 
lated to the (right-moving) reflection and transmission amplitudes by 

1 , , 

r=-; t=-. (5.2.15) 

a a 
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Definition ( Jost solutions) : If we were working in three dimensions 
with radial symmetry then we would study the solution of the radial 
Schrodinger equation |110] 

+ _ v{r)) f{r, k) = 0, (5.2.5) 

with the asymptotic conditions 

lim /(r ± k) exp{±ikr) = 1, (5.2.6) 

the so-called "Jost solutions". In one dimension we would consider a 
formally identical equation for the full wave-function, now with bound- 
ary conditions at both x +oo and x — oo. In addition, we now 
introduce the function 

E{r, z) = exp{±ikr)f{r, ±k). (5.2.7) 

which is also called a "Jost solution" , and which satisfies the differen- 
tial equation 

d^E(r,z) dE(r,z) , , , , , , 

- - - ^(0 ^(-, = 0, (5.2.8) 

with the following condition at infinity 



lim E{r,z) = 1, (5.2.9) 

7 >+00 

where 

z = ±2ik. (5.2.10) 



These conventions correspond to an incoming flux of right-moving particles 
(incident from the left) being partially transmitted and partially scattered. 
The left-moving Bogoliubov coefficients are just the complex conjugates of 
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the right-moving coefficients, however it should be borne in mind that the 
phases of j3 and (3* are physically meaningless in that they can be arbitrarily 
changed simply by moving the origin of coordinates. The phases of a and 
a* on the other hand do contain real physical information. In this chapter 

Phase integral technique: We shall be concerned with the solution 
of the differential equation [30], 131] 

^ + R{z)^Ij = 0, (5.2.16) 

where R{z) is an analytic function of z. This differential equation 
which could possibly result from separation of variables, describes 
large classes of important problems in various fields of physics, not 
only in quantum mechanics. 

The phase-integral functions, in terms of which the solution will be 
expressed, are of the general form 

/i(z) = q-'/\z)exp(^ + zj q{z) dz^ (5.2.17) 

f2{z) = q-'/\z) exp ( - I [ q{z)dz]. (5.2.18) 



These phase integral functions are very closely related to the WKB ap- 
proximation, and the lowest order approximation to the phase integral 
is essentially just the second-order WKB approximation. 



we will derive some very general bounds on |a| and \(3\, which also lead to 
general bounds on the reflection and transmission probabilities 

i?=|r|2; r=|tp. (5.2.19) 

The key idea is to re-write the second-order Schrodinger equation as a par- 
ticular type of Shabat-Zakharov [53] system: a particular set of two coupled 
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first-order differential equations for wliicli bounds can be easily established. 
A similar representation of the Schrodinger equation is briefly discussed by 
Peierls [51] and related representations are well-known, often being used with- 
out giving an explicit reference (see e.g. [55]). However an exhaustive search 
has not uncovered prior use of the particular representation of this chapter, 
nor the idea of using the representation to place bounds on one- dimensional 
scattering. 



Scattering matrix: In physics, the scattering matrix (or S-matrix) 
relates the initial state and the final state for an interaction of parti- 
cles. It is used in quantum mechanics, scattering theory and quantum 
field theory. In addition, it can also denote an infinite-dimensional 
matrix (or operator) that expresses the state of a scattering system 
consisting of waves or particles or both in the far future in terms of 
its state in the remote past; also called the S-matrix. In the case of 
electromagnetic (or acoustic) waves, it connects the intensity, phase, 
and polarization of the outgoing waves in the far field at various an- 
gles to the direction and polarization of the beam pointed toward an 
obstacle l33l l3^ 



Definition (Auxiliary functions): They are not a rigorously de- 
fined type of function. In contrast, these functions are either explic- 
itly constructed, or at least shown to exist, and are used to provide 
a formal solution to some assumed hypothesis, or otherwise prove the 
result in question. Devising an auxiliary function during the course of 
a proof in order to prove the result is not a technique exclusive to any 
particular branch of mathematics |103] . 
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Definition (Gauge conditions or Gauge fixing): Indicates a 
mathematical scheme for coping with redundant degrees of freedom 
in field variables, most commonly in the physics of gauge theories. In 
Abelian or non-Abelian gauge theories one represents each physically 
distinct configuration of the system as an equivalence class of detailed 
local field configurations. Any two detailed configurations in the same 
equivalence group are related by a gauge transformation, equivalent 
to a shear along unphysical axes in configuration space |104j . 
In the context of the current thesis things are somewhat simpler. 
There is a redundancy, which is eliminated by our "gauge condition" , 
but one does not have to deal with the full complexity of an actual 
"gauge theory" in the sense of particle physics. 



We start by introducing an arbitrary auxiliary function ip{x) which may 
be either real or complex, though we do demand that ip'{x) ^ 0, and then 
defining 

^Pix) = a{xf-^^^ + Kxf-^^^. (5.2.20) 

This representation effectively seeks to use quantities resembling the "phase 
integral" wavefunctions as a basis for the true wavefunction [72j. This rep- 
resentation is of course highly redundant, since one complex number iIj{x) 
has been traded for two complex numbers a{x) and b{x) plus an essentially 
arbitrary auxiliary function ip{x). In order for this representation to be most 
useful it is best to arrange things so that a{x) and b{x) asymptote to con- 
stants at spatial infinity, which we shall soon see implies that we should pick 
the auxiliary function to satisfy 

(p'{x) k±oo as X ±oo. (5.2.21) 

To trim down the number of degrees of freedom it is useful to impose a "gauge 
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condition" 



d 

dx 



d 

dx 



e-''^ = 0. 



(5.2.22) 



Subject to this gauge condition, 



dV' 
dx 



— i^J~{p' {a{x) exp(+i(^) — exp(— . 



(5.2.23) 



We now re-write the Schrodinger equation in terms of two coupled first-order 
differential equations for these position-dependent Bogoliubov coefficients. 
To do this note that 



d^v^ _ d ^' 

dx^ dx 



dx 



+ 



if' 
2iip' da 



'cp' dx sji^' 



2Vd6 ^" 
V dx 



dx 



(5.2.24) 



(5.2.25) 



(5.2.26) 



(5.2.27) 



(The last two relations use the "gauge condition" .) Now insert these formulae 
into the Schrodinger equation written in the form 



(5.2.28) 
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to deduce 



da 1 1 // 7 

- = +— 6exp(-2zv.) 



+i [k'^ix) - {if'y] {a + bexp{-2iip)) |, (5.2.29) 

db if,/ , ^ 

- = a exp(+2z^) 

-i [k^{x) - i^'f] (6 + aexp(+2i(^)) |. (5.2.30) 
It is easy to verify that this first-order system is compatible with the "gauge 



condition" (5.2.22), and that by iterating the system twice (subject to this 
gauge condition) one recovers exactly the original Schrodinger equation. 
These equations hold for arbitrary ip, real or complex, and when written 
in matrix form, exhibit a deep connection with the transfer matrix formal- 
ism 



5.3 Bounds 

To obtain our bounds on the Bogoliubov coefficients we start by restrict- 
ing attention to the case that (p{x) is a real function of x. (Since ip is an 
essentially arbitrary auxiliary function this is not a particularly restrictive 
condition). Under this assumption the probability current is 

^=Im{^*g} = {|ar-|6r}. (5.3.1) 

Now at X ~ +00 the wavef unction is purely right-moving and normalized to 
1, because we are considering one-dimensional Jost solutions Then for 
all X we have a conserved quantity 

\af-\b\^ = l. (5.3.2) 
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It is this result that makes it useful to interpret a{x) and b{x) as position- 
dependent Bogoliubov coefficients relative to the auxiliary function ip{x). Now 
use the fact that 

d\a\ 1 



dx 2 1 a I 



da da* 

a'- h a—— 

ax ax 



(5.3.3) 



and use equation (5.2.29) to obtain 
dial 1 1 



dx 



-[(p" [a*bexp{-2iip) + ab* exp{+2i(p)] 
2\a\ lip \ 

- [}{>'f\ [a*6exp(-2i(y9) - ab* exp(+2z(y9)] 



(5.3.4) 



That is 



d|a| 
dx 



1 1 

2H V 



yRe [(^" + z[A;2- (^')2]] [a*6exp(-2z(^)] 



(5.3.5) 



The right hand side can now be bounded from above, by systematically using 
Re(^ B) < \A\ \B\. This leads to 



d|a| ^ + 



dx 



2|(/?'| 



(5.3.6) 



It is essential that be real to have | exp(— 2zy9)| = 1 which is the other key 
ingredient above. Now define the non-negative quantity 



^ = ^[ip{x),k{x)] = 



y/(v.-)^ + [P(x)-(yO 
2|^'| 



and use the conservation law (5.3.2) to write 



d|a| 
dx 



1. 



Integrate this inequality 



(cosh ^ |a| I T'^ < 
95 



{} dx. 



(5.3.7) 



(5.3.8) 



(5.3.9) 



CHAPTER 5. FIRST DERIVATION OF THE BOUNDS 



Taking limits and Xf ^ +00 

/+00 

^dx. (5.3.10) 
-00 

That is ^ 

|a| < cosh ( / ^dx] . (5.3.11) 



00 



Which automatically implies 

\i3\<smh(^J i^dx^. (5.3.12) 

Since this result holds for all real choices of the auxiliary function ip{x), 
(subject only to (f' ^ and (f' k± ±00), it encodes an enormously 

wide class of bounds on the Bogoliubov coefficients. When translated to 
reflection and transmission coefficients the equivalent statements are 

T>sech.'^(^J i9dx^, (5.3.13) 

and ^ 

R<tanh^(^j i3 dx^ . (5.3.14) 

We shall soon turn this general result into more speciflc theorems in chapterlGl 



5.4 Transfer matrix representation 



The system of equations ( 5.2.29 )-( 5.2.30) can also be written in matrix form. 
It is convenient to deflne 



p^^" + z[k\x)-{^y]. 



(5.4.6) 



Then 



d 


a 


1 


zlm[|0] pexp{—2i{p) 




a 


dx 


b 


~ V 


p* exp{+2iip) —ilm[p] 




b 



(5.4.7) 
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Transfer matrix formalism: This is most commonly used in optics 
and acoustics to analyze the propagation of electromagnetic or acous- 
tic waves through a stratified (layered) medium. In those situations, 
the stack layers are normal to the z axis and the field within one layer 
can be represented as the superposition of a left and right traveling 
wave with wave number k [32], 

E{z) = Er exp{ikz) + Ei exp{-ikz), (5.4.1) 

We represent the field as the vector {E{z), F{z)), where 

F{z) = ikEr exp{ikz) — ikEi exp{—ikz). (5.4.2) 

The propagation over a distance L into the positive z direction is 
described by the matrix 



M 



cos{kL) I sin(A;L) 
-A;sin(A;L) cos{kL) 



(5.4.3) 



and 



" E{z + L) ' 


= M X 


' E{z) ' 


F{z + L) 




F{z) _ 



(5.4.4) 

The system transfer matrix is then built up by repeated multiplication 

Ms = Mj^ X ■ ■ ■ X M2 X Ml. (5.4.5) 

A very similar formalism applies to the one- dimensional Schrodinger 
equation, and the "path ordered" exponentials discussed in the body 
of the thesis are effectively transfer matrices built out of an infinite 
number of infinitesimally small steps. 
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This has the formal solution 



a{xf) 
b{xf) 



E{xf,Xi 



a{xi) 
h{xi) 



(5.4.8) 



in terms of a generalized position-dependent "transfer matrix" [73] 

E{xf,Xi) = 
V exp 



f 1 



zlm[p] pexp{—2i(p) 
p* exp{+2iip) — zlm[p] 





a 


dx^ 




b 





(5.4.9) 



where the symbol V denotes "path ordering". In particular, if we take 
Xi — i> — oo and xj ^ +oo we obtain a formal but exact expression for the 
Bogoliubov coefficients 



a (3* 
(3 a* 



E{oo, — oo) 



V exp 



oo 

-oo V 



p* exp{+2i{p) —ilm.[p] 



dx 



(5.4.10) 



The matrix E is not unitary, though it does have determinant 1. It is in fact 
an element of the group SU{1, 1). Taking 



+1 
-1 



(5.4.11) 



then (cr^)^ = +/, and defining E'^ = {E*Y , it is easy to see 



(5.4.12) 



This is the analog of the invariance of the Minkowski metric for Lorentz 
transformations in 5*0(3, 1). 
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Complex structure: There is a minor (non-propagating) error in 
the analysis of [88]. If we define the "complex structure" J by 



J 



1 
-1 



(5.4.13) 



then = —I and it is claimed in [88j that 

= JEJ. 



(5.4.14) 



This assertion is simply an error. Fortunately this is a side issue, and 
the error does not affect the rest of the argument. 



5.5 Discussion 

In this chapter we have re-cast, re-analyzed, and described the first derivation 
of scattering bounds as presented in [88]. The formalism as developed here 
works in terms of one free function ip{x). In other parts of this thesis we have 
established generalized bounds; some in terms of two arbitrary functions ip{x) 
and x{^)y cind some in terms of three arbitrary functions if{x), A(x), and 
x(x). The derivation of the present chapter is noteworthy because of its 
brevity and simplicity — and this chapter has acted as a "seed" , suggesting 
and hinting at generalizations that have ultimately become the content of 
the previous chapter. Of course, the "simple" calculation reported in this 
chapter is also the seed for the various published journal articles that have 
already arisen from this thesis, articles which are displayed in the appendices 
to the thesis. 
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Chapter 6 

Bounds: Special cases 



In this chapter we shall deal with some specific cases of these general bounds 
and develop a number of interesting specializations. We shall collect together 
a large number of results that otherwise appear quite unrelated, including re- 
flection above and below the barrier. We have divided the special case bounds 
we consider into five sub-cases: Special cases 1-4, and "future directions". 



6.1 Bounds: Special case 1 

We now reproduce the bounds of special case 1 in jSH] . 

Suppose now that the potential satisfies V"+oo = ^-oo- Also, choose the 
phase function ip{x) to be (y9 = k^oX and take x = 0. We also require k^o 7^ 0, 
that is E > V±oo. This is the special case discussed in a different context by 
Peierls [51]. Then the evolution equations simplify tremendously. 

We consider the quantity 



2\ip'\ 

which represents the generalization that we derived earlier in this thesis. 



in equation (4.5.13), of the bound reported in [88]. When ip = kooX, then 
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ip" = 0, and if we additionally choose x = O5 then this simplifies to 



= ^^''^"1, = (6.1.2) 

I "'oo I ^ I "^00 I 

Inserting this into our general bound now reproduces case 1 of |88] as desired. 
Furthermore, we can calculate 

\k^-kl\ \2miE-V)\ \2m{E-Voo)\ 



^[k{x),ip{x),x] 



m\V{x)-V^\ 



From the previous chapters |4] and |5} we derive 



T>sedi'^(^J i^[k,ip,x]dx^, (6.1.3) 



when hkoo = \/2m{E — Voo), then k^o = ^ ^ . . Now we can find 

n 



■^[A;(x),^(x),x] 



m\V{x)-V^\ \f 



m|y(x) — Vool 1 \ rh 



(6.1.4) 



Therefore, 



T > sech^ ^ y i!}[k, Lp, x] dx^ , 



2 , ;+°° 1 / m 



+00 



as desired. 

Similarly one can find R in terms of [k{x)^ip{x), x\- Note 
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Using {hkooY = 2m{E — Voo), the bounds become 



sech^l-^ / \V-V^\dx], (6.1.6) 



oo ) 

+ 00 



+ 00 



and 



R < tanhM l.T- ^^^ . f \V-V^\dx 



hV 2{E-V, 



oo I 



tanh^j^^^ J"^ \V-V^\dxy (6.1.7) 

These bounds are exact non-perturbative results, however for high energies it 
may be convenient to use the shghtly less restrictive (but analytically much 
more tractable) bounds obtained by simply taking the first non-trivial term 
in the Taylor series. 

When E > V±oo, using {hkoo)'^ = 2m{E — Voo), the bounds become 



T > 



ml f 

1-^ 



2{E - VooW 



5 



2 



(6.1.8) 



and 



m 



This version of the bounds also holds for all energies, but is not very re- 
strictive for low energy. (Somewhat better low energy bounds are developed 
in the chapter dealing with the Miller-Good transformation. See also the 
chapter [o] on black hole greybody factors.) Note that the bounds of this 
subsection make perfectly good sense for both scattering "over the barrier" 
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or "under the barrier" , there is no requirement on the presence or absence of 
classical turning points. 

The transfer matrices can be analyzed by checking that the evolution 
equations simplify to 

da —im{V — Voo) 



dx 



dh _ +im(y - 14 
dx 



{a exp(+2iA;ooa;) + h}. 



(6.1.10) 
(6.1.11) 



This can be written in matrix form as 



d 


a 


_ -im{V - Ko) 


1 exp{—2ikoox) 




a 


dx 


b 




— exp{+2ikoox) —1 




b 










(6.1 



.12) 

This version of the Shabat-Zakharov system [33] has a formal solution in 
terms of the relatively simple transfer matrix 



E{xf, Xi) = V exp 



—im 



(Vix) - Voo) 



1 e 



2'i Aj QT^ oc 



dx I , 



(6.1.13) 

The formal but exact expression for the Bogoliubov coefficients is now 



a p* 
(3 a* 



E(oo, —00 



h^k 



^ exp ( / (Vix) - 



00 J —00 



—e 



1 e 

+2ik 



dx 
(6.1.14) 



Furthermore, the form of the system (6.1.10 )-(6. 1.11 ) suggests that it might 
be useful to define 



a = a exp 



b = b exp 



+ 



hp knr 



im 

KPknc 



{V{y)-Vno)dy 
{V{y)~Vno)dy 



(6.1.15) 
(6.1.16) 
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Figure 6.1: This diagram (qualitatively) shows the potential V{x), and as- 
sociated function k'^{x), for the scattering problem considered in special 
case 1 1231. 



Then we can substitute equations (6.1.15) and (6.1.16) into equation (6.1.12). 
We derive 



da —imiy^x) 
dx K^kn 



b exp{—2ikoox) 



(6.1.17) 



and 



db +im(y{x) — Vo, 



a exp{+2ika 



.X 



(6.1.18) 
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respectively. This representation simplifies some of the results, for instance 
= £^(oo,— oo), 



a p* 



V exp 



—im 



{v{x) - Ko) 



— e 



e 



Also note that 



a — a exp 



P = P exp 



+ 



im 


r+oo 




/ — oo 




r+oo 


im 




' —oo 



{V{y) - V^) dy 



iViy)-Voo)dy 



dx 
(6.1.19) 



(6.1.20) 



(6.1.21) 



This can be used as the basis of an approximation scheme for p. Suppose 
that for all x we have \b{x) \ <S 1, so that \a{x)\ 1. Then 



db +im{V{x) - 1^ 



This may be immediately integrated to yield 

+im{V{x) - 14o) 



exp(+2iA;ooa;). 



(6.1.22) 



P 



exp(+2i/coox) dx. 



(6.1.23) 



This is immediately recognizable as the (first) Born approximation. If we 
instead work in terms of the original definition P we obtain a slightly different 
approximation 

»+oo 



p 



K^knr, 



exp 



X 



tm 



t9k 



{V{x) - V^) dx 



oo J — oo 



{V{x) - Voo) exp{+2ikoox) exp 



im 



{V{x)-Vo,)dy 



oo J —oo 



dx. 



(6.1.24) 



This is one form of the distorted Born wave approximation. 
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Definition (The Born approximation in scattering theory): 

The Born approximation consists of taking the unperturbed incident 
field, in place of the total field, as the driving field at each point in 
the scatterer. It is the perturbation technique most commonly applied 
to scattering by an extended body. It is a good approximation if the 
scattered field is small, compared to the incident field, everywhere 
inside the scatterer [27] . 

For our purposes, the Born approximation is a good approximation 
whenever the "reflected wave", characterized in our notation by b{x), 
is small everywhere, \b{x)\ ^ 1. (This is a stronger statement than 
merely saying <^ 1.) Similarly, one requires the reflection ampli- 
tude to satisfy |r| ^ 1. 



In brief, the analysis of this chapter has so far collected together a large 
number of results that otherwise appear quite unrelated. By taking further 
speciflc cases of these bounds, and related results, it is possible to reproduce 
many analytically known results, such as those for delta-function potentials, 
double-delta-function potentials, square wells, and sech^ potentials, as dis- 
cussed later in this chapter. 



6.2 Bounds: special case 2 

Having developed a good understanding of the bounds for special case 1, 
now we reconsider the bounds for special case 2 in ^B]: Suppose now we 
take k{x) = (p'{x) 7^ and again set x = 0. This implies that we are 
choosing our auxiliary function so that we use the WKB approximation for 
the true wavefunction as a "basis" for calculating the Bogoliubov coefficients. 

This choice is perfectly capable of dealing with the case V+00 7^ V-00, but 
by reason of the assumed reality of (p is limited to considering scattering over 
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the potential barrier. (This is the special case implicit in a different context 
in [55]). The evolution equations again simplify tremendously. 



When we substitute k(x) = (p'{x) and x = i^^to equation (4.5.55) and 



equation (4.5.56), we obtain: 



da 

dx 

db 

dx 



+ ^|<^"&exp(-2«(^)|, 
+ ^|<^"aexp(+2iy?)|. 



(6.2.1) 
(6.2.2) 



This form of the evolution equations can be related to the qualitative discus- 
sion of scattering "over a potential barrier" presented by Migdal [571158]. For 



this choice of auxiliary functions we consider the equation (6.1.1) to derive 



'&[k{x),(p{x),x] 



\k'\ 

" 2\ip'\ ~ 2\^\ ~ 2\k\'' 
2|(^'| 2\ky 



and the bounds become 



T > sech 

and 

R < tanh^ 
The relevant transfer matrix is now 



+~ \k'\ 

~\h\ 
\k'\ 

\k\ 



dx 



dx . 



E{xf, x-i) = V exp 
The Bogoliubov coefficients are now 
= E[oo, — oo) = V exp 

P a* 














_g+2*^ 



dx 



(6.2.3) 
(6.2.4) 

(6.2.5) 
(6.2.6) 

(6.2.7) 



-2iip 



\-2iip g 



dx 



(6.2.8) 
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6.3 Reflection above the barrier 



The system (6.2.1)-(6.2.2) can also be used as the basis of an approximation 
scheme for /3. Suppose that for all x we have \b{x)\ <C 1, so that |a(x)| ~ 1. 
Then 

db if" , , , , 

Tx^k ^'-'-'^ 

This may be immediately integrated to yield 

^ ^ 9 / ^ exp{+2zv) dx. (6.3.2) 

Or the equivalent 

This result serves to clarify the otherwise quite mysterious discussion of so- 
called "reflection above the barrier" given by Migdal [571158]. Even though 
the WKB wavefunctions are buried in the representation of the wavefunction 
underlying the analysis leading to this approximation, the validity of this 
result for does not require validity of the WKB approximation. 

If the shifted potential, V — V^, is "small" then we can recover the 
Born approximation in the usual manner. In that case k' = mV / {h^k) 
mV / {h'^ koo) , while exp(2z J k) ^ exp{2ikoox). A single integration by parts 
then yields 

/? ~ -^T^ / - Voo) exp(+22A;oox) dx. (6.3.4) 

This is now equivalent to the first Born approximation, in this particular 
context. 



6.4 Under the barrier? 



What goes wrong when we try to extend this analysis into the classically for- 



bidden region? Analytically continuing the system ( 6.2. !)-( 6.2.2) is trivial 
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Quantum effect: The reflection probability R = |rp is a function of 
the energy of the incoming particles. It can always in principle be eval- 
uated by numerically integrating the Schrodinger equation for different 
values of kinetic energy, E and potential energy, V{x). (However, the 
central idea of this thesis is to see just how much we can do ana- 
lytically, without resorting to numerical integration.) Furthermore, 
"quantum effects", such as over-barrier reflection and under-harrier 
transmission, are dominated by regions of the potential where the 
naive semiclassical treatment fails. Due to tunneling, the reflection 
probabilities at energies above the barrier are larger than zero, and 
the under- barrier reflection probabilities are smaller than one f28j. 



replace 



therefore 



_ 2m[E - V{x)] _ 2m[V{x) - E] 

- ^2 - ^ ' ^^-^-^^ 

j2m\V(x) - E] 
k = i, (6.4.2) 



and (as required) 



Now let 



and write 



,p'(x) = k^iK = i^/2m{V - E)/h . (6.4.3) 
K = yj2m{V - E)/h, (6.4.4) 



ip{x) = iptp + i K,{y)dy. (6.4.5) 
Jtp 

We will now use these results to re-analyze the Shabat-Zakharov system. 



X 
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Substitute the two above equations into equation (6.2.1), we obtain 
+ ^|<^"^exp(-2i(^)|, 



da 
dx 



+ £;&exp(-2iv9tp)exp ^ + 2^ K{y) dy 



+ —b exp{-2iiftp) exp 



K 



(6.4.6) 



Similarly to the equation (|6.2.2|), we obtain 
db 



dx 



+ ^|¥'"«exp(+2i^)| 



+ — aexp(+22V9), 
2k 



+ — aexp + 
2k, 



2i( 



ftp + i 



tp 



+ — a exp(+2iv5tp) exp 
2k 



K{y) dy 
K{y) dy 



tp 



K 



+ — aexp{+2i(ptp) exp 



(6.4.7) 



Thus we are now violating our previous condition that (f be real, though 
we still require ip' ^ 0. This is a perfectly good Shabat-Zakharov system 
that works in the forbidden region. But you cannot now use this to derive 
bounds on the transmission coefficient. The fly in the ointment resides in the 
fact that the formula for the probability current is modified, and that in the 
forbidden region the probability current is 

^ = Im |?^*^| = Im{a6* exp(+2i(^tp)}- (6.4.8) 

We can derive the above equation by combining 



Im < t/' 



dx 



(6.4.9) 
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and 

/ = {\a\'-m. (6.4.10) 

For a properly normalized flux in the allowed region (|ap — |6p = 1), we have 
in the forbidden region 

lm{ab*exp{+2iiptp)} = I. (6.4.11) 

While this does imply 2|a| |6| > 1, the inequality is unfortunately in the wrong 
direction to be useful for placing bounds on the transmission coefficient. It 
is for this reason that we have gone to the trouble of keeping track of the 
more general gauge condition represented by x — in the hope that we can 
use this to explore under the barrier. 

A good rigorous bound (not just a WKB approximation) on transmission 
under the barrier would be very useful. The best we have been able to 



do along these lines is presented in chapter [TO] discussing the Miller-Good 
transformation. 



6.5 Special case 2-a 



Suppose now that V{x) is continuous and monotonic increasing or decreasing, 
varying from V^oo = V"(— oo) to V+oo = V{+oo). 

Suppose E > max{V_ooj ^+00} so there is no classical turning point. Then 

»+oo |^,| 



dx 



In 



+00 



and the transmission and reflection probabilities satisfy 

Ah h 



and 



T > 



R < 



+00 



(^+00 k-oo)"^ 
{k+oo ~\~ k-00)'^ 



(6.5.1) 



(6.5.2) 



(6.5.3) 
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Figure 6.2: Sharp corners are guaranteed to maximize reflection, and abrupt 
transitions are guaranteed to maximize particle production 



Calculation: We consider 



1 \k'\ 
2jk\ 



and the bounds become 
T > sech^ 
> sech^ 
Remembering that 



+°° l\k'\ \ ( \ 

dx I = sech^ 



2 \k\ 



In 



In 



+ 00 



sech(x) 



exp(+x) + exp(— a;) 
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and noting 

exp 
we see that 

Therefore 
so we obtain 



—00/ 



max 



sech 



sech 



In 



■+00 



''+00 



^+00 


1 k-00 


^—00 


y ^+00 




■foo^— 00 


(^+00 


+ k-00) 




k 

00 "'—00 



(^+00 ~l" k^ 



T > 



Ak^c^k. 



{k 



as required. We now consider 



R < tanh^ 



< 



< 



tanh^ 

(^+00 



+00 



+00 



In 



w 



dx 



k-r 



(k 



+00 



k-oo)"^ 



(6.5.7) 

(6.5.8) 

(6.5.9) 
(6.5.10) 

(6.5.11) 



as required. These bounds are instantly understandable as the exact analytic 
results for a step-function potential |37l US, HU |88], and the result asserts 
that for arbitrary smooth monotonic potentials the step function provides 
upper and lower bounds on the exact result. If we are interested in physical 
situations such as a time-dependent refractive index (HHl EI] , or particle pro- 
duction due to the expansion of the universe [63], this technique shows that 
sudden changes in refractive index or size of the universe provide a strict 
upper bound on particle production. 



6.6 Special case 2-b 

Suppose now that V{x) has a single unique extremum (either a peak or 
a valley). Provided that E > max{l^_oo, Kxtremum, ^+00} so that there is no 
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classical turning point, then k{x) moves monotonically from k_oo to /^extremum 
and then back to fc+oo- Under these circumstances we consider 



\k'\ 



dx 



In 
In 
In 



k. 



extremum 



k_ 



+ 



In 



k 



extremum 



''+00 



''extremum 



X 



^extremum 



^extremum 
k k 



(6.6.1) 
(6.6.2) 
(6.6.3) 



as required. 



The bound (4.5.18) implies 

lal < cosh 



and which yields 



< sinh 



In 



In 



k 



extremum 



+00 



''extremum 



(6.6.4) 



(6.6.5) 



Numerous generalizations of these formulae are possible. For example, at the 
cost of a little extra notation, we also already have enough information to 
provide a bound on an asymmetric barrier or asymmetric well, as long as it 
has only a single extremum (maximum or minimum) we apply the previous 
equations to derive 



and 

Calculation: 

lal < cosh 



|a| < 



k"^ 

extremum 



+ k. 



h 

■CO '^—00 



2k 



extremum 



h 

+00 '^—00 



IP 

I extremum 



k-i-nok 



+00"' — 00 I 



2k 



extremum 



\/~k+ook 



(6.6.6) 
(6.6.7) 



00"'— 00 



In 



''extremum 



+00 



exp + 



In 



"-extremum 



exp 



In 



"-extremum 



^extremum 



2k 



extremum 



(6.6.8) 
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as required. 

In addition, we can derive 



1/3 1 < sinh 



In 



extremum 



exp + 



In 



+00 

^extremum 



exp 



In 



"'extremum 



— h h \ 

I "'extremum - 00 "-+00 | 



2k 



extremum 



^/k-ook 



(6.6.9) 



+CXD 



as required. 

Translated into statements about the transmission and reflection proba- 
bilities this becomes 

— 00 



and 

Calculation: 

T > sech^ 
= sech^ 



T > 



R < 



''+00"'— 00 ""extremum 

{^extremum "I" ^+00^-00}^ 



{^extremum ^+00^-00} 
{^extremum ~l~ ^+00^-00}^ 



(6.6.10) 



(6.6.11) 



2k 



dx I = sech 



In 



"'extremum 



, (6.6.12) 



In 



\/k-ook-\ 



+00 

extremum 



extremum 



+ ^+00^-00}^ 



Also 



R < tanh^ 



2k 



dx I — tanh 



In 



^extremum 



y/k-ook 



+00 



But 



tanh(x) 



(cxp(2,T) - 1) 
(exp(2a;) + 1) 
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therefore 
implying 

as required. 
Equivalently 

T > 



R < tanh^ 



In 



^extremum 



+00 



R < 



{^extremum ^+00 ^-00} 
{^extremum "I" ^+00^-00}^ 



- Kxtremum) " ^+00) " ^-00) 



and 



R < 



[{E - Kxtremum) + ^/ {E - V+^) {E - V^^)]^ 
[{E — Kxtrcmum) " — V+oo){E — V-oc 



[{E - V, 



extremum I 



y/{E-V+^)iE-V.^W' 



(6.6.15) 

(6.6.16) 

□ 

(6.6.17) 
(6.6.18) 



Calculation: We consider 



''it 00 



2m {E - V±c 



h 



and 



k 



extremum 



IfTl (yE V^xtremum) 



h 



(6.6.19) 



(6.6.20) 



we substitute the above equations into equation (6.6.10) to derive equation 



(6.6.17). Moreover, we substitute the above equation into equation (6.6.11) 



to derive equation (6.6.18). 



□ 



This can be compared, for example, with known analytic results for the 
asymmetric square well. To be more specific, if in addition V^(— oo) = 



V{+oo) 



Voo, so that k. 



— h 

OO "^+00 



fcoo, then we have 



and 



\a\ < 



extremum ' oo 

9h . h ' 

"'extremum "'OO 



\h'2 _ t.2 I 

I extremum oo I 

Oh h 

^ "'extremum "'oo 
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Translated into statements about the transmission and reflection probabili- 
ties this becomes 

^ - (F-v \(F-v ^ , in/ riT^' [^■^■^^) 

\^ ''oo)\^ cxtrcmum J ~r ^ v cxtrcmum ''ooj 



and 



{^E V^)(-E' ^xtremum) ~l~ 4(^xtremum ^'^)^ 

Equivalently 

(^xtrcmum 
i^E VJ^xtrcmum ^ 

and 2 

(^xtremum ^'^) ^ f-><-\ 

^ - r9F-T/ -V ^2 - (6.6.26) 

^^ZJZ/ f extremum ^00) 

For low energies, these results are weaker than the bounds derived under 



T>1- (6.6.25) 



special case 1, (6.1.6, 6.1.7) and (6.1.8, 6.1.9), but have the advantage of 
requiring more selective information about the potential. For high energies, 

E > ^'(^extremum- V^o)' g 27) 

2m{j:^-\V{x)-V^\dxY 

the present result (when it is applicable) leads to tighter bounds on the 
transmission and reflection coefficients. 



6.7 Special case 2-c 



Suppose now that V{x) has a number of extrema, (both peaks and val- 
leys). We allow V{+oo) ^ V{—oo), but demand that for all extrema E > 
max{V^_oo, V+oo, Kxtremum} SO that there is no classical turning point. 

For definiteness, suppose the ordering is: — oo peak — > valley . . . valley 

peak ^ +00. Then 



+00 



\k'\ 
~k 



dx 



In 



^1 

"'peak 



In 



"'valley 
"'peak 



+ 



In 



"'peak 
'^valley 



+ 



In 



''peak 



(6.7.1) 
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Defining 





— ^peak' 




peaks 




— 11 "'valley) 




valley 




J_ J_ extremum 



Up{k) U,{k), 



extrema 



we see 



+00 



k'\ 



dx 



In 



k-ook+oo^vik) 



k.^k+ooIll{k) + Ul{k) 



This bounds the Bogoliubov coefficients as 
and 



k^^k^^nlik)-ulik)\ 



Then the transmission and reflection probabihties satisfy 

4k+^k^^Ul{k) 



T > 



and 



R < 



{Ulik) + k^^k_^Ul{k)} 
^Ilp{k) — k^ook-oo^vik)} 



2 ' 



2 • 



(6.7.2) 
(6.7.3) 
(6.7.4) 

(6.7.5) 



(6.7.6) 



(6.7.7) 



(6.7.8) 



(6.7.9) 



{Ulik) + k^^k_^Ul{k)y 
In these formulae, peaks and valleys can be interchanged in the obvious way, 
and by letting the initial or flnal peak sink down to Kj-00 as appropriate we 
obtain bounds for sequences such as: — 00 ^ valley peak . . . valley — » peak 
— » +00, or: —00 peak valley . . .peak valley +00. In the case 
of one or zero extrema these formulae reduce to the previously given results. 



[Equations (6.6.10)-(6.6.11 ).] Further modiflcations of these formulae are 
still possible, the cost is that more speciflc assumptions are needed to derive 
more speciflc results. 
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6.8 Bounds: Special case 3 



In the following we will consider the bounds in special case 3. In particular, 
the most outstanding features of this case is: 
Let X = 0, /co > and pick 



cp' — ^ max{k'^{x), /cq}, 



(6.8.1) 



with Xq defined by k'^ix^) — /cq. Then we have 



1 \k'\ 
2M 



1 kl - e 



I 2 ko 



e < ki 



(6.8.2) 



Note that there are step function discontinuities at Xq, but no delta- function 
contribution. It now follows that 



= - In 
2 



k. 



ko 



+ ^/ [kl-k']dx + lln 
^ko Jk2<B 2 



k 



+ 00 

ko 



(6.8.3) 



that is 



= - In 



k—cok 



0O"'+0O 



+ 



2k 



Jfe2<fe2 



[kl - dx. (6.8.4) 



Consider 



so 



/ {-k^)dx= I K^dx, 



(6.8.5) 



jr [ [kl-k'^]dx < I K^dx + klLk2<kl], 



(6.8.6) 



So collecting terms we have 

n^— oo'^+oo 



In 



A-2 



1 /■ 1 

+ ^ / «^ + - /co i-fe2<fe2 ■ (6.8.7) 
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Now note 

T > sech^ 



> sech I - In 



k—nck. 



OO "^+00 



+ 



^ y + ^ /co ivfc2<fc2^ , (6.8.8) 



so that 
T > 

where 



exp(S(x)) + ^ ^° exp(-S(x)) 



+0O 



2 ' 



(6.8.9) 



If 1 
(2^) = ^ / «^ + - /co i-fe2<fe2, (6.8.10) 

when /co is still an adjustable parametric, though we definitely need < /cq < 
min{/c+oo, ^-oo}- 



6.9 Bounds: Special case 4 

For the fourth special case we want to derive something that looks simi- 
lar to the WKB approximation, but is a strict bound instead of being an 
(uncontrolled) estimate. 

Choose X = 0, /co > and pick 

if' = ^max{|A;2(x)|,A;2}. (6.9.1) 

Even for a potential with only a single hump there are now five regions to 
analyze, the two allowed regions, the forbidden region, and two transition 
regions enclosing the two the classical turning points. As usual we shall 
define our notation so that in the forbidden region: 

k'^{x) = -k^{x)] K > 0. (6.9.2) 

In the two allowed regions 

k':,kl (6.9.3) 
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In the two transition regions 

2 ko '■ 
Finally in the forbidden region 



So in the forbidden region by the triangle inequahty 



(6.9.4) 



(6.9.5) 



(6.9.6) 



Collecting these 



< 



1? < -In 



+ 



2A;, 



1.2 

decreasing 



Jkl>k-^>-kl 



'^extremum 

ko 

''inclreasing 



If 1 
- / K{x)dx+-ln 



'^extremum 

ko 



1 r inclreasing T 

f— / [^o^-^^]dx+-ln 
^/^o J-kl<k-^<kl ^ 



Now in each of the transition regions 



■'+00 



SO 



[fcg-fc^jda; < 2/c^ / Idx, 

i-itg<ik2<it2 



— 2/Co L_fc2<;,2<fc2. 



(6.9.7) 

(6.9.8) 

(6.9.9) 
(6.9.10) 



Here L_^2<fc2<fc2 is the combined length of the transition regions — /cq < k'^ < 
/cq. That is, the integral coming from each transition region is bounded (up 
to a constant) by the physical width of that transition region. Furthermore 
in the forbidden region 



/ k{x) dx < k{x) dx, 



(6.9.11) 
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so collecting terms we have 



?9 < - In 
- 2 



oo ^+ oo ^extremum 



+ ^0 -^-fc2<fc2<fc2 + 



1 



«;2>0 



n{x)Ax, (6.9.12) 



where L now denotes the combined total width of the two transition regions. 
Now note 



T > sech^ (b '§> exp 



so that 
T > 



exp[— 2A;o I/(A;o)] exp 



(6.9.13) 



(6.9.14) 



^-oo^+oo'^extremum L J k'^>G 

The bound is considerably weaker than could have been derived by making 
more restrictive hypotheses, but has the advantage of elegance and looking 
very similar to the WKB bound. Note the key requirements: We must be 
dealing with a single-hump potential, and must be in the range 



< /Cq < min{/c_oo, /c_|_oo, ft^extremum}- 



(6.9.15) 



The parameter k^ is otherwise arbitrary, and so can be chosen to maximize 
the pref actor. 

One thing we could do is to choose a different value of ko at each tran- 
sition, call them ki and ^2, and repeat the analysis keeping careful track of 
the K integration near the turning points. Then 



j, < 



In 



^-oo ^+oo /^cxtrcmum 



k\ ^2 



+ 
+ 
+ - 



2^1 J-kl<k'^<kl 
1 



[kl - A;2] dx 



2^2 J-kl<k'^<kl 
1 



[kl - A;2] da; 



2 Jkl>K^>o 
1 



nix) dx 



n{x) dx 



«;2>0 



n{x) dx. 



So far this is a rigorous bound; now we are going to adopt a linear approxi- 
mation near the turning points, so that near the first turning point 

A;^(x) si(x-xi), (6.9.16) 
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with k{x) reaching the values ±kl at the positions x — Xi± kf/si. Therefore 



Li = 2kf/si. Then 

/ 

Similarly 



[/ci - A;^] dx fa / [kf - si{x - xi)]dx = — -. 

-kl<k'^<kl Jxi-kf/si *i 



/ K{x)dx ~ / ^/si{x — xi) dx ^ 

Jkl>K^>0 Jxi 



(6.9.17) 



(6.9.18) 



'XI 3 Si 

so that combining, and assuming the linear approximation is a good one 



< - In 
- 2 



fe-oo ^+oo ^extremum 
21,2 



+ 



+ 



+ 



2) S\ 3 §2 2 J 



ti{x)dx. (6.9.19) 



Extremize with respect to /ci, then 

ki si 
and similarly for A;2, then 



so kl 



(6.9.20) 



1 



1? < - In 
- 2 



oo '''+(X} '^extrcmum 
((l/2)5i)2/3 ((1/2)^2)2/3, 



2 1 

+ 7: + 7: I K{x)dx, (6.9.21) 



and hence 

((l/2)s0^/^ ((1/2)52)2/^ 



T > 



exp[— 4/3] exp 



k2>0 



«;(x) dx 



(6.9.22) 



- oo "'+00 "'extremum 

Although such results are certainly a major advance in our standing of the 
rigorous bounds, this particular bound is not 100% rigorous due to the linear 
approximation. This suggests that further exploration of these ideas, with a 
view to obtaining a 100% rigorous bound, might be profitable. 



6.10 Bounds: Future directions 

Prom the general definition 

= - ' (6-10-1) 



■d[k{x),(p{x),xix)] 



2\ip' 
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and the bound 

T > sech^ 

there are many other special cases you could in principle derive. The possi- 
bilities seem endless and the art is in finding something useful. 



dx 



(6.10.2) 



6.11 Discussion 



In this chapter we dealt with some specific cases of the general bounds and 
developed a number of interesting specializations. In addition, we collected 
together a large number of results that otherwise appeared quite unrelated, 
including refiection above and below the barrier. The special case bounds 
were divided into five topics: special cases 1-4, and "future directions". In 
addition, all special cases were chosen for their directness and simplicity. 



Furthermore, special case 1, as presented in equations (6.1.6 )-(6. 1.7) and 



(6.1.8)-(6.1.9), has the advantage that it applies to both scattering over the 



barrier and under the barrier. On the other hand, special presented 



in equations (6.2.5)-(6.2.6) and their specializations, applies only to scatter- 



ing over the barrier but has the advantage of being much more selective in 
how much information is needed concerning the scattering potential. 
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Chapter 7 



Parametric oscillations 



7.1 Introduction 



In this chapter we shall present the basic concept of a "parametric oscilla- 
tor". This is a simple harmonic oscillator whose parameters (its resonance 
frequency uj and damping time (3) vary in time. The other interesting way of 
understanding a parametric oscillator is that it is a device that oscillates when 
one of its "parameters" (a physical entity, like capacitance) is changed |119j . 

We shall re-cast and represent the general bounds in terms of the spe- 
cific mathematical structure of parametric oscillations. This time-dependent 
problem is closely related to the spatial properties of the time-independent 
Schrodinger equation. 

Although the discussion so far has been presented in terms of the spatial 
properties of the time-independent Schrodinger equation, the mathematical 
structure of parametrically excited oscillations is identical, needing only a few 
minor translations to be brought into the current form. For a parametrically 
excited oscillator we have 



Just map t — > x, uj{t) k{x), and (p ^ ip. In the general analysis of 




di2 
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equation (6.1.1) the quantity should be replaced by 



%(t),^(t)]E 

This is often written as 



2|(^'| 



(7.1.2) 



2\ip\ 



(7.1.3) 



but this is purely a convention, a change in notation. (Physicists typically 
use dots for time derivatives and primes for space derivatives.) The analysis 
then parallels that of the Schrodinger equation. Some key results are given 
below. 



7.2 Special case 1 

If ti;(— oo) = ujq = ci;(+oo) 7^ 0, then by choosing the auxiliary function to be 



(f = Uot we can use equations (5.3.7)-(5.3.14) to deduce 



lal < cosh 



1 



+00 



and 



2uJo 



(7.2.1) 



(7.2.2) 



Calculation: We now consider 



|a| < cosh|^ j i9[ip{t),uj{t)] dt 



(7.2.3) 



We substitute equation (7.1.2) into equation (7.2.3), now we derive 



lal < cosh 



2\^'\ 



dt 



(7.2.4) 
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A 








bp - m'' = 1 





Figure 7.1: For a parametric oscillator incoming "ground state" fluctuations 
from the past, (left, coefficient 1), can in the future be amplified, (right, 
coefficient a). Quantum mechanically, the coefficient (3 is related to particle 
excitation from the ground state [23] ■ The quantities a and P are typically 
referred to as the Bogoliubov coefficients. 
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Furthermore, under the stated assumptions we can simphfy the above integral 
to derive 

lal < cosh( / \uj^{t) -ulldt), (7.2.5) 



2uJo 



which automatically implies 
1/5 1 < sinh 

as required. 



(7.2.6) 



7.3 Special case 2 

If uj{—oo) and ti;(+oo) ^ are both finite so that suitable asymptotic states 
exist, and assuming co'^(t) > so that the frequency is always positive, then 



applying equations (6.2.5)-(6.2.6) to the case of parametric resonance yields 

-iK(t)| 



and 



la I < cosh 



\(3\ < sinh 



2 \uj{t)\ 

l\um\ 
2 \uj{t)\ 



dt 



dt 



(7.3.1) 



(7.3.2) 



7.4 Special case 2-a 

Suppose now that uj^it) is positive semidefinite, continuous, and monotonic 
increasing or decreasing, varying from uj^oo = uj{—oo) 7^ to some distinct 
value uj+00 = uj{+C!c) 7^ 0. The Bogoliubov coefficients satisfy 



and 



lal < 



< 



uj-00 + ^ 
2^a;_ooto'_|_c 



\L0- 



(7.4.1) 



(7.4.2) 
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Calculation: We now consider 



\a\ 



< 



cosh 


1 lu'it) 


J-oo 2 \uj{t)\ 


cosh 






2 \0U-ooJ 


<^-oo 




2^u; 


) 



dt 



(7.4.3) 



as required. Similarly 



\(3\ < sinh 
< sinh 



2 |cj(t)| 



dt 



< 



\UJ-r 



ilnf^ 



2^a;_oo<^+o 



(7.4.4) 



as required. 



7.5 Special case 2-b 

Under the restriction a;(— oo) — cuq — a;(+oo) ^ 0, with the additional 
constraint that u!{t) has a single unique extremum (either a maximum or a 
minimum but not both), and provided that i^gxtremum > so that we do not 
encounter complex frequencies (no classical turning point), the Bogoliubov 
coefficients satisfy 

|a| < ^0 + ^extremum ^ ^7 5 -^^ 

2ti^0^extremum 



and 



ltd a - cu,2 



\P\ < ^ -extremum I _ (7.5.2) 

2^<^0^extremum 
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Calculation: To prove the above equations, we consider 

f+OO 



/" 

J — ( 



\0J I 



dx 





In 


'^extremum 
















- .2 






In 


extremum 
















" ,2 






In 


^extremum 








. ^0 . 





+ 



In 



^extremum 



(7.5.3) 
(7.5.4) 



This imphes 



\a\ 



< cosh 

< cosh 



In 



2 
In 



extremum 







(7.5.5) 



'-'-'extremum 



< 



2 2 
^0 "I" "^extremum 

2u!oU!, 



extremum 



Furthermore, 



\P\ < sinh 



< sinh 



In 



2 
In 



extremum 







(7.5.6) 



< 



^extremum 



'^extremum 



extremum 



Suppose now that a;^(t) has a single unique extremum (either a peak or a 
valley), but that we allow the two asymptotic frequencies to differ a;(+cxo) ^ 
C(;(— oo), and suppose further that uj^ii) > so that there is no classical 
turning point. The Bogoliubov coefficients satisfy 



and 



lal < 



< 



'^-oo^-foo ~l~ ^extremum 



2 -y/ ^—oo^+oo ''-^extremum 



I'^-oo'^+oo 



extremum I 



2 \J ^—oo^-foo ''-^extremum 



(7.5.7) 



(7.5.8) 
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Calculation: We now show how to derive the above equations: For the 
moment we shall consider 



\a\ < cosh 



< 



In 



'^extremum 

^extremum ~l~ ^-oo^+oo 
2 00^*^+00 ''^extremum 



But this now yields 



(7.5.9) 



< 



extremum 



2 ^'^—oo '^+00 '^extremum 



(7.5.10) 



as required. 



7.6 Special case 2-c 

Suppose now that uj{t) has a number of extrema (both peaks and valleys). 
We allow a;(+oo) 7^ a;(— 00), but demand that for all extrema (^extremum > 
so that there is no classical turning point. 

For definiteness, suppose the ordering is: —00 — > peak — > valley . . . valley 

peak — > +00. Define 

(7.6.1) 
(7.6.2) 
(7.6.3) 

extrema 

The Bogohubov coefficients satisfy 

a < -rj r , (7.6.4) 

^CU+oo^-oonel'^) 

and 

xAW^^ne(a;) 
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peaks 


n^(c^) 


= H '^valley 




valleys 




— II ''^extremum 
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In these formulae, peaks and valleys can be interchanged in the obvious way, 
and by letting the initial or final peak sink down to a;±oo as appropriate we 
obtain bounds for sequences such as: — oo — > valley —>■ peak . . . valley 
peak +00, or: —00 — > peak valley . . . peak — > valley +00. In the 
case of one or zero extrema these formulae reduce to the previously given 
results. 

Again, further specializations of these formulae are still possible. As 
always there is a trade-off between the strength of the result and its generality. 

Calculation: We see 




In 




(7.6.6) 



The Bogoliubov coefficients in this case are bounded by 



a < 



o;_ooW+oon^.(cj) + Ul 



(7.6.7) 



2^cc;+oo^_oone(a;) 



and 



m< 



HI 



(7.6.8) 



7.7 Bounds: Special case 3 



We let X = 0, cuo > 0, then we can choose 




(7.7.1) 



with defined by uj'^{t^) — uj'{ 



n. Then we have 




(7.7.2) 
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Note that there are step function discontinuities at t^, but no delta-function 
contribution. It now follows that 



1?= - In 
2 



That is 



(b -d = - In 




/ 2 


^0 



As usual we shall define our notation so that in the forbidden region: 



(7.7.3) 
(7.7.4) 

(7.7.5) 



So 



(7.7.6) 



Collecting terms we have 



< -In 
- 2 



+ 2^ / + l"^' ^^-^-^^ 



and the bound on the Bogohubov coefficients in this case become 
\a\ < cosh I (p Ti}], 



< - — : exp 



2u;o 



2^a;_oo<^+o 



2uJo 
exp I - 



n dt + -uJoL^2^^2 



2u;o 



1 

-c 
2 



n dt- -a;oL^2<^2 ^ 



(7.7.8) 



and 



\p\ < sinh 



(/4 



< 



c 



2u;o 



exp 



2a;n 



1 

2^ 



dt + -uo L^2^^2 



(7.7.9) 
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which gives us our third special case bound. Then the transmission and 
reflection probabihties satisfy 

T > -^—^^ ^ ^, (7.7.10) 

^/^^^^^^ exp(Z(t)) + exp(-Z(t)) 



^0 ^a;_oo<^+o 



and 



where 



{a;_ooa;+ooexp(Z(t)) -a;gcxp(-Z(t))f ('7 7 11) 

-{a; exp{Z{t)) + uJq exp(-Z(i))}2 ' 

m ^^Jn'dt + ^cuo L.2<,2. (7.7.12) 

7.8 Discussion 

Though the discussion in previous chapters has been presented in terms of 
the spatial properties of the time-independent Schrodinger equation, we have 
seen in this chapter that the mathematical structure of parametrically excited 
oscillations is essentially identical, needing only a few minor translations to 
be brought into the current form. 

In summary, the bounds presented in this chapter are useful in establish- 
ing qualitative analytic properties of parametric oscillators, and as such are 
complementary to both explicit numerical investigations and the guidance 
extracted from exact analytic solutions. 
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Bounding the Bogohubov 
coefficients 

8.1 Introduction 

In this chapter we will again be considering (from a somewhat different point 
of view) the problem of finding approximate solutions for wave equations in 
general, and quantum mechanical problems in particular. It appears that 
as yet relatively little work seems to have been put into the complementary 
problem of founding rigorous bounds on the exact solutions. We have in 
mind either bounds on parametric amplification and the related quantum 
phenomenon of particle production (as encoded in the Bogoliubov coeffi- 
cients), or bounds on transmission and refiection coefficients. 

In the last section of appendix B, we introduce and discuss the time or- 
dering and give some more details of time-ordered exponentials — these are 
a very convenient trick for formally solving certain matrix differential equa- 
tions. Practising physicists and applied mathematicians will all have seen 
the WKB approximation for barrier penetration probability. Unfortunately, 
the WKB approximation is an example of an uncontrolled approximation, 
and in general we do not know if it is an over-estimate or an under-estimate. 

As part of the main work, we modify and improve an approach first 
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developed in [88]. We shall examine this question by developing a formal 
but exact solution for the appropriate second-order linear ODE, in terms of 
a time-ordered exponential of 2 x 2 matrices, then relating the Bogoliubov 
coefficients to certain invariants of this matrix. By bounding the matrix in 
an appropriate manner, we can thereby bound the Bogoliubov coefficients. 



8.2 The second-order ODE 

We would first like to present "the second-order ODE" techniques developed 
in [88], that are applicable to numerous physical situations; situations which 
are both extremely interesting and important. Consider the ODE 

4>{t) + uj^t) (j){t) = 0, (8.2.1) 

or its equivalent in the space domain [HE] 

0"(x) + A;2(x)0(x) = 0. (8.2.2) 



It is easy to see that equation (8.2.1) can be viewed (in terms of the time 



domain) as an example of parametrically excited oscillation; it arises for 
instance when a wave propagates through a medium whose refractive index 
is externally controlled to be a function of time (though remaining spatially 
invariant)]^ In contrast, the spatial version of this equation as presented 



in (8.2.2) arises classically in situations where the refractive index is spatially 
dependent (so called "index gradient" situations), or in a quantum physics 
context when considering the Schrodinger equation for a time-independent 
potential: 

-—(l)"{x) + Vix)(P{x) = E(l){x), (8.2.3) 
2m 



^For instance, situations of this type have been used to model sonoluminescence |112j . 
and more recently both quasiparticle production in analogue spacetimes |113j and analogue 
signature change events |114j . In all these situations it is extremely useful to have rigorous 
and largely model-independent bounds on the amount of particle production that might 
reasonably be expected. 
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as long as one makes the translation 

2m[E - V{x)] 



(8.2.4) 



However they arise, equations (8.2.1) and (8.2.2) are central to the study of 



both quantum physics and wave phenomena generally. As the result of this 



Refractive index (or index of refraction): This index is a measure 
for how much the speed of light (or other waves such as sound waves) 
is reduced inside a medium. The refractive index, n, of a medium is 
defined as the ratio of the phase velocity, c, of a wave phenomenon (for 
instance light or sound) in some reference medium (such as vacuum, 
or air at standard temperature and pressure) to the phase velocity, Vp, 
in the medium itself |105j : 



n 



c 



1.2.5) 



central significance, over the last century or more a vast body of work has 



gone into the question of finding approximate solutions to equations (8.2.1) 



and (8.2.2 ). Most of these approximations are typically based on JWKB tech- 
niques and their variants (phase integral techniques, etc.) |115j . In contrast 
very little work seems to have gone into the physically important question of 
finding explicit bounds on the relevant Bogoliubov coefficients and/or reflec- 
tion and transmission coefficients 188). 



Index-gradient methods: So-called index gradient optics is the 
branch of optics covering optical effects produced by a gradual spatial 
variation of the refractive index of a material |102] . 
One can analogously speak of index gradient acoustics when the speed 
of sound is slowly varying as a function of position. 
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In this chapter we shall modify and streamline the analysis of [88]; pre- 
senting an alternative proof that is considerably more direct and focussed 
than that in [SS]- We can make this discussion appear to be so simple and 
straightforward by assuming that uit) luq (equivalently k{x) ko) out- 
side some region of compact support (equivalently That is, 



concentrating on the time-domain formulation of equation (8.2.1), the quan- 
tity co'^(t) — c<Jq is a function of compact support]^ Because of this compact 
support property we know that everywhere outside the region the ex- 



act solution of the wave equation (8.2.1) is given by linear combinations of 



exp(±zu;ot), and that the central question to be investigated is the manner 
in which exact solutions on the initial domain (— oo, ti) "connect" with exact 
solutions on the final domain {tf^+oo). Our approach will be to focus on 

Comment: Describing and characterizing this "connection" between 
the (known) plane wave solutions to the left and right of the barrier 
is exactly what the Bogoliubov coefficients are designed to do. 
In order to explain these issues in more detail, in the next section we 
shall introduce the time-ordered exponentials, otherwise the rigorous 
bounds we shall derive on the transmission, reflection, and Bogoliubov 
coefficients will be almost impossible to understand. 



using the above deflnition as a guide to the appropriate starting point. We 
can now systematically develop a formal but exact solution for the appropri- 
ate second-order linear ODE in terms of a time-ordered exponential of 2 x 2 
matrices, then relating the Bogoliubov coefficients to certain invariants of 
this matrix. 

^Of course, this "compact support" condition is not strictly necessary, and at the cost of 
a little more analysis one can straightforwardly extend the comments below to a situation 
where there is a finite limit u){t) Wqo as f — > ±oo |88j . At the cost of somewhat more 
tedious additional work, there are also useful things that can be said of the situation where 
uj{t) uJ±oo, with ui-oo (jJ+oo, as t ^ ±00 j55] . 
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We are interested in solving, exactly but possibly formally, the second- 
order ODE 

0(t) +cj2(t)0(t) = 0. (8.2.6) 
One way of proceeding is as follows: Define a momentum 



TT 



(8.2.7) 



and then rewrite the second-order ODE as a system of first-order ODEs 

(8.2.8) 



vr; 

,2/ 



TT = -iu\t) 0; (8.2.9) 

or in matrix notation (where we have carefully arranged all matrix elements 
and vector components to carry the same engineering dimensions) 



d 







dt 







Uq 

-uj^/ujq 



tt/ujo 



.2.10) 



This matrix ODE always has a formal solution in terms of the so-called 
"time-ordered exponential" 



uo 
-uj^{t)/uJo 








")1 


71 /Uo 



to 

(8.2.11) 

The meaning of the time-ordered exponential is somewhat tricky, but ulti- 
mately is just a 2 X 2 matrix specialization of the operator- valued version 
of the "time-ordered exponential" familiar from developing quantum field 
theoretic perturbation theory in the so-called "interaction picture" |116j . 
Specifically, let us partition the interval {tQ,t) as follows: 



to < ti < t2 < t^... < tn^s < tn-2 < tji-l < tn 

and define the "mesh" as 



(8.2.12) 



M = max {ti — ti_i}. 

ie(l,n) 



(8.2.13) 
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Then define the time-ordered exponential as 



T{t) 







At 



hm 

M^O, (n^oo) 



n-1 

JJexp 

i=0 



UJQ 

-u/{tn-i)/ujQ 



{tn—i tn—i—l] 



.2.14) 



Comment: In this matrix product "late times" are always ordered 
to the left, and "early times" to the right. Now we can extract all 
the interesting physics by working with this time-ordered matrix. (If 
we work in the space domain then the equivalent matrix T is "path- 
ordered", and is closely related to the so-called "transfer matrix".) 

• Since all of the "complicated" physics takes place for t e {ti,tf), 
it is also useful to define 



T ^ £r { cxp 












1= 


a b 









c d 



(8.2.15) 



We are guaranteed that det[T] = 1, that is ad — be = 1. This fol- 
lows from the fact that det[T] — exp{tr(ln[T])}, and the explicit 
formula for T above. 



• Another particularly nice feature is that with the current defini- 
tions the transfer matrix T is manifestly real. This is relatively 
rare when setting up scattering or particle production problems, 
so we shall make the most of it. 
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8.3 Bogoliubov coefficients 



We have already been introduced to the concept of Bogohubov coefficients, 
and also some ways to calculate them, in chapter [5j Let us now extract the 
Bogoliubov coefficients in the present situation. Before tj, and after tj, the 
wave-function is just linear combinations of exp{±iu!ot). We can prepare 
things so that before ti the wavefunction is pure exp{+iuQt), 

iP{t <U) = exp{+iiUot); (8.3.1) 

in which case after tf the wavefunction will be a linear combination 

"^{t > tf) = aexp{+iuQt) + /5exp(— icuo t), (8.3.2) 

where the Bogoliubov coefficients a and (3 are to be calculated. That is, we 
have 








exp{+iuJo ti) 




U 


i exp{+iuo ti) 



(8.3.3) 



and 










if 



aexp{+iuotf) + Pexp^—iuotj) 
i {a exp{+iuQtf) — [3 exp(— icjo t/)} 



But we also have 



implying 






= T 





Tx/uJo 







(8.3.4) 



(8.3.5) 



aexp{+iuQtf) + /5exp(— zcuq t/) 
i {a exp{+iuQtf) — f3 exp(— icjo t/)} 



a exp(+ico'o ti) + bi exp{+iuJo ti) 
cexp{+iu!oti) + diexp{+iu!Qti) 

(8.3.6) 



Solving these simultaneous linear equations we find 

1 
2 
1 
2 



a 



13 



[a + d + i{b- c)] exp(-ztuo [tf - ti]), 
[a - d + i{b + c)] exp{-iuJo [t/ + ti]), 
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so that the Bogohubov coefficients are simple hnear combinations of elements 
of the matrix T. Then (remember the matrix T is real) 

1 



\a 



^ = -{{a + df + {b-cf}, (8.3.9) 

1 
4 



'2- ^ {(a-d)' + {b + c)^}, (8.3.10) 



and so 



= 5^1^^^£±M^ = ad-6c=l, (8.3.12) 

thus verifying that, (thanks to the unit determinant condition), the Bogoh- 
ubov coefficients are properly normalized. Particle production is governed 
by the /3 coefficient in the combination 

- \{{a-df + {b + cf}, (8.3.13) 

= ^{a"^ + d'^ -2ad + b^ + + 2bc} , (8.3.14) 

= ^{a^ + d^ + b^ + c''-2}, (8.3.15) 

= ^tr{TT^-/}. (8.3.16) 

Note that the transpose is now time- anti- ordered. 
Similarly, we have 

lap 



\{{a+dr+{b-cr}, 


(8.3.17) 


^ {a^ + + 2ad + 6^ + - 26c} , 


(8.3.18) 


^{a' + d' + b' + c' + 2}, 


(8.3.19) 


^-tr{TT^ + I}. 


(8.3.20) 


it is now obvious that 




i«r - i/^r = ^ tr{7} = 1. 


(8.3.21) 
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Comment: We can always formally solve the relevant ODE, either 
equation (8.2.1) or its equivalent equation (8.2.2), in terms of the 
time-ordered exponential, and we can always formally extract the Bo- 
goliubov coefficients in terms of traces of the form tr{T T^}. We shall 
now use these formal results to derive rigorous bounds on the Bogoli- 
ubov coefficients. 



8.4 Elementary bound: 



Now consider the quantity 



X{t) = T{t)T{tf 
= ^ < exp 



^ < exp 



UJQ 

-iu^{i)/iUo 



dt 



Uq 



-I T 



dt 



(8.4.1) 



(8.4.2) 



This object satisfies the differential equation 



dX 



uo 
-uj^{t)/uJo 



X{t)+X{t) 



-uj^{t)/uJo 
uo 



(8.4.3) 



with the boundary condition 



XiU) = I. 



(8.4.4) 



Now note 



tr(X) = tr{TT^} = + b"^ + + 



14.5) 
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Furthermore 



dX 



+ 



uo 
—up-jijjQ 

a} + 6^ ac + hd 
ac + bd + 



+ ac + hd 
ac + bd + 

—u^/cuo 
ujo 



2u;o{ac + hd) 

uJo{^ + d^) - (a;Va;o)(a' + {-2uj^ / uJo){ac + M) 



^.4.6) 



and so we see 
tr 



LOq 

up-jojQ 



Therefore 



dtr[X] 
dt 



—uj'^/ujo 

UJQ 



2(ac + 



a; 



2{ac + bd) 



Uo - 



(8.4.7) 
(8.4.8) 



Using this key result, and some very simple analysis, we shall now derive our 
first elementary bound on the Bogoliubov coefficients. 



For any 2 real numbers, using (a; + y)^ > and (a; — y)^ > 0, we have 

+ y'^ >2\xy\. (8.4.9) 
In particular, for any 4 real numbers this implies 

a' + 6^ + c2 + ^2 > 2y^(a2 + lP){c^ + d"). (8.4.10) 

But we also have 

\ac + hd\'^ + \ad - hc\^ = a^c^ + 2abcd + b'^d'^ + a^d^ - 2abcd + b'^c^, 

(8.4.11) 

= {a^ + b^){c^ + d^), (8.4.12) 
thus, for any 4 real numbers 

a'^ + b'^ + c^ + d^>2 y/\ac + bd\'^ + \ad-bc\'^. (8.4.13) 
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• For the particular case we are interested in we additionally have the 
unit determinant condition ad — 6c = 1, so the above implies 



whence 



a^ + b^ + c^ + d^> 2^y\ac + bd\^ + 1, 



2\ac + bd\ < v/(a^ + b^ + c'^ + d^f - 4. 



Then, collecting these results, we see 



dtr[X] 



dt 



= 2{ac + bd) 



UJQ 



< 2\ac + bd\ 



whence 



dtr[X] 
dt 



< + 62 + c2 + d^f - 4 



<^0 



v/tr[X]2 _ 4 



Wo - 



whence 



This imphes 



whence 



We now have 



so that 



1 dtr[X] 



^tr[X]2 - 4 dt 

dcosh-Ur[X/2] 
dt 

tr[X] < 2 cosh 
1 



< 



< 



uo - 







dtj 


{[ 


UjQ 









W < - { cosh 



LOq 



^0 



dt - u , 



= sinh 









dtj 




i: 


CJo 











(8.4.14) 
(8.4.15) 

(8.4.16) 

(8.4.17) 
(8.4.18) 

(8.4.19) 

(8.4.20) 
(8.4.21) 



1/3^ = - {tr [T T^] - 2} = - {tr {X} - 2} , (8.4.22) 



(8.4.23) 
(8.4.24) 
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So finally 











{i 




LUq 




Jti 







(8.4.25) 



and consequently 



lal < cosh 









dtj 






UJQ 




hi 







(8.4.26) 



These bounds are quite remarkable in their generality. A version of this 
result was derived in [88| but the present derivation is largely independent 
and has the virtue of being completely elementary — in particular, the use 
of complex numbers has been minimized, and we have absolutely eliminated 
the use of the "auxiliary functions" and "gauge conditions" that were needed 
for the derivation in [HE]. If one translates this to the space domain, then 
the equivalent barrier penetration coefficient is Ttransmission ^ VI'^P; the 



equivalent reflection coefficient is R ^ 
translations 



\a\ 



Making the appropriate 



transmission 



> sech 











ko 


dx j> 







(8.4.27) 



and 



R < tanh^ 





k\x) 


da; I 


m: 


ko , 
ko 







(8.4.28) 



Comment: For completeness we mention that reference [88] provides 
a large number of consistency checks on these bounds by comparing 
them with known exact results |117] . 
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8.5 Lower bound on 



To obtain a lower bound on the Bogoliubov coefficient, consider any real 
valued parameter e. Then since the matrix T is itself real, 

tr{(T-eT^)^ (T-eT^)} > 0, (8.5.1) 

so that 

(1 + e^) tr(r T^) - 2e ti{T^) > 0, (8.5.2) 

whence 

tr(r^r)>^^tr(r2), (8.5.3) 

This bound is extremized for e = ±1, whence 

tr(r^r) > |tr(r2)| , (8.5.4) 

and so 

W>\{\tT{T')\-2}. (8.5.5) 

This is certainly a bound, but it is not as useful as one might hope. It is 
useful only if tT[T^] > 2. But 

ti[T^] = a^+d'^+2bc = a^+(f+2{ad-l) = = (tr[T])^-2. (8.5.6) 

So using the unit determinant condition, tr[T^] > 2 can be seen to require 
|a + d| > 2, that is, tr[T] > 2. But when does this happen? For the real 
matrix 

a b 
c d 

with unit determinant the eigenvalues are 



(8.5.7) 



X='-±^± ^S+SEI. (8.5.8) 

The condition a + d> 2 is thus equivalent to the condition that the eigenval- 
ues are real. Unfortunately there seems to be no simple way to then relate 
this to the properties of the function u!{t). 
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8.6 A more general upper bound 



Now let Q{t) be an arbitrary everywhere real and nonzero function of t 
with the dimensions of frequency. Then we can rewrite the Schrodinger 



ODE (8.2.1) as: 



d 










' 















(8.6.1) 



Again all the matrix elements have been carefully chosen to have the same 
engineering dimension. Again we can formally solve this in terms of the 
time-ordered product: 



The new T matrix is 
T = ^ J exp 



dt 



dt 



to 

8.6.2) 



(8.6.3) 



Note that the matrix T is still real, and that because 



tr 



-uj\i)/Q -Utt/Q) 



0. 



(8.6.4) 



it still follows that T has determinant unity: 



a b 
c d 



ad — be = 1. 



(8.6.5) 



This means that much of the earlier computations carry through without 
change. In particular as long as at the initial and final times we impose 
fl{t) — i> cuo as t — > and t — >• tt, we still have 



1 



a 



[a + d + i{b- c)] exp{-iLUo[tf - tj]), 



(8.6.6) 
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f3 ^ -[a-d + i{b + c)] exp(-ia;o[t/ + ti\), 
m' = \tr{TT^-l}, 
lap = ]tr{TT^ + 1]. 



(8.6.7) 
(8.6.8) 
(8.6.9) 



Now consider the quantity 

x{t) = T{t) T{tf = 



X 



^ < cxp 



-a;2(f)/f] - 1(0/1^) 



dt 



^.6.10) 



This now satisfies the differential equation 



dX 



x + x 



with the boundary condition 



and 



X{ti) = I, 



tr[X] = + 6^ + + (f. 
A brief computation yields 

+ 6^ ac + fed 



(8.6.11) 

(8.6.12) 
(8.6.13) 



dX 



+ 



|(n/n) n 

-u}^{i)/n -^in/Q) 

0? + fe^ ac + fee/ 

ac + fed c^ + d^ 



ac + fee? c^ + c?^ 



(8.6.14) 



{SllQ){c? + fe2) + 2^{ac + fed) l](c2 + d2) _ (o;7l])(a2 + fe2) 

-(a;Vl])(a2 + fe^) + ^(c^ + d^) -(2a;Vfi)(ac + fed) - (r2/l])(c2 + d^) 

(8.6.15) 
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Then taking the trace, there is now one extra term 



dtr[X] 
dt 



2{ac + bd) 



UJ 



(8.6.16) 



Note that if Q{t) — > loq then Q — > and we recover the ODE of the "ele- 
mentary" bound. In this more general setting we now proceed by using the 
following facts: 



As previously we note 
\ac + bd\^ + \ad-bc\^ 



a^c^ + 2abcd + b'^d^ + a^d^ - 2abcd + b'^c^, 
{a' + b'){c' + d% 

(8.6.17) 



which implies 



that is 



\ac + bd\ = ^/(^^2:rP)pT^2y3i^ 



2\ac + bd\ = V4(a2 + 62)(c2 + d2) _4_ 
• Additionally, we use 



(8.6.18) 
(8.6.19) 



|a' + 6' - - d^\ = ^J\a^ + 62 + c2 + rf2|2 _ 4(^2 + 52)^^2 + d^), 

(8.6.20) 

implying 

1^2 + 6^ _ c2 - d2|2 ^ (2|ac + bd\f = la^ + 6' + + - 4- (8.6.21) 

In particular, combining these observations, this means that we can find an 
angle 9 (which is in general some complicated real function of a, b, c, d) such 
that 

2(ac + 6d) = V|a2 + 62 + c2 + d2|2_4 gin^, (8.6.22) 
+ 6^ - - = V|a2 + 62 + c2 + d2|2 _4 cos 9, (8.6.23) 
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whence 

dtr[X] 
dt 



^y\a^ + b^ + c^ + d^\^-A {sine 



cos 6 



(8.6.24) 

But for any real 9 we certainly have by the Cauchy-Schwartz inequality 



sva.9 

implying 

dtr[X] 
dt 

Therefore 



+ cos 6 



< 



+ 



n 



1l 



< V|a2 + &2 + c2 + d2|2_4 







2 
















+ 











dtr[X] 
dt 



< Vtr[X]2 _ 4 







2 




2 






+ 










S2. 





implying 



whence 



so that 



1 dtr[X] 
^tr[X]2 - 4 dt - \ 

dcosh-^(tr[X]/2) 









2 


r 


2 








+ 


n- — 








Q 


n 





dt 



< 







2 




2 






+ 








Q 


n 





ir[X\ =tr[rr^] < 2 cosh 



+ 



(8.6.25) 



(8.6.26) 



(8.6.27) 



(8.6.28) 



(8.6.29) 



dt } . (8.6.30) 



Using the general formulae for jap and in terms of trlTT^}, and sim- 
plifying, we see 



< 



sinh^ji^ ' _1_^Q2 + [Q2_^2]2^^|^ (8.6.31) 
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and 



|a|'<cosh2|^^^*'py^S^Ti^?^^dt|. (8.6.32) 

This result is completely equivalent to the corresponding result in |Blj ; though 
again note that the derivation is largely independent and that it no longer 
requires one to introduce any "gauge fixing" condition, nor need we intro- 
duce any WKB-like ansatz. The current proof is much more "direct" , and at 
worst uses simple inequalities and straightforward ODE theory. If we work 
in the space domain instead of the time domain and make the translations 
Q{t) — >• (p'{x), uj{t) — >• k{x), we see 

lap <cosh2|^ f ^ J^^(<^")2 + [(^/)2_^2]2^2;|, (8.6.33) 

and 

< sinh^ 1^ \^\\/^^"^^ + t^*^')' ~ "^"^l • ^^■^■2^) 

This is perhaps physically more transparent in terms of the equivalent trans- 
mission and refiection coefficients 

-j-J {^"y + [{^'y - k^f dx\ , (8.6.35) 



-^transmission ^ SCch 

and 



1 1 

2 A, b 



i?<tanh^|^ f ^ -^^J{ip"y + [{ip'y -k^f dx\ . (8.6.36) 



Comment: For completeness we mention that reference |i88j provides 
a number of consistency checks on these more general bounds by com- 
paring them with known exact results |117] . 
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8.7 The "optimal" choice of Q(t)? 

What is the optimal choice of Q{t) that one can make leading to the most 
stringent bound on the Bogohubov coefficients? The bound we have just 
derived holds for arbitrary f2(t), subject to the two boundary conditions 
Q{ti) = ujq = Q{t f) and the overall constraint Q(t) ^ 0. Since both sinh and 
cosh are convex functions, finding the most stringent constraint on \f3\ and 
I q; I is thus a variational calculus problem equivalent to minimizing the action 

S= r -L Jq2^^^2_ ^2]2 (871) 

Jti l^'l 

The relevant Euler-Lagrange equations are quite messy, and progress (at 
least insofar as there is any practicable progress) is better made by using an 
indirect attack. The Lagrangian is 

1 



,'^' + [^'-ujr, (8.7.2) 
and so the corresponding canonical momentum can be evaluated as 

- , " (8.7.3) 

Prom the boundary conditions we can deduce 

n{t,)^^^n{tf). (8.7.4) 

The Hamiltonian is now 



|Q| ;/Q2 + [Q2_^2]2 |Q| ^Q2 + [Q2_^2]2 

(8.7.5) 



Unfortunately the Hamiltonian is explicitly time-dependent [via u!{t)] and 
so is not conserved. The best we can say is that at the endpoints of the motion 

H{ti)^0^ H{tf). (8.7.12) 
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Hamiltonian (classical mechanics): For a closed system the sum 
of the kinetic and potential energy in the system is represented by a 
quantity called the Hamiltonian, which leads to a set of differential 
equations known as the Hamilton equations for that system [TU7J. 
The Hamilton equations are generally written as follows: 

p ^ (8.7.6) 

, - f . (3.7.7) 

the dot denotes the ordinary derivative with respect to time of the 
functions p = p{t) (called generalized momenta) and q = q{t) (called 
generalized coordinates), and 7i = T-C{p,q,t) is the so-called Hamilto- 
nian. Thus, a little more explicitly, one can equivalently write 

^p(t) = -^nipit),q{t),t), (8.7.8) 

= ^n{p{t),q{t),t), (8.7.9) 

and specify the domain of values in which the parameter t ("time") 
varies. 



By solving for f2 as a function of vr and Q we can also write 
and 

(8.7.14) 

Note that Q at the endpoints cannot in general be explicitly evaluated in 
terms of the boundary conditions. 

An alternative formulation which slightly simplifies the analysis is to 
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Hamiltonian (quantum mechanics): The Hamihonian H is the 
quantum mechanical operator corresponding to the total energy of the 
system. It generates the time evolution of quantum states. If \ip{t)) 
is the state of the system at time t, then |1U6] 



where h is the reduced Planck constant. This equation is known as 
the (time dependent) Schrddinger equation. (This is the same form 
as the Hamilton- Jacobi equation, which is one of the reasons H is also 
called the Hamiltonian.) Given the state at some initial time (t = 0), 
we can integrate it to derive the state at any subsequent time. If H is 
independent of time, then 



In this thesis we are sometimes using quantum Hamiltonians, and 
sometimes classical Hamiltonians. Whether we are in a quantum or 
classical situation will have to be determined from context. 



change variables by writing 



i/|^(t)) = zn-|^(t)), 






n{t) = ujQ exp[^(t)] 




where the boundary conditions are now 





and the action is now rewritten as 



S 
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Then, in terms of this new variable we have 



-2(9 



with (dimensionless) conjugate momentum 



TT 



e 



=2e _ 



-2$ 



and boundary conditions 



7r(t0 = l = 7r(i/). 
The (non-conserved) Hamiltonian is 



e2e_ 4 e-2e 



:,26» _ tid p-2e 



subject to 



Inverting, we see 



TT 



a;o 



.261 



a; 



-26 



and 



29 _ ^ -20 
, ,2 



(8.7.18) 



(8.7.19) 



(8.7.20) 



(8.7.21) 



(8.7.22) 



(8.7.23) 



(8.7.24) 



This has given us a somewhat simpler variational problem, unfortunately the 
Euler-Lagrange equations are still too messy to provide useful results. 

Overall, we see that while solving the variational problem would indeed 
result in an optimum bound, there is no explicit general formula for such a 
solution. In the tradeoff between optimality and explicitness, we will have to 
accept the use of sub-optimal but explicit bounds. 
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8.8 Sub-optimal but explicit bounds 

From our general bounds 



'2 < Sinh^ l^j^' |^^fi2+ [f]2 _^2]2 



and 



lap < cosh^ 



2 A Ifil 

the following special cases are of particular interest: 



18.2) 



Q = uq: In this case we simply obtain the "elementary" bound considered 
above. 

Q = u: This case only makes sense if cj^ > is always positive. (Otherwise 
uj and hence Q becomes imaginary in the "classically forbidden" region; 
the matrix T then becomes complex, and the entire formalism breaks 
down). Subject to this constraint we find 



< sinh^ 



and 



|ap < cosh^ 
This case was also considered in 







Co 


dtj 










hi 










UJ 






1" 




dtj 




Jti 


UJ 





(8.8.3) 
(8.8.4) 



^ = uj'" uj\ ^: This case again only makes sense if cj^ > is always positive. 
Subject to this constraint we find 



< sinh^ 



and 



la I < cosh 



2 



Cj2 UJ 



2e 



+ 



UJ^ 



UJ, 



2-2e 




dt 



15) 



cl;2 cu2^ 



+ 



ujI"^' - a;2-2. 



UJ^ 



OJ, 



2-2e 




dt \ . 



This nicely interpolates between the two cases given above, which cor- 
respond to e = and e = 1 respectively. 
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Comment: The WKB approximation technique provides an approx- 



imate solution to the Schrodinger equation (8.2.3) of a quantum me- 
chanical system; however this technique fails at the classical turning 
points of the system of the potential energy function, V{x) 

Standard WKB theory provides a set of "joining conditions" for "step- 
ping over" the turning points and penetrating into the classically for- 
bidden region. Unfortunately we do not seem to have any analogue of 
these "joining conditions" in the formalism we set up in this thesis. 



Triangle inequality: Since y/ x"^ + y"^ < \x\ + \y\ we see that 
and 



lap < cosh^ 



These bounds, because they are explicit, are often the most useful quantities 
to calculate. 





/■*/ 


(I 






/•*/ 












dt ^ 








dt j 




Jti 




' 2 


Jti 
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8.9 The "interaction picture 



Interaction picture (sometimes called the Dirac picture): This 
is an intermediate between the Schrodinger picture and the Heisenberg 
picture. Whereas in the other two pictures either the state vector or 
the operators carry time dependence, in the interaction picture both 
carry part of the time dependence of observables. Equations that in- 
clude operators acting at different times, which hold in the interaction 
picture, do not necessarily hold in the Schrodinger or the Heisenberg 
picture. This is because time-dependent unitary transformations re- 
late operators in one picture to the analogous operators in the oth- 



For our purposes the interaction picture is useful because it lets us 
develop a perturbation theory. 



If we split the function uj(t)'^ into an exactly solvable piece uJe(tY and a 
perturbation uj/^{t)'^ then we can develop a formal perturbation series for the 
transfer matrix T, in close analogy to the procedures for developing quantum 
field theoretic perturbation theory in the interaction picture. Specifically let 
us write 



ers [inn]. 




and 



dT(t) 





dt 



Now defining 



T(t)=Te(t)TA(t), 




we shall develop a formal solution for T^it). Consider 



dT(t) 





dt 
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and compare it with 

dT{t) dTe(t) 



dt 

Therefore 
whence 



dt 



TaW + Te{t) = Qe{t) %{t) T^it) + Te{t) "^^^^^^ 



dt 



dt 



dT^jt) 
dt 



[Teit)-^ QA{t)Te{t)} Ta, 



^.9.5) 
^.9.6) 



TAit) = ^exp ^ {Teit)-' QAit) Teit)}dt^ . (8.9.7) 



For the fuU transfer matrix T we have 

ft 



T{t) = Te{t) X ^exp ^ {Te{t)-' QA{t) Te{t)}di 



and we have succeeded in sphtting it into an exact piece Tg{t) plus a distortion 
due to Qa {t) . This can now be used as the starting point for a perturbation 
expansion. (The analogy with quantum field theoretic perturbation theory 
in the interaction picture should now be completely clear.) 

To develop some formal bounds on the Bogohubov coefficients it is useful 
to suppress (currently) unnecessary phases by defining 



a — -[a + d + i{b — c)], 
/3 = J [a - d + i (6 + c)] . 



(8.9.9) 
(8.9.10) 



The virtue of these definitions is that for T = Ta they satisfy a simple 
composition rule which can easily be verified via matrix multiplication. From 
T = Te Ta we have 



a b 




c d 





ttg a A + be Ca Oe ^A + &e C^A 
Cg Oa + de Ca Ce &A + de dA 



Then some simple linear algebra leads to 

P ^ die Pa + Pe <5a, 



(8.9.12) 
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(8.9.13) 



But then 



<5e + /3e 



< 



+ 



(8.9.14) 



that is 



< loiel |/3a| + \Pe\ |q;a| , 



or the equivalent 



m<\/l + \Pef \PA\ + \Pe\ VI + I/SaI'. 



(8.9.15) 



^.9.16) 



Similarly 



\P\ = \p\ = aePA+Pe&l 

that is 



> 



Q!e /^A 



/3e Q!a 



|/?| > I l^el \Pa\ - \Pe\ I^aI I , 



Pa\ - \Pe "aI I , 

(8.9.17) 
(8.9.18) 



or the equivalent 



> 



^l+\Pef \Pa\ - \Pe\ ^1 + \Pa\ 



(8.9.19) 



The benefit now is that one has bounded the Bogohubov coefRcient in terms 
of the (assumed known) exact coefficient Pe and the contribution from the 
perturbation Pa- Suitably choosing the split between exact and perturbative 
contributions to o;^, one could in principle obtain arbitrarily accurate bounds. 



8.10 Conclusion 

In this chapter we again considered rigorous bounds on transmission, re- 
flection, and Bogoliubov coefficients. In particular, the most outstanding 
features of this chapter are: 
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• We have re-considered the general bounds on the Bogohubov coeffi- 
cients developed in [88]. Additionally, we have seen how to extend the 
bounds in |SB] in many different ways. Moreover, we do not need to 
"gauge fix" , nor do we need to appeal to any "WKB-like ansatz" to get 
the discussion started. 

• We have formulated some rigorous bounds that we can place on barrier 
penetration probabilities, or equivalently on the Bogoliubov coefficients 
associated with a time-dependent potential. Furthermore, we have not 
seen anything like these bounds anywhere else. 

• In addition, probably there are "optimal" bounds still waiting to be 
discovered. 

• It is apparent that the current bounds are not the best that can be 
achieved, and we strongly suspect that it may be possible to develop 
yet further extensions to the current formalism. 

Considering the fundamental importance of the questions we are asking, 
it is remarkable how little work on this topic can currently be found in the 
literature. 

Possible extensions might include somehow relaxing the reality constraint 
on Q(t) without damaging too much of the current formalism, a better under- 
standing of the variational problem defining the "optimal" bound (thus hope- 
fully leading to an explicit form thereof), or using several "probe functions" 
[instead of the single function ^l{t)] to more closely bound the Bogoliubov 
coefficients. 
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Chapter 9 



Bounding the greybody factors 
for Schwarzschild black holes 

9.1 Introduction 

Black hole greybody factors are important because they modify the spectrum 
of Hawking radiation seen at spatial infinity [71], so that it is not quite 
Planckian [75]. (That is, it is no longer exactly "blackbody radiation", which 
is why [with slight abuse of language], it is called "greybody".) There is a 
vast scientific literature dealing with estimates of these black-hole greybody 
factors, using a wide variety of techniques [90] . 

Unfortunately, most of these calculations adopt various approximations 
that move one away from the physically most important regions of parameter 
space. Sometimes one is forced into the extremal limit, sometimes one is 
forced to asymptotically high or low frequencies, sometimes techniques work 
only away from (3+1) dimensions, sometimes the nature of the approximation 
is uncontrolled. As a specific example, monodromy techniques fail for s = 1 
(photons) [77] , which is observationally one of the most important cases one 
would wish to consider. 

Faced with these limitations, in this chapter we ask a slightly different 
question: Restricting attention to the physically most important situations 
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(Schwarzschild black holes, (3+1) dimensions, intermediate frequencies, un- 
constrained spin and angular momentum) is it possible to at least place 
rigorous (and hopefully simple) analytic bounds on the greybody factors? 

By now considering the Regge- Wheeler equation for excitations around 
Schwarzschild spacetime, and adapting the general analysis discussed in ear- 
lier sections of this thesis, and published in references [SHI [^, we shall 
demonstrate that rigorous analytic bounds are indeed achievable. While it is 
certainly true that these bounds may not answer all the physical questions 
one might legitimately wish to ask, they are definitely a solid step in the 
right direction. 

Before starting the detailed calculations for this chapter, we should stress 
some important issues related to the greybody factors in Schwarzschild black 
holes: 

• Hawking radiation was originally derived in the geometric optics ap- 
proximation where it can be shown to be described by ideal black body 
radiation — a black body is an object that absorbs all light that falls 
on it. No electromagnetic radiation passes through it and none is re- 
flected. Because no light is reflected or transmitted, the object appears 
black when it is cold |121j . This derivation led to a calculation of the 
temperature and entropy of a black hole. 

• In this chapter, we shall try to rigorously solve for bounds on the trans- 
mission probabilities for waves moving through the region of space out- 
side of a Schwarzschild black hole — Schwarzschild black holes are the 
simplest type, one that is described as a spherically symmetric body 
with no (electric or magnetic) charge or angular momentum (no rota- 
tion). 

• The most important issue is to realise that "greybody factor" is actually 
a synonym for "transmission probability". Indeed, the phrase "grey- 
body factor" is used more in the thermodynamics and spectroscopy 
communities, while the phrase "transmission probability" is used more 
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in the quantum mechanics community, but these communities are re- 
ferring to the same concept. 

• The greybody factor describes the emissivity of the black hole which is 
not that of a perfect blackbody. 

• As a black hole radiates energy by Hawking radiation, energy conser- 
vation implies that it will lose mass. 



Schwarzschild black hole or static black hole: This black hole 
has no charge or angular momentum. A Schwarzschild black hole is 
described by the Schwarzschild metric, and cannot be distinguished 
from any other Schwarzschild black hole except by its mass. This 
black hole is characterized by a spherical surface, called the event 
horizon, which encloses the central singularity. The event horizon is 
situated at the Schwarzschild radius, often called the radius of a black 
hole. 

Any spherically symmetric mass distribution whose radius is smaller 
than the Schwarzschild radius forms a black hole |120j . 



9.2 Hawking radiation 

A black body is an object that absorbs all light that falls on it. No electro- 
magnetic radiation passes through it and none is reflected. Because no light 
is reflected or transmitted, the object appears black when it is cold |121j . 
The light emitted by a black body is called black body radiation. 

Hawking radiation is an approximately thermal radiation with an approx- 
imately black body spectrum predicted to be emitted by black holes due to 
quantum effects. It was discovered by Stephen Hawking who provided the 
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Evaporating black holes: There are several significant consequences 
that appear as a result of considering evaporating black holes: 

(1) We first note that black hole evaporation generates a more consis- 
tent view of black hole thermodynamics, by showing how black holes 
interact thermally with the remainder of the universe. 

(2) Secondly, the temperature of a black hole increases as it radiates 
away mass. In fact, the rate of temperature rises as explosively, with 
the most plausible of the endpoint scenarios being the complete disso- 
lution of the black hole in a violent explosion of gamma rays. To make 
a complete description of this assumed dissolution process, one would 
need a model of quantum gravity, which would be needed in the fi- 
nal stages of the evaporation process, when the black hole approaches 
Planck mass and Planck radius. (The precise details of what is go- 
ing on here is the subject of much continued, and sometimes heated, 
debate.) 

(3) Finally, the simplest models of black hole evaporation lead to the 
black hole information paradox. For instance, the "information con- 
tent" of a black hole appears to be lost when it evaporates, as under 
many of these models the Hawking radiation is purely random. It 
is felt by many (not all physicists) that the Hawking radiation must 
be somehow perturbed to contain the missing information, that the 
Hawking evaporation somehow carries the missing information to in- 
finity. This would imply that information is not allowed to be lost 
under these conditions |122j . (Technically, in this model the Hawking 
evaporation process is asserted to be "unitary" . The precise details of 
what is going on here is the subject of much continued, and sometimes 
heated, debate [931191].) 
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Entropy: Black holes are truly thermodynamic systems with an ac- 
tual and precisely calculable temperature and entropy. However, in 
terms of statistical mechanics, the entropy should be the logarithm of 
the number of independent states of the black hole. Understanding 
how to count these states would constitute a significant progress in 
the probe to understand quantum gravity [T24J. (In this regard, some 
string-based models have led to at least partial success.) 



theoretical argument for its existence in 1974, and is closely related to work 
by the physicist Jacob Bekenstein who predicted that black holes should have 
a finite, non-zero entropy |122j . 

Greybody factors in black hole physics modify the naive Planckian spec- 
trum that is predicted for Hawking radiation when working in the limit of 
geometrical optics. We shall consider the Schwarzschild geometry in (3+1) 
dimensions, and analyze the Regge-Wheeler equation for arbitrary particle 
spin s and wave-mode angular momentum i, deriving rigorous bounds on the 
greybody factors as a function of s, i, wave frequency uj, and the black hole 
mass m. 

9.3 Regge— Wheeler equation 

The well-known Regge-Wheeler equation describes the axial perturbation 
of Schwarzschild metric in the linear approximation [93]. In terms of the 
tortoise coordinate the Regge-Wheeler equation (Gjv ^ 1) is 



[u' - V{r)]ij, 



(9.3.1) 



where for the specific case of a Schwarzschild black hole 



dr 
dr^ 



1 



2m 



r 



(9.3.2) 
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Comment: The most standard interpretation of the Hawking radia- 
tion theory states that imphcit "virtual" particle-antiparticle pairs are 
occasionally created outside the event horizon of a black hole. There 
are three things that can happen to a pair of particles just outside the 
event horizon [125j : 

• Both particles are pulled into the black hole. 

• Both particles escape from the black hole. 

• One particle escapes while the other is pulled into the black hole. 
In particular, the particle that has escaped becomes real and 
can consequently be observed from Earth (or by any outside ob- 
server). The particle that was pulled into the black hole remains 
virtual, and to satisfy conservation of energy must have negative 
mass-energy. The black hole absorbs this negative mass-energy 
and as a result, loses mass and appears to shrink. The total 
rate of power emission is (to an excellent approximation) pro- 
portional to the inverse square of the black hole's mass. 

This is only a visual picture, at best an aid to understanding, and many 
relativists would argue that the true situation can only be begun to 
be understood by calculating the renormalized stress-energy tensor in 
the vicinity of the horizon. Such calculations agree that the key point 
is that negative energy is flowing into the black hole. 
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Quantum field theory in curved spacetime: Ordinary quantum 
field theories are defined in flat Minkowski space, which is an outstand- 
ingly good approximation for expressing the physics of microscopic 
particles in weak gravitational fields like those seen on Earth. 
In contrast, one can formulate quantum field theories in curved space- 
time to describe situations in which gravity is powerful enough to 
infiuence quantum matter, while, however, it is not powerful enough 
to require quantization itself. Importantly, these theories depend on 
classical general relativity to express a curved background spacetime. 
A generalized quantum field theory is used to describe the behavior 
of quantum matter within that spacetime. In addition, one can use 
this formalism to show that black holes emit a blackbody spectrum of 
particles known as Hawking radiation, leading to the possibility that 
they evaporate over time. It is believed that this radiation represents 
a significant part of the thermodynamics of black holes 
For our purposes, we will not be dealing with curved-space quantum 
field theory directly. However, the general techniques we develop in 
this thesis can be adapted to answer specific questions in curved-space 
quantum field theory, such as placing limits on the transmission co- 
efficients in black hole scattering, and bounding cosmological particle 
production due to the expansion of the universe. 



and the Regge-Wheeler potential is 



Here s is the spin of the particle and i is the angular momentum of the 
specific wave mode under consideration, with i > s. Thus V{r) > outside 
the horizon, where r G (2m, oo). The greybody factors we are interested in 
are just the transmission probabilities for wave modes propagating through 




(9.3.3) 
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this Regge- Wheeler potentiaL 



Tortoise coordinate: In a Schwarzschild spacetime the so-called 
"tortoise coordinate" is defined by |126] : 



r + 2G'Mln 



2GM 



1 



(9.3.4) 



The tortoise coordinate r* approaches — oo as "r" approaches the 
Schwarzschild radius r = 2GM. It satisfies 

dr* 2GM\'^ 



dr 



1 . (9.3.5) 

r / 



If we watch an object fall towards the Schwarzschild radius, then in 
terms of the r coordinate the object seems to slow down and (asymp- 
totically) stop at the horizon. In contrast, 

dr* 

- (9.3.6) 

remains finite. In some vague analogy with the fable of the "tortoise 
and the hare", the r* coordinate has come to be known as the "tor- 
toise" coordinate. 



Despite comments often encountered in the literature, one can explic- 
itly solve for r as a function of the tortoise coordinate — in terms 
of the Lambert W function we have the exact result 



r 



2m [1 + Wie^''" 



whereas 



r*(r) = r + 2m In 



-2m] /2m 



2m 



2m 



(9.3.7) 



(9.3.8) 



Unfortunately this formal result, while certainly correct and exact, is 
less useful than one might suppose. (We have not been able to turn 
this observation into any useful calculation.) 
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• Despite other comments often encountered in the hterature, one can 
also exphcitly solve the Regge-Wheeler equation — now in terms of 
Heun functions [95J. Unfortunately this is again less useful than one 
might suppose, this time because relatively little is known about the 
analytical behaviour of Heun functions — this is an area of ongoing 
research in mathematical analysis ^B] . 



Heun functions: In mathematics, the Heun differential equation is 
a second-order linear ordinary differential equation of the form |127] 



+ 



7 6 

- + 

z z — 



z — d 



dw 
dz 



a(3z — q 



w 



0. (9.3.9) 



z{z — l){z — d) 

Note that the constraint e = Q; + /5 — 7 — 5 + lis needed to ensure 
regularity of the singular point at 00. Every second-order linear ODE 
in the complex plane (or on the Riemann sphere, to be more accurate) 
with four regular singular points can be transformed into this equation. 
This standardized equation has four regular singular points: 0, 1, d, 
and 00. 

Studying the properties of these functions is an area of ongoing re- 
search in mathematical analysis 
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9.4 Bounds 

The general bounds developed earlier in this thesis, and published in refer- 
ences [HSl |89], can, in the current situation, be written as 

,2 



r>sech'|y ■(9dr,|. (9.4.1) 

Here T is the transmission probability (greybody factor), and is the func- 
tion 

^^^'^''^ + 1-^;-'^-'^, (9.4.2) 

where, h is some positive function, h{r^:) > 0, satisfying the limits h{—oo) = 
h{+oo) = uj, which is otherwise arbitrary. Two different derivations of this 
general result, and numerous consistency checks, can be found in earlier 
chapters [s] and [s] of this thesis, and in references [HHl ES] . 

(These bounds were originally developed as a technical step when study- 
ing the completely unrelated issue of sonoluminescence |112] , and since then 
have also been used to place hmits on particle production in analogue space- 
times [H3] and resonant cavities [SI] , to investigate qubit master equations [H5] , 
and to motivate further general investigations of one- dimensional scattering 
theory [86].) For current purposes, the most useful practical results are ob- 
tained by considering two special cases: 



1. If we set /i = c<j then 



T > sech^ <! ^ / V{n) dn !> , (9.4.3) 



whence 



T > sech^ 



1 

2^ 



i{e + l) 2m(l-s2) 

^2 .y^3 



dr } . (9.4.4) 

1 2m L J 
Therefore, since the remaining integral is trivial, we obtain our first 
explicit bound: 

T>=ecl42^(i±ii±(i^|. (9.4.5) 
I 8tom I 



174 



9.4. BOUNDS 



That is: 

r>.ech^((i±I)!±(^l^|. (9.4.6) 

\ Sum J ^ ' 

Note that this bound is meaningful for all frequencies. This is sufficient 
to tell us that at high frequencies the Regge- Wheeler barrier is almost 
fully transparent, while even at arbitrarily low frequencies some nonzero 
fraction of the Hawking flux will tunnel through. A particularly nice 
feature of this first bound is that it is so easy to write down for arbitrary 
s and i. 

2. If we now set h = \JuP' — V, which in this case implicitly means that 
we are not permitting any classically forbidden region, then 



T > sech 



dr, \ . (9.4.7) 



Since for arbitrary s and i the Regge-Wheeler potential is easily seen 
to have a unique peak at which it is a maximum, this becomes 

T > sech^/lnf^H, (9.4.8) 



CO 



= sech^ jln 1^ V^!IE^ ) ^ , (9.4.9) 

which is easily seen to be monotonic decreasing as a function of V^eak- 
However calculating the location of the peak, and value of the Regge- 
Wheeler potential at the peak is somewhat more tedious than evaluat- 



ing the previous bound (9.4.5). Note that the present bound fails, and 
gives no useful information, once u"^ < V^eak, corresponding to a clas- 
sically forbidden region. More explicitly, the bound can be rewritten 
as: 

(9.4.10) 



The most interesting point of the study of black hole greybody factors 
and (once one moves into the complex plane), the closely related problem of 
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Monodromy theorem: The monodromy theorem is a significant re- 
suh regarding analytic continuation of a complex-analytic function to 
a larger set. The basic concept is that one can extend a complex ana- 
lytic function along curves — the process can be started in the original 
domain of the function and can be ended in the bigger set. In addi- 
tion, this theorem gives sufficient conditions for analytic continuation 
to give the same value at a given point regardless of the specific curve 
used to get there, so that the resulting extended analytic function 
would be well-defined and single- valued (128|. 



locating the quasinormal modes [771 EH EH], is a subject that has attracted 
a wide amount of interest. In particular, quasinormal modes are poles of 
the transmission coefficient and refiect the black hole's ringdown reaction 
to a perturbation f77]. Unfortunately, as a specific example, monodromy 
techniques fail for s = 1 (photons) [77], which is observationally one of the 
most important cases one would wish to consider. It is for this reason that 
we resort to our general bounds to extract as much information as possible. 
Let us now consider various sub-cases: 

• For s = 1 (ie, photons) the situation simplifies considerably. (Remem- 
ber, this is the case for which monodromy techniques fail [77].) For 

s = 1 we always have rpeak = 3m and 

V^.. ^ (9.4.11, 



Consequently, from (9.4.10) 



[54^W-£(£+l)]2 ■ ^^-^-'^^ 
In almost the entire region where this bound applies (tu^ > Vpeak) it is 



in fact a better bound than (9.4.5) above 
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For s = (ie, scalars) and i = (the s-wave), we have rpeak = 8m/3 
and 

27 



Consequently 



^ 4096u;W[1024^W-27] 
^^=°'^=° - [2048c.2^2 _ 27]2 • 



In a large fraction of the region where this bound applies it is in fact a 



better bound than (9.4.5) above. 

For s = but £ > 1 it is easy to see that throughout the black hole 
exterior, Vr G (2m, oo), we have 



/ N , 2m \ 
Vs=o,i>i{r) < ( 1-— j 



re2 



(9.4.15) 



which is the s = 1 potential with the replacement 1) — > £^ + £ + 1. 
This bound on the potential has its maximum at rpcak = 3m, implying 

V'peak,.=0/>1 < ^^^2 • (9.4.16) 

Therefore the monotonicity of the bound on the greybody factor implies 

^ 108^W[27a;W-(£^ + £+l)] 
^'='''^' ^ [54^W-(£2 + £+l)]2 ' ^9.4.17) 

(for uj, m, and i held fixed, and subject to s < £). 

For s > 1 it is easy to see that throughout the black hole exterior, Vr G 
(2m, oo), keeping i held fixed, we have V^>i(r) < Vs=i{r). Therefore 

^eak,s>l < ^eak,s=l- (9.4.18) 

Therefore the monotonicity of the bound on the greybody factor implies 

108.;W[27.;W-£(£+l)] 

[54a;2m2-£(^+l)]2 ' ^""^ 

(for u, m, and £ held fixed, and subject to s < £). 
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More generally, it is useful to define 



+ 1 



(9.4.20) 



Excluding the case (s, £) = (0, 0), which was explicitly dealt with above, 
the remainder of the physically interesting region is confined to the 
range e G (— l,+l/2]. Then a brief computation yields 

rpeak = 3m|l- ^ + 0(e2)|, (9.4.21) 

and 

£{£ + !) 2e 2^ 



In fact one can show that 

over the physically interesting range. (This bound on Vpeak is tightest 
for (s, i) = (0, 1), corresponding to e = +1/2, where it provides a better 
than 1% estimate, and becomes progressively weaker as one moves to 
e = —1.) This then implies 



J-{s/mm > r.n. nfn , i " (9.4.24) 



As always there is a trade-off between strength of the bound and the 
ease with which it can be written down. 

While this second set of bounds has required a little more case by case analy- 
sis, we should in counterpoint observe that this second set of bounds provides 
much stronger information at very high frequencies, where in fact 

T> l-0[\/peakW"1. (9.4.25) 

Unfortunately this second set of bounds is (because of details in the deriva- 
tion, see earlier chapters in this thesis, and |88l|89]) not capable of providing 
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information once the frequency has dropped low enough for the scattering 
problem to develop classical turning points — in other words a scattering 
problem with a classically forbidden region is not amenable to treatment 
using bounds of the second class considered above. For sufficently low fre- 



quencies, bounds of the form (9.4.5) are more appropriate, with 



T > C(exp{-l/cu}). (9.4.26) 
What we have not done, at least not yet, is to use the full generality implicit 



in equation (9.4.2). Subject to rather mild constraints, there is a freely 



specifiable function h{r^:) available that can potentially be used to extract 
tighter bounds. Work along these lines is continuing. 



9.5 Discussion 

The study of black hole greybody factors [20], and (once one moves into 
the complex plane), the closely related problem of locating the quasinormal 
modes [771 EH |98] , is a subject that has attracted a wide amount of interest 
in both the general relativity and particle physics communities. In addition, 
we wish emphasize some specific features in this chapter: 

• In this chapter, we have developed a complementary set of results — 
we have sought and obtained several rigorous analytic bounds that can 
be placed on the greybody factors. 

• Even though these bounds are not necessarily tight bounds on the 
exact greybody factors they do serve to focus attention on general and 
robust features of these greybody factors. Moreover they provide a new 
method of extracting physical information. 

• In the current formalism, (as opposed to, for instance, monodromy 
techniques [77]), it is obviously clear that one does not have to know 
anything about what is going on inside the black hole in order to ob- 
tain information regarding the greybody factors. This is as it should 
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be, since physically the greybody factors are simply transmission co- 
efficients relating the horizon to spatial infinity, and make no intrinsic 
reference to the nature of the central singularity. 

• Looking further afield, here should be no intrinsic difficulty in extending 
these results to Reissner-Nordstrom black holes, dilaton black holes, or 
to higher dimensions — all that is really needed is an exact expression 
for the Regge-Wheeler potential. 

Finally, it is perhaps more interesting to see if one can significantly improve 
these bounds in some qualitative manner, perhaps by making a more strategic 
choice for the essentially free function /i(r*). 
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The Miller-Good 
transformation 

10.1 Introduction 

In this chapter, we shall consider further topics of general interest in quantum 
physics, such as transmission through a potential barrier — the potential 
being a region in a field of force where the force exerted on a particle is such 
as to oppose the passage of the particle through that region |129] — and 
the (formally) closely related issue of particle production from a parametric 
resonance. 

We have already developed a rather general bound on quantum trans- 
mission probabilities, and in the previous chapter (chapter [o]) have applied 
it to bounding the greybody factors of a Schwarzschild black hole. In this 
current chapter we shall take a different tack — we shall use the Miller- 
Good transformation (which maps an initial Schrodinger equation to a final 
Schrodinger equation for a different potential) to significantly generalize the 
previous bound. Moreover, we shall see that the Miller-Good transformation 
is an efficient method whereby to to generalize the bound, to make it more 
efficient and powerful. 
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Parametric resonance: This effect called "parametric resonance" 
is related to the phenomenon of mechanical excitation and oscillation 
at certain frequencies (and the associated harmonics). This effect is 
different from regular resonance because it often exhibits instability 
phenomena. 

Parametric resonance occurs in a mechanical system when a system is 
parametrically excited and oscillates at one of its resonant frequencies. 
Parametric excitation differs from forcing since the action appears as 
a time varying modification on a system parameter. The classical 
example of parametric resonance is that of the vertically forced pen- 
dulum [TT8] . 

For our purposes we will always view parametric resonance in terms 
of a time-dependent modifcation of the oscillation frequency. 



where k{xy = 2'm[E -~V{x)]/h'^ . As long as V{x) leads to finite (possibly dif- 
ferent) constants V±oo on left and right infinity, then for E > ma.x{V+oo, V-oo} 
one can set up a one-dimensional scattering problem in a completely standard 
manner — see for example [371 SHI SSI ED EHl ED EH OT- The scattering 
problem is completely characterized by the transmission and reflection am- 
plitudes {t and r), though the most important aspects of the physics can 
be extracted from the transmission and reflection probabilities (T = |tp and 
R = |rp). Relatively little work has gone into providing general analytic 
bounds on the transmission probabilities, (as opposed to approximate esti- 
mates), and the only known result as far as we have been able to determine 
is this: 



10.2 Setting up the problem 



Consider the Schrodinger equation. 
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Theorem 10.1. Consider the Schrddinger equation (10.2.1). Let h{x) > 
be some positive but otherwise arbitrary once-differentiable function. Then 
the transmission probability is bounded from below by 



- + 00 



T > sech^ 



2h 



dx 



^10.2.2) 



To obtain useful information, one should choose asymptotic conditions on 
the function h{x) so that the integral converges — otherwise one obtains the 
true but trivial result T > sech^cxo = 0. (There is of course a related bound 
in the reflection probability, R, and if one works with the formally equiva- 
lent problem of parametric oscillations, a bound on the resulting Bogoliubov 
coefficients and particle production.) 

This quite remarkable bound was first derived in [HH!, with further dis- 
cussion and an alternate proof being provided in chapter [s] (and published 
in [89]). These bounds were originally used as a technical step when study- 
ing a specific model for sonoluminescence |112] . and since then have also 
been used to place limits on particle production in analogue spacetimes ^3] 
and resonant cavities |83j, to investigate qubit master equations [85], and 
to motivate further general investigations of one-dimensional scattering the- 
ory [HB] . Most recently, these bounds have also been applied to the greybody 
factors of a Schwarzschild black hole (see previous chapter, and the related 
publication [90]). 

A slightly weaker, but much more tractable, form of the bound can be 
obtained by applying the triangle inequality. For h{x) > 0: 



T > sech 



ln(/i)'| 



da; 



(10.2.3) 



Five important special cases are: 



1. If we take h = k^o, where k^o = lim^^^-i-oo k{x), then we have [88] [89] 



T > sech 



2k 



k'\ dx \ . 



10.2.4) 
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2. If we define fc+oo = ^^t^x-^+oo k{x) ^ k_oo = lim^^^-oo ^(a;), and take 
h{x) to be any function that smoothly and monotonically interpolates 
between k_oo and A;+oo5 then we have 



T > sech 



In 



+00 



+ 



+00 



h 



dx 



10.2.5) 



This is already more general than the most closely related result pre- 
sented in [HHIIHH]. 



3. If we have a single extremum in h{x) then 



T > sech 



In 



h h 

'^+00 '^—00 

"■ext 



+ 



+00 



dx 



10.2.6) 



This is already more general than the most closely related result pre- 
sented in 



4. If we have a single minimum in fc^(x), and choose h!^ = max{/c^, A^}, 
assuming fc^j^ < < k'^^, (but still permitting fc^j^ < 0, so we are 
allowing for the possibility of a classically forbidden region), then 



T > sech 



2 



A;. 



A2 



+ 



2A 



\A^-k^\ dx 



(10.2.7) 



A2>fe2 

This is already more general than the most closely related result pre- 
sented in IS 



5. If fc^(x) has a single minimum and < fc^j^ < k^^, then 



T > sech 



In 



'^+00^—00 

'^min 



10.2.8) 



This is the limit of (10.2.7) above as A ^ k^[^ > 0, and is one of the 
special cases considered in [88] . 



In this chapter we shall not be seeking to apply the general bound ( 10.2.2 ), 



its weakened form (10.2.3), or any of its specializations as given in (10.2.4)- 
(10.2.8) above. Instead we shall be seeking to extend and generalize the 
bound to make it more powerful. The tool we shall use to do this is the 
Miller-Good transformation ISj 
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10.3 The Miller— Good transformation 



Consider the Schrodinger equation ( 10.2.1 ), and consider the substitution [87i 



= ^^jf=, U{X{x)). (10.3.1) 

We will want X to be our "new" position variable, so X{x) has to be an 
invertible function. This implies (via, for instance, the inverse function the- 
orem) that we need dX/dx 7^ 0. In reality, since the argument will be 
smoothest if we arrange things so that the variables X and x both agree as 
to which direction is left or right, we can without loss of generality assert 
dX/dx > 0, whence also dx/dX > 0. 

Now compute (using the notation Ux = dU/dX): 

u'ix) = Ux{X) 4X' - i^^l^ f/(X), (10.3.2) 



and 



1 X'" 3 iX''^^ 

u"{x) = UxxiX) {X'f" - 2 U + -^f^,U. (10.3.3) 



Insert this into the original Schrodinger equation, u{x)" + k{x)'^u{x) = 0, to 
see that 

' e _1 X'" 3 (X"^2 
{Jcy ~ 2(X')3 ^ IJX'Y 
which we can write as 

Uxx + K^U = 0, (10.3.5) 



Uxx + \y^,- ^y^, + 7 Wt^ }U = 0, (10.3.4) 



with 

(X'y [ 2 X' 4 (X')2 

That is, a Schrodinger equation in terms of u{x) and k{x) has been trans- 
formed into a completely equivalent Schrodinger equation in terms of U{X) 
and K{X). We can also rewrite this as 
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The combination 




1 X'" 3 {xy 



2 X' A (X')2 ' 



(10.3.8) 



shows up in numerous a priori unrelated branches of physics and is sometimes 
referred to as the "Schwarzian derivative" . 



Definition (Schwarzian derivative): The Schwarzian derivative 
is an absolute operator that is constant under all linear fractional 
transformations. For this reason, it commonly occurs in the theory of 
the complex projective line, and especially, in the theory of modular 
forms and hypergeometric series. 

The Schwarzian derivative also occurs in numerous other situations, 
including quantum field theory calculations related to stellar collapse 
and black hole formation, moving mirror calculations related to Unruh 
radiation, and much more. 

Explicitly, the Schwarzian derivative of a function / of one complex 
variable x is described by |130j 



The alternative notation {f,x} = {Sf){x) is frequently used. 



• As previously commented, to make sure the coordinate transformation 
X ^ X is well defined we want to have X'{x) > 0, let us call this 
j{x) = X'{x) with j{x) > 0. We can then write 




(10.3.9) 




(10.3.10) 
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Let us suppose that \im^^±^j{x) = j±oo 7^ 0; then = k±oo/j±oo, 
so if k'^{x) has nice asymptotic behaviour allowing one to define a scat- 
tering problem, then so does K'^{x). 

Another possibly more useful substitution (based on what we saw with 
the Schwarzian derivative) is to set J(a;)^^ = X'{x) with J{x) > in 
equation (10.3.6). When we let = X', we can then also write this 
as: 



We can then write 



X" 
X'" 



-2J'J-^ 
-2J" 6(J' 



/\2 



J3 



J4 



r 

y 



(10.3.11) 
(10.3.12) 

(10.3.13) 



Let us suppose that lima-^-too J{x) = J±oc 7^ 0; then K±oo = ^±oo<^±oo5 
so if fc^(x) has nice asymptotic behaviour allowing one to define a scat- 
tering problem, so does K'^{x). 



Calculation: We now substitute equations ( 10.3.11 )-( 10.3.12) above into 



(10.3.6), then we have 



1 X'" 3 (X 



/f\2 



iX') 



+ 



2 X' 4 (X') 

1 / -2J^^ Q{J'Y 
2V J3 + J2 



-2J' 
4V 1~ 



.'U-l(^.«^).^.?(^ 



J^lk^ 



.J" 3J'2 3J'2 

+ 



J 
r_ 

y 



J2 J2 



10.3.14) 



as required. 
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Two theorems: These observations about the behaviour at spatial infinity 
lead immediately and naturally to the result: 

Theorem 10.2. Suppose j ±00 = 1, {equivalently, J±oo = !)■ Then the "po- 
tentials" k'^{x) and K'^{X) have the same reflection and transmission ampli- 
tudes, and same reflection and transmission probabilities. 



This is automatic since K±oo = fc±oo, so equation (10.2.1) and the trans- 



formed equation (10.3.5) both have the same asymptotic plane-wave solu- 



tions. Furthermore the Miller-Good transformation (10.3.1) maps any linear 



combination of solutions of equation (10.2.1) into the same linear combina- 



tion of solutions of the transformed equation (10.3.5). QED. 



Theorem 10.3. Suppose j ±00 7^ 1, {equivalently, J±oo 7^ 1). What is the re- 
lation between the reflection and transmission amplitudes, and reflection and 
transmission probabilities of the two "potentials" fc^(x) and K'^{X)? This is 
also trivial — the "potentials" k'^{x) and ii'^(X) have the same reflection and 
transmission amplitudes, and same reflection and transmission probabilities. 

The only thing that now changes is that the properly normalized asymptotic 
states are distinct 

exp{ikoo x) exp{iKao x) 



;i0.3.15) 



but map into each other under the Miller-Good transformation. QED. 



10.4 Improved general bounds 



We already know 

2 r /'+°° 1 

T > sech < / i9 dx > . 
Here T is the transmission probability, and is the function 



(10.4.1) 



2h 



(10.4.2) 
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with h{x) > 0. But since the scattering problems defined by k{x) and K{X) 
have the same transmission probabihties, we also have 



T > sech 



+00 



^ dX 



with 



and 



dX = X' dx = J dx, 



2h 



-J(- 

2hy \x' 



1 



2hj 



J \ 2 J 4 p 



(10.4.3) 
(10.4.4) 

(10.4.5) 
(10.4.6) 

(10.4.7) 



That is: \/h{x) > 0, Vj(a;) > we now have (the first form of) the improved 
bound 



T > sech 



1 

2h 



ir,,_i/: + 3(^ 

J I 2 J 4 j2 



n 2 



dx 



(10.4.8) 

Since this new bound contains two freely specifiable functions it is definitely 



stronger than the result we started from, (10.2.2). The result is perhaps a 



little more manageable if we work in terms of J instead of j. We follow the 
previous logic but now set 



dX = X' dx = dx, 



(10.4.9) 



and 



7? 



^{hxY + [K^ - h^Y _ 1 



2h 



2h 



h' 
~X' 



+ 



J" 
T 



(10.4.10) 
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That is: V/i(x) > 0, \fJ{x) > we have (the second form of) the improved 
bound 



T > sech 



2h 



{h'Y + 



.P I A;2 



J" 
T 



/l2 
J2 



dx 



(10.4.11) 



A useful further modification is to substitute h = HJ^, then \/H{x) > 
0, VJ(x) > we have (the third form of) the improved bound by the follow- 
ing argument: When we let h = HJ^, then we derive 

h' = H'J^ + 2HJJ'. (10.4.12) 



We now substitute above equations into (10.4.11), implying that T is greater 
than 

"+CXD -j^ 



sech 



2HP 



{H'P + 2HJJ'y + 



J" 
T 



mp 



(10.4^13) 



We can simplify the above equation, and so we now also derive 

r>00 



T > sech 



1 
2H 



H' + 2H 



J' 
7 



J" 
u 



n 2 



dx > . (10.4.14) 



Equations (10.4.8), (10.4.11), and (10.4.14), are completely equivalent ver- 



sions of our new bound. 



10.5 Some applications and special cases 



We can now use these improved general bounds, (10.4.8), (10.4.11), and 



(10.4.14), to obtain several more specialized bounds that are applicable in 



more specific situations. 



10.5.1 Schwarzian bound 



First, take h = (constant) in equation (10.4.11), then 



T > sech^ 



J^l,2 J" 



h 

J2 



dx 



(10.5.1) 
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In order for this bound to convey nontrivial information we need the hmit 
hm^j^-i-oo J'^k^ = h"^, otherwise the integral diverges and the bound triviahzes 
to T > 0. The further speciahzation of this resuh reported in 



equation (10.2.4) above corresponds to J = (constant) 



^ and 

\/hjKx>, which 

clearly is a weaker bound than that reported here. In the present situation 
we can without loss of generality set h —>■ k^o it^ which case 

— -\- — \ - — 
koo \ J J 

We now need \imx^±oo J = 1 in order to make the integral converge. If 
k"^ > 0, so that there is no classically forbidden region, then we can choose 
J = ^/koo/k, in which case 



T > sech 



dx 



(10.5.2) 



T > sech 



Vk \Vk 



dx 



(10.5.3) 



Comment: This is a particularly elegant bound in terms of the 



Schwarzian derivative, [equation (10.3.8)], which however unfortu- 
nately fails if there is a classically forbidden region. This bound is 
also computationally awkward to evaluate for specific potentials. Fur- 
thermore, in the current context there does not seem to be any efficient 
or especially edifying way of choosing J{x) in the forbidden region, and 



while the bound in equation (10.5.2) is explicit it is not particularly 
useful. 



10.5.2 Low-energy improvement 



We could alternatively set H = (constant) in equation (10.4.14), to derive 

+ TTTTT k'^ + — -H'^ dx}. (10.5.4) 



T > sech 




J' 



J 



4iJ2 
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In order for this bound to convey nontrivial information we need to enforce 
linijj^-i-oo k"^ = k"^ = H'^i while hm^^^-i-oo J' = 0, and hma;^-|-oo J" = 0. Oth- 
erwise the integral diverges and the bound trivializes to T > 0. Thus 



T > sech 




J' 
7 



+ 



J" 



1 2 



dx 



Again, the further specialization of this result reported in 



(10.5.5) 



and equa- 



tion (10.2.4) above corresponds to J = (constant), which clearly is a weaker 
bound than that reported here. To turn this into something a little more 
explicit, since J{x) > we can without any loss of generality write 



J{x) = exp 



x{x) dx 



[10.5.6) 



where x{^) is unconstrained. This permits is to write 



T > sech 



'x^ + [k'^ + - X' - kl^f dx 



^10.5.7) 



Then by the triangle inequality 



T > sech 



\x\ + 



2K. 



2 I 

oo 



dx 



^10.5.8) 



A further application of the triangle inequality yields 



T > sech 



\X' 



+ 



k' + x'-k. 



dx 



(10.5.9) 



Now if k"^ < k'^, (this is not that rare an occurrence, in a non-relativistic 
quantum scattering setting, where k'^—k'^ = 2mV/ h'^ and we have normalized 
to Voo = 0, it corresponds to scattering from a potential that is everywhere 
positive), then we can choose = k'^ — k"^ so that 



T > sech 



2k^ 



\X 



dx 



(10.5.10) 
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Assuming a unique maximum for x (again not unreasonable, this corresponds 
to a single hump potential) this implies 



T > sech 



/ a/A;^ — /c^ da; > . 

J — oo I 



(10.5.11) 



This is a new and nontrivial bound, which in quantum physics language, 
where k'^ = 2m{E — V)/h'^, corresponds to 

' V2mV 



T > sech 



+ 



dx > . 



(10.5.12) 



If under the same hypotheses we choose x = 0, then the bound reported 



m 



and equation (|10.2.4|) above corresponds to 



T > sech 



1 



2VE J 



dx 



(10.5.13) 



Thus for sufficiently small E the new bound in equation (10.5.12) is more 



stringent than the old bound in equation (|10.5.13|) provided 



1 V2m V 

' max ^ o / ^ ■ 

2 h 



^10.5.14) 



Comment: Note the long chain of inequalities leading to these results 



this suggests that these final inequalities (10.5.11) and (10.5.12) 



are not optimal and that one might still be able to strengthen them 
considerably. 



10.5.3 WKB-like bound 



Another option is to return to equation (10.5.9) and make the choice = 
max{0, — /c^} = K^, so that n = \k\ in the classically forbidden region k"^ < 0, 
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while K = in the classically allowed region k"^ > 0. But then equation 



(10.5.9) reduces to 



T>sech2<; / «:dx + ^ + ^+ / dx } . (10.5.15) 

J l^oo 2 J 2k^ 

Key points here are the presence of /^2<o barrier penetration inte- 

gral that normally shows up in the standard WKB approximation to barrier 
penetration, Uraax the height of the barrier, and L the width of the barrier. 
These is also a contribution from the classically allowed region (as in general 
there must be, potentials with no classically forbidden region still generically 
have nontrivial scattering). Compare this with the standard WKB estimate: 

TwKB-sech^J j /«dx + ln2l. (10.5.16) 

I fc2<0 J 

This form of the WKB approximation for barrier penetration is derived, for 
instance, in Bohm's classic textbook [82], and can also be found in many other 
places. Under the usual conditions applying to the WKB approximation for 
barrier penetration we have /;j,2<o'^ dx ^ 1, in which case one obtains the 
more well-known version 

iwKB-exp^ -2 j Kdx). (10.5.17) 

fc2<0 



The bound in equation (10.5.15) is the closest we have so far been able to get 



to obtaining a rigorous bound that somewhat resembles the standard WKB 



estimate. Again we do not expect the bound in equation (10.5.15) to be 
optimal, and are continuing to search for improvements on this WKB-like 
bound. 

10.5.4 Further transforming the bound 



In an attempt to strengthen the inequalities (10.5.11) and (10.5.12), we again 



use the fact that J(x) > to (without any loss of generality) write J(x) 
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The Classically Forbidden Region: In regions where the total en- 
ergy is less than the potential energy, a region where in classical physics 
there would be zero probability of finding a particle, the amplitude of 
the wave function decreases. Further away from the boundary where 
the potential energy changes, the probability of the object being lo- 
cated there decreases. However, unlike classical physics, in the region 
where the total energy is less than the potential energy, the probability 
of finding the object is not always zero (131]. 



exp [J xi^) dx] , where x{^) is unconstrained. The general bound in equation 



(10.4.14) can then be transformed to: For all H{x) > 0, for all xi^) 



T > sech 



H' 



+ 2x 



(10.5.1. 



This leaves us with considerable freedom. Regardless of the sign of k'^{x), 
we can always choose to enforce /c^ + — = 0, and so eliminate either x 
or H, obtaining 



T > sech 



H' 



+ 



(V^2 _ 



dx 



;i0.5.19) 



(subject to H{x) > and H'^{x) — k'^{x) > 0), and 



T > sech^ 



Vx^T¥ 



-\ 2 



+ 2x 



X^ + k"^ 



^10.5.20) 



(subject to x^i.^) + fc^(x) > 0), respectively. Finding an explicit bound is 
now largely a matter of art rather than method. For example if we take 



= max{A;% A^} or x = max{0, - k^} 



(10.5.21) 
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then from either equation (10.5.19) or equation (10.5.20), again under the 
restriction that we are deahng with a single-hump positive potential, we 
obtain 

T > sech^ I ^ In + ~ ^'^"^^ + J VA^^ dx 

(10.5.22) 

Note that A is a free parameter which could in principle be chosen to optimize 
the bound, however the resulting integral equation is too messy to be of 
any practical interest. This bound is somewhat similar to that reported in 



equations (10.2.7) and (10.5.11), but there are some very real differences. 



10.6 Summary and Discussion 

The bounds presented in this chapter are generally not "WKB-like" — apart 



from the one case reported in equation (10.5.15) there is no need (nor does 
it seem useful) to separate the region of integration into classically allowed 
and classically forbidden regions. In fact it is far from clear how closely these 
bounds might ultimately be related to WKB estimates of the transmission 
probabilities, and this is an issue to which we hope to return in the future. 

We should mention that if one works with the formally equivalent problem 
of a parametric oscillator in the time domain then the relevant differential 
equation is 

u{t) + k{tY u{t) = 0, (10.6.1) 

and instead of asking questions about transmission amplitudes and proba- 
bilities one is naturally driven to ask formally equivalent questions about 
Bogoliubov coefficients and particle production. The key translation step is 
to realize that there is an equivalence 



1 1 - T 
tt; A^^^^. (10.6.2) 
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This leads to bounds on the number of particles produced that are of the 
form T > sech^{(some appropriate integral)}, thereby implying 



N < sinh^{(some appropriate integral)}. 



(10.6.3) 



To be more explicit about this, our new improved bound can be written in 
any of three equivalent forms: 

1. For all H{x) > 0, for all J{x) > 0, 



T > sech 



2H 



H' + 2H 



J' 
1 



+ 



J" 
u 



da; J . 
(10.6.4) 



2. For all h{x) > 0, for all J{x) > 0, 

f OO 1 



T > sech^ 



2h 



{h'f + 



J2 U2 + 



J" 1 



H 2 



J J2 



dx > . (10.6.5) 



3. For all h{x) > 0, for all j{x) > 0, 

"H-oo 



T > sech^ 



2h 



dx j . 
(10.6.6) 



The equivalent statements about particle production are: 
1. For all H(t) > 0, for all J{t) > 0, 



N < sinh^ 



J-oo '^H 



H' + 2H 



J' 
7 



+ 



I" 

tJ 



dt 



(10.6.7) 



2. For all h{t) > 0, for all J{t) > 0, 



N < sinh^ 



2h 



J" 1 /l2 



J J2 



dt > . ('10.6.^ 
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3. For all h{t) > 0, for all j(t) > 0, 



< sinh^ 




3(jO! 
2 J 4 f 



dt 
{10.6.9) 



In closing, we reiterate that these general bounds reported in equations 



(10.4.8), (10.4.11), and (10.4.14), their specializations in equations (10.5.2) 



(10.5.3), (10.5.11), (10.5.12), (10.5.15), and (10.5.22), and the equivalent par 



tide production bounds in equations (10.6.7)-( 10.6.9), are all general purpose 



tools that are applicable to a wide variety of physical situations [831 EH ESI 
861 1901ITT2] . Furthermore we strongly suspect that further generalizations of 
these bounds are still possible. 
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Chapter 11 

Analytic bounds on 
transmission probabilities 



11.1 Introduction 



In this chapter, we shall develop some additional and novel analytic bounds 
on transmission probabilities (and the related reflection probabilities and 
Bogoliubov coefficients) for generic one-dimensional scattering problems. We 
shall review the basic concepts underlying this fascinating topic by rewriting 
the Schrodinger equation for some complicated potential whose properties we 
are trying to investigate in terms of some simpler potential whose properties 
are assumed known plus a (possibly large) "shift" in the potential. Doing 
so permits us to extract considerable useful information without having to 
exactly solve the full scattering problem. 

In earlier chapters of this thesis, and in several published papers [SSI ESI 
[90l |9T] , we have derived a number of rigorous bounds on transmission prob- 
abilities (and reflection probabilities, and Bogoliubov coefficients) for one- 
dimensional scattering problems. The derivation of these bounds generally 
proceeds by rewriting the Schrodinger equation in terms of some equivalent 
system of ffist-order equations, and then analytically bounding the growth 
of certain quantities related to the net flux of particles as one sweeps across 
the potential. 
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CHAPTER 11. ANALYTIC BOUNDS ON TRANSMISSION 
PROBABILITIES 



In this chapter we shall obtain significantly different results, of both the- 
oretical and practical interest. While a vast amount of effort has gone into 
studying the Schrodinger equation and its scattering properties [23 SHI SSI 
EB EOl [m |82l [991 1100], it appears that relatively little work has gone into 
providing general analytic bounds on the transmission probabilities, (as op- 
posed to approximate estimates). The only known results as far as we have 
been able to determine are presented in [88], in the earlier chapters of this 
thesis, and the related publications [891 Effl El] based on work reported in 
this thesis. Several quite remarkable bounds were first derived in [HH], with 
further discussion and an alternate proof being provided in [89] . 




These bounds were originally used as a technical step when studying a 
specific model for sonoluminescence |112j . and since then have also been used 
to place limits on particle production in analogue spacetimes [83] and reso- 
nant cavities [M] , to investigate qubit master equations [HS] , and to motivate 
further general investigations of one-dimensional scattering theory [86]. Re- 
cently, these bounds have also been applied to the greybody factors of a 
Schwarzschild black hole [HD]- Most recently, significant extensions of the 
original bounds have been developed by adapting the Miller-Good transfor- 
mations [91] . 

In this chapter we shall return to this problem, developing a new set 
of techniques that are more amenable to the development of both upper 
and lower bounds. For technical reasons the new techniques are also more 
amenable to investigating behavior "under the barrier". The basic idea is 
to re-cast the Schrodinger equation for some complicated potential whose 
properties we are trying to investigate in terms of some simpler potential 
whose properties are assumed known, plus a "shift" in the potential. 
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11.2 From Schrodinger equation to system of 
ODEs 

We are interested in the scattering properties of the Schrodinger equation, 

i;"{x) + k{xYij{x) = 0, (11.2.1) 

where k{xy = 2'm[E — V{x)]/h'^. As long as V{x) tends to finite (possibly dis- 
tinct) constants V±oo on left and right infinity, then for E > max{V+oo, V^oo} 
one can set up a one- dimensional scattering problem in a completely stan- 
dard manner — see, for example, standard references such as [371 SOI SSI EH 
[80l EH |82l |99l 1100] , and the background discussion presented in earlier chap- 
ters of this thesis. The scattering problem is completely characterized by the 
transmission and reflection amplitudes (denoted t and r), although the most 
important aspects of the physics can be extracted from the transmission and 
reflection probabilities (T = |tp and R = |rp). 

11.2.1 Ansatz 

The idea is to try to say things about exact solutions to the ODE 

^"(x) + P{x) ij{x) = 0, (11.2.2) 

by comparing this ODE to some "simpler" one 

^;'(x) + A;2(x) V^o(x) = 0, (11.2.3) 

for which we are assumed to the know exact solutions ipoi^)- ^ manner 
similar to the analysis in references [88t 189], we will start by introducing the 
ansatz 

i/j{x) = a{x) ipoix) + b{x) iP*q{x). (11.2.4) 

This representation is of course extremely highly redundant, since one com- 
plex number iIj{x) has been traded for two complex numbers a{x) and b{x). 
This redundancy allows us, without any loss of generality, to enforce one 
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auxiliary constraint connecting a{x) and b{x). We find it particularly useful 
to enforce the auxiliary condition 

^,„+£,;^0. (11.2.5) 

Subject to this auxiliary constraint on the derivatives of a{x) and b{x), the 
derivative of ipi^x) takes on the especially simple form 

^ = a^', + bro- (11-2.6) 

(This ansatz is largely inspired by the techniques of references [HSIES], where 
JWKB estimates for the wave function were similarly used as a "basis" for 
formally writing down the exact solutions.) 



11.2.2 Probability density and probability current 

For the probability density we have: 

p = r^, (11-2.7) 

= \a{x)^o + b{x)^*o\\ (11.2.8) 
= M + \bmiJo\' + 2Re{ab*ijl}, (11.2.9) 

= {|a|2 + |6|2|}po + 2Re{a6*V'o}- (11.2.10) 

Furthermore, for the probability current: 

^ = Im|v^*^|, (11.2.11) 



Im I [a*ro + b*H Ho + &V^o1 |> (11.2.12) 

Im I \a\'roi^', + |6| VoV-o' + «&*^o^o + a*6^*^*'| , (11.2.13) 

{\a\'-mim{roi^o}, (11.2.14) 
{\a\'-m^o. (11.2.15) 
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Under the conditions we are interested in, (which correspond to a time- 
independent solution of the Schrodinger equation), we have p = 0, and so 
= 0. (And similarly po = 0, so d^^o = 0.) That is, ^ and are 
position-independent constants, an observation which then puts a constraint 
on the amplitudes \a\ and Applying an appropriate boundary condition, 
which we can take to be a(— oo) = 1, 6(— oo) = 0, we then see 

\a\^-\b\^ = l. (11.2.16) 

This observation justifies interpreting a{x) and b{x) as "position-dependent 
Bogoliubov coefficients" . Furthermore without any loss in generality we can 
choose the normalizations on ip and ipo so as to set the net ffuxes to unity: 
^ = ^0 = 1. 



11.2.3 Second derivatives of the wavefunction 

We shall now re-write the Schrodinger equation in terms of two coupled 
first-order differential equations for these position-dependent Bogoliubov co- 
efficients a{x) and b{x). To do this, evaluate d^'i/'/dx^ making repeated use 
of the auxiliary condition (11.2.5) 

g = ^{a^o + bro), (11.2.17) 
= a>; + 6>*' + aV^;' + 6V^*", (11.2.18) 

= a'^'Q-a'^ro'-akl^Po-bklro, (11.2.19) 
ro 

= ^{ro^'o-^oro'}-kl[a^o + bro], (11.2.20) 



-k'^[a^Q + br,i (11.2.21) 







2ia' 

= -^-kliaipo + bro]- (11.2.22) 

Where in the last line we have finally used our normalization choice = 1. 
This is one of the two relations we wish to establish. Now use the gauge 
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condition to eliminate da/dx in favour of db/dx to obtain a second relation 
for d'^ip/dx'^. This now permits us to write d'^tp/dx'^ in either of the two 
equivalent forms 

d'^ip 2ia' 



dx"^ 



2ib' 

-00 



- kl [a^o + bipl] 



;il.2.23) 
;il.2.24) 



11.2.4 SDE as a first-order system 

Now insert these formulae for the second derivative of the wavefunction into 
the Schrodinger equation written in the form 

d^^lj 



dx'^ 



+ k{xY ^ = 0, 



to deduce the pair of first-order ODEs 
da i 



dx 
db 

dx 



11.2.25) 

11.2.26) 
11.2.27) 



It is easy to verify that this first-order system is compatible with the auxiliary 



condition (11.2.5), and that by iterating the system twice (subject to this 
auxiliary condition) one recovers exactly the original Schrodinger equation. 
We can re-write this Ist-order system of ODEs in matrix form as 



d 

dx 



i[k 



l^ol 





a 




b 



11.2.28) 



11.2.5 Formal (partial) solution 

Define magnitudes and phases by 

a = \a\ e*"^"; b = \b\ e'^'; ?/^o = |^o| e'^° . (11.2.33) 
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Comment: Matrix ODEs of this general form are often referred to 
as Shabhat-Zakharov or Zakharov-Shabat systems [HB]. This matrix 
ODE can be used to write down a formal solution to the SDE in terms 
of "path-ordered exponentials" as in references [HHl EH] . We choose not 
to adopt this route here, instead opting for a more direct computation 
in terms of the magnitudes and phases of a and b. 



Calculate 

a' = \a\' e'^^ + i\a\ e*"^" 0'^ = e^"^" {\a\' + i\a\ <p'^} , (11.2.34) 

whence 

\a\' + i\a\ 01 = ^[fc' - kl] iV^ol' + \h\ e-'^^--^^+^^^^}. (11.2.35) 
Similarly we also have 

+ t\h\ 0; = -^[P - Q iV^ol' {\h\ + \a\ e-(<A'>-<A-2</,o)|, (11.2.36) 
Now take the real part of both these equations, whence 

l«r = + - kl] \h\ iV-or sin(0, - 0fe + 20o); (11.2.37) 



2' 

^1 



= - kl] \a\ |V^orsin(0,-0fe + 20o). (11.2.38) 



Therefore 



a 



2 [k^ - kl] l^ol' sin(0, - 0fc + 20o) V|a|2 - 1. (11.2.39) 



That is 

-l== = \[k'' -kl] |V'orsin(0,-0f. + 20o), (11.2.40) 

whence 

{cosh-i \a\]l[ = ]^j\e- kl] |^o|'sin(0, - 0^ + 20o) dx. (11.2.41) 
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The path-ordered exponentials (or the ordered exponential): 

In a general formal context, this mathematical object is defined in 
arbitrary non-commutative algebras, and it is the closest equivalent 
possible to the exponential function of the integral in the commutative 
algebras. 

In practice the values lie in matrix and operator algebras. For the ele- 
ment A{t) from the algebra {g* ) (the set g with the non-commutative 
product *), where t is the "time parameter" the ordered exponen- 
tial [132] . 

OE[A]{t):=exp( [ dt'A{t')]. (11.2.29) 







of A can be defined via one of several equivalent approaches: 

• As the limit of the ordered product of the infinitesimal exponen- 
tials: 



OE[A]{t) = lim <^ exp(ey4(itAr))xexp(eA(itjv_i))x- ■ ■xexp(eA(i:o)) 

(11.2.30) 

where the time moments {to,ti, . . . ,1^} are defined as tj = j x e 
for j = 0, . . . iV, and e = t/N. 

• Via the initial value problem, where the OE[A]{t) is the unique 
solution of the system of equations: 

^^3p)=A(t)x0i^[A](t), (11.2.31) 
where OE[A]{0) = 1. 

• Via an integral equation: 

OE[A]{t) = 1+ / dtiA(ti)+ / dti / dt2A{t,) X A{t2) 
Jo Jo Jo 

+ dti dt2 / dhAiti) X A{t2) X A{h) + ... 



^0 ^0 



(11.2.32) 
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Now apply the boundary conditions: At x — — oo we have both a(— oo) — 1, 
and 6 (—00) — 0. Therefore 

cosh-^ \a{x)\ = \ f [k^- kl] |^orsin(0„ -(P^ + 20o) dx, (11.2.42) 
^ J -00 



and so 



|a(a;)| =cosh|^y" [A;' - fco'] iV^ol' sin(0„ - 0^ + 20o) da;| . (11.2.43) 
In particular 

1 r+oo 

cosh-1 |a(oo)| = - / - kl] l^'ol' sin((/.„ - 06 + 2(/)o) dx, (11.2.44) 

or equivalently 

|a(oo)| ^coshj^y [A;'-A;g] |V'o|'sin(0„-0b + 2(/)o) dxj. (11.2.45) 

Of course this is only a, formal solution since (t)a{x) and (t)b{x) are, (at least at 
this stage) , "unknown" . But we shall argue that this formula still contains 
useful information. In particular, in view of the normalization conditions 
relating a and 6, and the parity properties of cosh and sinh, we can also 
write 



(11.2.46) 
(11.2.47) 



a(oo) 


= cosh 


2 . 


/•+00 

/ [k' 
1 —00 


-kl] 




(pb + 200 ) dx 








r+00 








|6(oo) 


= sinh 


2. 


/ [k' 
'—00 


-kl] 




06 + 200 ) da; 



11.2.6 First set of bounds 



To determine the first elementary set of bounds on a and h is now trivial. 
We just note that 

|sin(0„-06 + 20o)| < 1. (11.2.48) 



Therefore 



|a(oo)| <cosh|^y"^ |A;2 - iV'olMxj ; (11.2.49) 
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|6(oo)| < sinhj^y jfc^ - fc^j |7/>o|Mx| . (11.2.50) 
What does this now tell us about the Bogoliubov coefficients? 

11.2.7 Bogoliubov coefficients 

The slightly unusual thing, (compared to our earlier work in this thesis and in 
references [HHl ESI EI] ) , is that now the "known" function ipo may also have 
its own Bogoliubov coefficients. Let us assume we have set our boundary 
conditions so that for the "known" situation 

ipo{x ~ — oo) ~ exp{ik{—oo)x}, (11.2.51) 

and 

ipo{x ~ +oo) ~ aoexp{ik{+oo)x} + j3o exp{—ik{+oo)x} . (11.2.52) 
Then the way we have set things up, for the "full" problem we still have 

iIj^x ~ — oo) ~ exp{iA;(— oo)a;}, (11.2.53) 

whereas 

^{x^+oo) ~ a(oo) ^/'o(a;) + 6(oo) ^/'o(x), (11.2.54) 
~ [ao a(oo) + /3q b{oo)] exp{iA;(+oo)x} 

+ [(3o a(oo) + ao 6(oo)] exp{-2A;(+oo)x}. (11.2.55) 

That is, the overall Bogoliubov coefficients satisfy 

a = aoa(oo) +/3o&(oo); (11.2.56) 

p = Poa{oo) + alb{oo). (11.2.57) 

These equations relate the Bogoliubov coefficients of the "full" problem 
{ipix), k{x)} to those of the simpler "known" problem {i/jq^x), ko^x)}, plus 
the evolution of the a{x) and b{x) coefficients. Now observe that 

|a| < |ao| |a(oo)| + \Po\ |6(oo)|. (11.2.58) 
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But we can define 

|q;o| = coshOo; \Po\ — sinhOo; ["^(^o)! — cosh©; |&(oo)| = sinhO; 

(11.2.59) 

in terms of which 

\a\ <cosheocoshe + sinheosinhe = cosh(eo + e). (11.2.60) 
That is: Since we know 

2 r+oo 

e < Gbound =7, \k^- kl\ IM^dx, (11.2.61) 

we can deduce 

\a\ < cosh jcosh"^ |q;o| + - k^] iV'opdxj ; (11.2.62) 

\(3\ < sinh jsinh-^ \(3o\ + ^ J - kl\ \M'^dx^ . (11.2.63) 

11.2.8 Second set of bounds 

A considerably trickier inequaUty, now leading to a lower bound on the Bo- 
goliubov coefficients, is obtained by considering what the phases would have 
to be to achieve as much destructive interference as possible. That implies 

\a\ > \ao\ \a{oo)\ - \Po\ \b{°°)\, (11.2.64) 

whence 

|a| >cosh|eo-e|. (11.2.65) 

Therefore, using © < ©bound, it follows that as long as ©bound < ©o, one can 
deduce 

\a\> cosh {©0 - ©bound} ■ (11.2.66) 

(If on the other hand ©bound ^ ©o, then one only obtains the trivial bound 
\a\ > 1.) Another way of writing these bounds is as follows 



> cosh jcosh-^ laol - ^J \k'^ - kl\ |V'o|^dx| ; (11.2.67) 
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\(3\ > sinh jsinh-^ |/3o| J \k'^ - Q \ilJo\'^dx^ ; (11.2.68) 

with the tacit understanding that the bound remains vahd only so long as 
argument of the hyperbolic function is positive. 



11.2.9 Transmission probabilities 

As usual, the transmission probability (barrier penetration probability) is 
related to the Bogoliubov coefficient by 



r=^, (11.2.69) 



whence 



T > sech^ jcosh-^ laol + \J^ 1^^ - ^ol l^ol^d^^ ■ (11.2.70) 



That is 



T > sech=^ I cosh" ^ (To" ^/^) + - kl\ |^orda;| , (11.2.71) 



or even 

"+00 



T > sech^ |sech-^(ro/^) + - kl\ {il^o^dx'^ . (11.2.72) 

Furthermore, as long as the argument of the sech is positive, we also have 
the upper bound 

T < sech^ |sech-^(ro/^) -\J^ \k^ - ^ol iV'ol^dxj . (11.2.73) 

If one wishes to make the algebraic dependence on Tq clearer, by expanding 
the hyperbolic functions these formulae may be recast as the statement that 
T is greater than 

n 

cosh |i |P - IM^dx^ + vT^sinh \k^ - ki\ l^olMx}]" 

(11.2.74) 
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and (as long as the numerator is positive before squaring), that T is less than 

To 

cosh||/+^|A;2-fc2| l^olMx}- v^r^sinh|i/+^|P-A;2| |^o|Mx}]' 

(11.2.75) 

11.3 Consistency check 

There is one special case in which we can easily compare with the previous 
results of chapters [s] and [s] of this thesis, and references [HHIIHH]. Take ko = 
k{±oo) to be independent of position, so that our comparison problem is a 
free particle. In that case 



^0 = — 7= — ; m\ = -r'^ = 1? "0 = 1; /5o = o. 



11.3.1 



Then the bounds derived above simplify to 



1 



+ 00 



|a|<cosh<^— / |A:^ - A;^| da; ^ , (11.3.2) 

|/?| < sinhj-^ /"^ |fc2_^2|^^\ (11.3.3) 

This is "Case I" of reference [HE], and the "elementary bound" of refer- 
ence [M], which demonstrates consistency whenever the formalisms overlap. 
(Note that it is not possible to obtain "Case 11" of reference [HFJ or the 
"general bound" of reference [88t [89] from the present analysis — this is 
not a problem, it is just an indication that this new bound really is a dij- 
ferent bound that only partially overlaps with the previous results of refer- 
ences [HHIIMIEI].) 

A second (elementary) check is to see what happens if we set iIj{x) 
ipo{x), effectively assuming that the full problem is analytically solvable. In 
that case T — > Tq, (and similarly both a — >• ao and j3 — > Pq), as indeed they 
should. 
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11.4 Keeping the phases? 

We can extract a little more information by taking the imaginary parts of 



equations (11.2.35) and (11.2.36) to obtain: 



\b\ 

1 + f4 COS 

a 



% + 200, 



Subtracting 



<p,y = [e-kl] i^ol 



2\a\ \b\ 



cos 



h + 20o; 



:ii.4.i) 



fll.4.2) 



^11.4.3) 



This is now a differential equation that only depends on the difference in 
the phases — the overall average phase {(pa + 4>b)/'^ has completely decou- 
pled. (Moreover, in determining the transmission and reflection probabilities, 
this average phase also neatly decouples). To see how far we can push this 
observation, let us now define a "nett" phase 

A = 0, - 0b + 200- (11.4.4) 

Furthermore, as per the previous subsections, we retain the definitions 

\a\ = coshB; |6| = sinhG. (11.4.5) 



Then equation (11.2.43|) becomes 
e(a;) = 



[k - ko] |^orsin(A(x)) dx 



;ii.4.6) 



;il.4.7) 



while the "nett" phase satisfies 

A(a:)' = {[e- kl] l^oP + 20'o} + cos[A(a:)]. 

We can even substitute for 0(x) and thus rewrite this as a single integro- 
dijjerential equation for A(x): 

A(x)' = {[k'-kl] iV^oP +20'o} 

cos[A(x)]. (11.4.8) 



r - ki] 


2 


sinh ( ^ 


^Lik'-k'o] 1^0 


2 sin[A(x)] dx ) 
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This equation is completely equivalent to the original Schrodinger equation 
we started from. Unfortunately further manipulations seem intractable, and 
it does not appear practicable to push these observations any further. 

11.5 Application: Small shift in the potential 

Let us now consider the situation 

V{x) = Vo{x) + e6V{x), (11.5.1) 

for e "sufficiently small" . 

11.5.1 First-order changes 



To be consistent with previous notation (1.3.14) let us define 



e = kl + e S^^II^'^^kl + edv. (11.5.2) 

Using equation (11.2.43) we obtain the preliminary estimates 

\a{x)\ = l + 0{e^), (11.5.3) 

and similarly 

\b{x)\ = 0{e). (11.5.4) 

It is now useful to change variables by introducing some explicit phases so 
as to define 

a = a exp j [k^ - k^] \iPq\ dx^ ; (11.5.5) 

b = b exp (^-^ j [k^ - kl] lipH dx^ . (11.5.6) 



Doing so modifies the system of differential equations (11.2.26, 11.2.27) so 
that it becomes 

da i 

dx 

db i 
dx 



+ l[k' - kl] b exp I [e - kl] l^o'l dx) ; (11.5.7) 
-^[k^ - kl] a i^l exp j [k'^ - kl] {ipH dx^ . (11.5.8) 
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The advantage of doing this is that in the current situation we can now 
estimate 

0{e'), (11.5.9) 



da 
dx 



: 6v{x) ^l{x) exp [+i I e6v {ip'^l dx ] + 0{e'). (11.5.10) 



dx 

Integrating 

b{oo) = — ^ J 6v{x) ipl{x) exp J ^ I'^ol + ^(^^)- 

(11.5.11) 

This is not the standard Born approximation, though it can be viewed as an 
instance of the so-called "distorted wave Born approximation" [27] . In terms 
of the absolute values we definitely have 

|6(cx))| = |6(oo)| < - / \6v{x)\\^l{x)\dx + 0{e^). (11.5.12) 
11.5.2 Particle production 

When it comes to considering particle production we note that 

P = Po a{oo) + al b{oo) = po + b{oo) + C(e^), (11.5.13) 
so, since N = the change in the number of particles produced is 

SN = 6\(3^\ = Re{2(3*6(3} = Re{2al(3ob{oo)} + 0{e^). (11.5.14) 
In particular 

/ + 00 

\Sv{x)\ \ipl{x)\ dx + 0{e^). (11.5.15) 
■oo 

Since 

IttoPol = VT+W\Po\ = ViVo + 1 V^= v/iVo(iVo + l), (11.5.16) 
can also write as 

r+oo 

I 5Ar I < e v/iVo(iVoTT) / \6v{x)\\^l{x)\dx + 0{e^). (11.5.17) 



Note that one will only get an order e change in the particle production if the 
"known" problem {ipQ, ko} already results in nonzero particle production. 
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11.5.3 Transmission probability 

To see how a small shift in the potential affects the transmission probability 
we note 

1 _ 1 _ 1 

~W~ |«oa(oo)+/?o*&M|2 ~ |ao + /5o*&M + (^(e')r ^^^'^-^^^ 



But then 

1 



|aoP \l + {P*ob{oo)/ao} + 0{e^ 



implying 



So the change in the transmission probability is 



ttO 



Taking absolute values one obtains 

r>+00 



;il.5.19) 



T = To <; 1 - 2Re <; ^i^i^ K + O(e^) [> . (11.5.20) 



5T = -To <; 2Re ^i^M k + Qie') } . (11.5.21) 



/+00 
\6v{x)\ \tPl{x)\dx + 0{e^). (11.5.22) 
-oo 

Note that one will only get an order e change in the transmission probability 
if the "known" problem {tpQ, ko} already results in nonzero transmission (and 
nonzero reflection). 



11.6 Discussion 

We wish emphasize the advantages of the particular bounds derived in this 
chapter: 

• They are very simple to derive — the algebra is a lot less complicated 
than some of the other approaches that have been developed in earlier 
chapters of this thesis, and published in several papers [HHl EHl I901I9T] . 
(And a lot less complicated than some of the blind alleys we have 
explored.) 
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• Under suitable circumstances the procedure of this chapter yields both 
upper and lower bounds. Obtaining both upper and lower bounds is in 
general very difficult to do — see in particular the attempts in |89j . 

• All of the other bounds we have developed in earlier chapters of this 
thesis, and published in several papers [HHl EHl EOl EI] needed some 
condition on the phase of the wave-function, (some condition similar 
to v?' ^ 0), which had the ultimate effect of making it difficult to make 
statements about tunnelling "under the barrier". There is no such 
requirement in the present analysis. (The closest analogue is that we 
need ^ 0, which we normalize without loss of generality to = 1-) 
In particular this means that there should be no particular difficulty in 
applying the bound in the classically forbidden region — the "art" will 
lie in finding a suitable form for ^pQ which is simple enough to carry out 
exact computations while still providing useful information. 

In closing, we reiterate the fact that generic one-dimensional scattering prob- 
lems, which have been extensively studied for close to a century, nevertheless 
still lead to interesting features and novel results. 
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12.1 What we have achieved 

In this chapter we shall discuss the overall concept and achievements of this 
thesis. To a large extent much of the work has already has summarised 
at the end of each section. This thesis has been written with the goal of 
making it fully accessible to people with a basic background in non-relativistic 
quantum physics, especially in the physics of transmission, reflection, and 
Bogohubov coefficients. Mathematically, the key feature is an analytic study 
of the properties of certain second-order linear differential equations, and the 
derivation of analytic bounds on the growth of solutions of these equations (as 
a function of position and/or time). In this thesis we divided our efforts into 
analyzing four separate problems relating to rigorous bounds on transmission, 
reflection, and Bogohubov coefficients — these are considered under four 
separate themes: 

1. Bounding the Bogohubov coefficients, 

2. Bounding the greybody factors for Schwarzschild black holes, 

3. Transmission probabilities and the Miller-Good transformation, and 

4. Analytic bounds on transmission probabilities. 
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In addition, all four of these separate themes which are reported in this 
thesis, are also the seeds for the various published journal articles (plus one 
submitted article) that have already arisen from this thesis. 

This thesis is divided into twelve main chapters. In the following, we shall 
summarise and anaylse the main work we derived in each chapter. 

12.2 The main analysis: 

Structure of the thesis 

We first provided sufficient context for the reader to appreciate the role played 
by the various topics to be discussed in this thesis, and to place them into a 
wider perspective. Firstly, we introduced the Schrodinger equation — a spe- 
cific partial differential equation used in the development of the "new" (1925) 
quantum theory. Secondly, we provided the basic theory underlying the 
WKB approximation., and the concept of the time-independent Schrodinger 
equation — both are foundations for all our subsequent analyses. We have 
presented a very general introduction to these concepts first — so that the 
bounds we derived on transmission, refiection, and Bogoliubov coefficients 
were easier to understand. 

We mainly considered the scattering theory in one space dimension — 
because it is mathematically simple and physically transparent. 

In particular, it is interesting to show how to derive the basic ideas of 
transmission and refiection directly by using scattering theory. In addition, 
we have just seen an important connection between refiection and transmis- 
sion amplitudes. We called the probability that a given incident particle 
is refiected as the "refiection coefficient". While the probability that it is 
transmitted is called the "transmission coefficient" . 

In chapter [3} we collected many known analytic results in a form amenable 
to comparison with the general results we subsequently derived. In addition, 
we also introduced the concept of quasinormal modes [QNM] . We used these 
tools for comparing the bounds with known analytic results. Moreover, we 
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reproduced some of the analytically known results, and showed (or at least 
sketched) how to derive their scattering amplitudes, and so calculate quanti- 
ties such as the tunnelling probabilities and quasinormal modes. We did this 
explicitly for the delta-function potential, double-delta-function potential, 
square potential barrier, tanh potential, sech^ potential, asymmetric square- 
well potential, the Poeschl-Teller potential and its variants, and finally the 
Eckart-Rosen-Morse-Poeschl-Teller potential. 

We also obtained a number of significant bounds, considerably stronger 
than those in [88], of both theoretical and practical interest. Even though 
the calculations we have presented are sometimes somewhat tedious, we feel 
however, they are more than worth the effort — since there is a fundamen- 
tal lesson to be learnt from them. Technically, we demonstrated that the 
Schrodinger equation can be written as a Shabat-Zakharov or a Zakharov- 
Shabat system, which can then be re-written in 2 x 2 matrix form. 

In chapter [5| we have again moved our attention back to a Shabat- 
Zakharov system of ODEs by re-casting and describing the first derivation 
of scattering bounds as presented by Visser in reference |BB]- The formal- 
ism as developed here works in terms of one free function if{x). In other 
parts of this thesis we have established generalized bounds; some in terms of 
two arbitrary functions ip{x) and xi^)^ ^"^^ some in terms of three arbitrary 
functions ip{x), A(x), and x{^)- The derivation of this chapter is noteworthy 
because of its brevity and simplicity. 

All of the above techniques from chapter [5] are important to develop a 
number of interesting bounds in chapter |6j We dealt with some specific cases 
of these bounds and develop a number of interesting specializations. We 
have collected together a large number of results that otherwise appear quite 
unrelated, including reflection above and below the barrier. In addition, we 
have divided the special case bounds we considered into five special cases: 
special cases 1-4, and "future directions". Finally, we took further specific 
cases of these bounds and related results to reproduce many analytically 
known results. 
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Consequently, we have re-cast and represented these bounds (from chap- 
ter |6| in terms of the mathematical structure of parametric oscillations. This 
time-dependent problem is closely related to the spatial properties of the 
time-independent Schrodinger equation. 

In chapter [8| we re-assessed the general bounds on the Bogoliubov co- 
efficients developed in earlier chapters of this thesis, and published in ref- 
erence [SB], providing a new and largely independent derivation of the key 
results, one that short-circuits much of the technical discussion in chapter [5} 



After this investigation about bounding the Bogoliubov coefficients and 
their techniques we have moved to study the greybody factors in Schwarzs- 
child black hole. The "greybody factor" is actually a synonym for "transmis- 
sion probability" . Indeed, the phrase "greybody factor" is used more in the 
thermodynamics and spectroscopy communities, while the phrase "transmis- 
sion probability" is used more in the quantum mechanics community, but 
they are referring to the same concept. In this thesis, we developed a com- 
plementary set of results — we derived several rigorous analytic bounds that 
can be placed on the greybody factors. Even though these bounds are not 
necessarily tight bounds on the exact greybody factors, they do serve to focus 
attention on general and robust features of these greybody factors. Moreover 
they provide a new method of extracting physical information. Furthermore, 
we considered the greybody factors in black hole physics, which modify the 
naive Planckian spectrum that is predicted for Hawking radiation when work- 
ing in the limit of geometrical optics. 

We used the Miller-Good transformation (which maps an initial Schrod- 
inger equation to a final Schrodinger equation for a different potential) to 
significantly generalize the previous bound. Moreover, we shall see that the 
Miller-Good transformation is an efficient process to generalize the bound, 
to make it more efficient and powerful. 

Finally, we have again shifted our attention back to the analytic bounds 
and transmission probabilities context. We developed a new set of tech- 



and published in reference 
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niques that are more amenable to the development of both upper and lower 
bounds. Moreover, we derived significantly different results (a number of 
rigorous bounds on transmission probabilities for one dimensional scattering 
problems), of both theoretical and practical interest. Instead of explaining 

Several ways to derive bounds for arbitrary wave phenomena 
(including Schwarzchild black hole greybody factors) 



Method 


Characteristics and properties of bounds 


[[)). btanclarcl Wlvr) 


— uncontrolled approximation. 

- do not know if the approximation 
is high or low. 


(1). WKB "basis" 


- rigorous bounds. 

- use "gauge fix" method. 


(2). Bogoliubov coefficients 


- rigorous bounds. 

- do not use "gauge fix" . 

- short-circuits the technical details of method (1). 


(3). general bounds on the 
greybody factors 


- obtained several rigorous analytic bounds 
that can be placed on the greybody factors. 


(4). the Miller-Good transformation 


- generally not "WKB-like" . 

- no need to separate the region into 
classically allowed and forbidden regions. 


(5). the analytic bounds 


- very simple to derive the bounds. 

- procedure yields both upper and lower bounds. 



Table 12.1: This table shows several ways to derive bounds for arbitrary wave 
phenomena (including the greybody factors of Schwarzchild black holes), and 
also the key properties of these bounds. (See chapters [sj [9| 10, and 11 for 
more details.) 



the details of the analysis yet again, 

we would like to stress a few points that we believe are useful to under- 
stand the overall concept of the thesis: 
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• The Schrodinger equation describes the space -and time- dependence of 
quantum mechanical systems, and its apphcation to the wave function 
are the basic idea that describes the wavehke properties of a subatomic 
system. 

• The probabihty current express the reflection and transmission coef- 
ficients. The probability current is based on the assumption that the 
intensity of a beam is the product of the speed of its particles and their 
linear number density. It is then a mathematical theorem that this 
probability current is conserved. 

• The WKB approximation is generally applicable to problems of wave 
propagation in which the frequency of the wave is very high or equiva- 
lently, the wavelength of the wave is very short. 

• The ideas of reflection and transmission of waves in both unbound and 
bound states are important. By considering reflection and transmission 
of waves in unbound states, we have seen that in principle they are 
completely specified by the potential function V{x). 

• The quasinormal modes [QNM] are the modes of energy dissipation 
of a perturbed object or field. In particular, the most outstanding 
and well-known example is the perturbation of a wine glass with a 
knife: the glass begins to ring, it rings with a set, or superposition, 
of its natural frequencies - its modes of sonic energy dissipation. As 
previously explained, when the glass went on ringing forever, we can 
call these modes normal. For instance, here the amplitude of oscillation 
decays in time, so we call its modes quasi-normal [H] . 

• The Schrodinger equation can be written as a Shabat-Zakharov system, 
which can then be re-written in 2 x 2 matrix form. We rearranged 
this formation in terms of a generalized position-dependent "transfer 
matrix" involving the symbol V which denotes "path ordering" . 
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• A "parametric oscillator" is a simple harmonic oscillator whose param- 
eters (its resonance frequency uj and damping time (3) vary in time. The 
other interesting way of understanding a parametric oscillator is that 
it is a device that oscillates when one of its "parameters" (a physical 
entity, like capacitance) is changed. 



12.3 Further interesting issues 

There are some interesting ways this thesis could be extended in future work. 
We would like to wrap-up by providing a list of things that we believe are 
interesting to continue to analyze: 

• This present research certainly deserves more work on how to extend 
the bounds in many different ways. While we have already established 
several powerful techniques to derive rigorous bounds on transmission, 
reflection, and Bogoliubov coefficients, we feel, however, that there are 
probably "optimal" bounds still waiting to be discovered. 

• In particular, it is apparent that the current bounds in chapter [8] are 
not the best that can be achieved, and we strongly suspect that it may 
be possible to develop yet further extensions to the current formalism. 
It is in fact possible that the "more general" bounds are close to being 
discovered and will have further development in the near future. 



The bounds presented in chapter 10 are generally not "WKB-like 



apart from the one case reported in equation (10.5.15) there is no need 
(nor does it seem useful) to separate the region of integration into 
classically allowed and classically forbidden regions. In fact it is far 
from clear how closely these bounds might ultimately be related to 
WKB estimates of the transmission probabilities, and this is an issue 
to which we hope to return in the future. 

Finally, we have seen that even though the topic considered in this 
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thesis is ultimately a quantum mechanics subject, dating back to 1925, 
this does not mean that everything has already been done. It is con- 
ceivable that the new techniques in this project help us to derive more 
rigorous and tighter bounds for the barrier penetration probability. 

All the above suggestions would be interesting and feasible, although some 
of them would be more tedious to work on than others. 

In summary, this thesis provides a platform for better understanding the 
rigorous bounds that one can place on the Bogoliubov coefficients associated 
with a time-dependent potential, and the several rigorous analytic bounds 
that can be placed on the greybody factors. This thesis developed a way of 
looking for nice and accurate bounds. Furthermore, one primary goal of this 
thesis was to explore the best way of finding barrier penetration probability. 
In conclusion, we can say that this project will be another step to improv- 
ing our understanding of quantum mechanics, in particular, non-relativistic 
quantum physics, especially in regard to transmission, reflection, and Bogoli- 
ubov coefficients. 
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